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1 Introduction 

Up to the first half of fifties of 20th century all known elementary particles were assumed to 
be structureless, i.e. point-like. The latter property is reflected also in the principles of local 
quantum field theory (QFT) (unifying consistently quantum theory, the concept of the field and 
the relativistic invariance), which is considered to be a dynamical theory of elementary particle 
mutual interactions. 

Only Hofstadter experiments [1] on elastic scattering of electrons on protons at SLAC have 
clearly demonstrated disagreement between theoretical expression for the cross-section calcu- 
lated in the framework of the quantum electro-dynamics (QED) and the obtained experimental 
results. This phenomenon have revealed the non-point-like nature of the proton, which later on 
have been extended also to all other existing strongly interacting particles. 

As a result at the calculation of the matrix element of the elastic scattering of electrons on 
hadrons in one-photon-exchange approximation one does not know explicitly (unlike the elec- 
tron) the electromagnetic current of considered hadron to be written only symbolically and prac- 
tically using various symmetries it is decomposed into maximal number of linearly independent 
co-variants constructed from momenta of incoming and outgoing hadrons and their spin parame- 
ters. The corresponding coefficients are scalar functions (the electromagnetic (EM) form factors 
(FFs)) of one variable to be chosen the squared momentum t = q 2 = —Q 2 transferred by the 
virtual photon. Their number depends on the spin of the considered hadron. 

Similarly can be introduced the weak FFs of hadrons, representing the contribution of the 
weak structure of hadrons into the dynamical quantities describing a weak interaction of hadrons 
in various weak processes. 

A natural explanation of the electro-weak structure of hadrons was obtained only after the 
discovery of quark-gluon structure of strongly interacting observable particles. 

The behavior of the EM FFs in the whole interval of their definition from — oo to +oo is 
expected to be theoretically predicted by the quantum chromo-dynamics (QCD), the gauge in- 
variant local QFT describing mutual interactions of colored quarks and gluons. However, as it 
is well known, merely at sufficiently small distances (thanks to its asymptotic freedom) QCD 
becomes a weakly coupled quark-gluon theory to be amenable to a perturbative expansion in 
the running coupling constant a s and predicts just the asymptotic behavior of FFs. In low mo- 
mentum transfer region, where a s becomes large, the quark-gluon perturbation theory breaks 
down and non-perturbative methods in QCD are not well worked out to give interesting results 
on the FFs of hadrons. The same is valid also for the low energy time-like region where FFs are 
complex functions of their variable and acquire the most complicated, resonant, behavior. 

In such situation a phenomenological approach based on the analyticity of FFs starts to be 
very useful, which compensates above-mentioned problems to some extent and renders possible 
to achieve a line of interesting results. 

In an interpretation of experimental data on EM FFs appears to be useful utilization of the 
analytic properties in the form of integral (so-called dispersion) relations together with the uni- 
tarity condition of FFs, which have brought the investigated FFs into relations with other FFs 
and amplitudes of various processes of strongly interacting particles. Such approach in the case 
of nucleons have led to a prediction of an existence of isoscalar and isovector vector mesons, 
and subsequently to the vector- meson-dominance (VMD) model [2]. The latter model is based 
on the assumption (to be later on experimentally confirmed in electron-positron annihilation into 
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hadrons), that the interaction of the virtual photon with hadron is realized by a transformation 
of the photon to a vector meson with the same quantum numbers and then this vector meson is 
strongly interacting with considered hadron as in any other hadron collision. 

The VMD model was revealed before the discovery of the quark model of hadrons. Despite 
of this fact the latter is consistent with VMD model. Really, at the energy of the photon nearly to 
the mass of the vector meson the latter is changed to the quark-antiquark pair, which is as a result 
of the confinement effect immediately bound into the vector meson with the photon quantum 
numbers. 

Though the VMD model from the point of view of the global analysis of existing FF data has 
been in the past very frequently applied, it suffers from a lot of shortcomings. It does not take 
into account the instability of vector mesons, the unitarity condition of FFs and also the analytic 
properties of FFs, which could lead to a more realistic behavior of FFs in the time-like resonant 
region. Other serious shortcoming is the same asymptotic behavior for FFs of all hadrons, which 
is in contradiction with the predictions of quark model for baryons and atomic nuclei. 

A solution from this situation is the universal Unitary and Analytic (U &cA) model of electro- 
weak FFs, which is a unification of the experimental fact of a creation of unstable vector-meson 
resonances in the electron-positron annihilation into hadrons, the analytic properties of FFs in 
the complex t-plane and the correct asymptotic behavior of FFs as predicted by the quark model 
of hadrons. Its applications to many hadrons and nuclei have led to a lot of interesting results to 
be verifiable experimentally. 
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2 Electromagnetic form factors of strongly interacting particles and their properties 

In this section general properties of EM FFs of strongly interacting particles, like number of FFs 
of the hadron under consideration, their most usual parametrization and its shortcomings, will be 
reviewed and subsequently resolved. 



2.1 Number of electromagnetic form factors of a considered hadron 

The elastic electron scattering e~h — > e~h and the annihilations e + e~ «-> hh, where h means 
an arbitrary strongly interacting object (including also the atomic nucleus), are the most usual 
processes, in which the concept of EM FF appears. The cross-sections of these processes are 
proportional to the absolute value squared of the corresponding scattering amplitudes, which 
one knows formally to calculate in the framework of the quantum electrodynamics (QED) per- 
turbation expansion according to the fine structure constant a ~ 1/137. Since the value of 
the fine structure constant a<C 1, those scattering amplitudes are considered practically in the 
one-photon-exchange approximation as follows 

M(e~h -> e~h) ~ M™(«,t) = e 2 S(fc 2 ) 7 ^(fc 1 )^(^| J? M \h) (2.1) 

and 

M(e+e- -> hh) ~ M^(t,s) = e 2 i;(fc 2 )7 Al u(fc 1 )^(0| J? M \hh), (2.2) 

where g^/q 2 is the photon propagator and (h\J^ M \h) (resp.(0| M \hh)) is a matrix element 
of the hadron EM current, which, however, due to the non-point-like nature of the hadron h, is 
unknown. Therefore, in practice it is decomposed according to a maximal number of linearly 
independent relativistic covariants constructed from the four-momenta and spin parameters of h 
as follows 

{h\J* M \h) =Y,Vi{t) (2-3) 

i 

or 

(M|J / f M |0)=^X^ i (t), (2.4) 

i 

where the scalar coefficients Fj (t) are the EM FFs of the hadron h as functions of one invariant 
variable t- the momentum squared to be transferred by a virtual photon. 

The number of Fi (t) depends essentially on the spin S of h. 

Let us consider the most topical cases. 

We start with the nonet of pseudoscalar mesons 7r + , tt°, ir~, K + , K°, K°, K~, 77, 77'. Since 
they possess spin to be zero, for the construction of covariants (p 2 — and (p 2 + p\)n only 
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two four-momenta, p\ and pi, are available. By an application of the gauge invariance of the EM 
interactions one comes to the following final parametrization 

(P2|J m (0)|pi) = F P (t)( Pl + P2 )^ (2.5) 

only with one EM FF Fp(t) completely describing the EM structure of any member of the nonet 
of pseudoscalar mesons. Moreover, making use of transformation properties of the EM current 
operator Ju{x) and the one-particle state vector with regard to the all three discrete C, P, T 
transformations simultaneously, one finds the relation between particle and antiparticle FFs 

F P (t) = -Fp(t) (2.6) 

from where it follows that for the true neutral pseudoscalar mesons 7r°, rj, rj the EM FFs are 
identically equal zero for all values from — oo < t < +oo. 

The pseudoscalar meson EM FFs are normalized at t = to the charge of the considered 
meson. 

A consideration of the nonzero value of the isospin of the pion does not enlarge a number of 
EM FFs and both charged pions are described by the same FF. 

A completely different situation is with kaons. The K + and K° belong to the same iso- 
multiplet with 1 = 1/2. Therefore instead of the positively charged and neutral kaons one can 
introduce the EM current of the kaon and to investigate what isotopic structure it has. One can 
show, that it splits on a sum of isotopic scalar and isotopic vector. In connection with the latter 
the isoscalar F^ (t) and isovector F'^ (t) FFs of the kaon are introduced to be expressed by F£ (t) 
and (t) as follows 

F s K (t)= 1 -[F+(t) + F° K (t)] (2.7) 



F v K {t) = \[F+{t)-F%{t)], (2.8) 
from which immediately the normalization condition 

Fk(0) = Fk(0) = \ (2-9) 
is obtained. 

In principle, there is no problem of obtaining of the experimental information on |i*p(i)| in 
t < and t > regions as da/d^l — \F P (t)\ 2 . However, the data on nuclei (e.g He 4 , C 12 , O 16 ) 
exist only for t < up to now and there is no concept of the EM FF of nucleus for t > to be 
known by nuclear physicists. 

As a consequence of a compound nature of nuclei so-called diffraction minima appear in 
t < region at the absolute value of their charge FF |-F C (*)| > which are interpreted as zeros of 
F c (t) on the real axis of the complex t- plane. It is observed, that if more compound nucleus is 
investigated more diffraction minima emerge in the same range of momentum transfer values. 

In the case of the octet 1/2+ baryons p, n, A, E + , S°, S°, S~ and e.g. He 3 , H 3 nuclei 
covariants R fl (pi,p2) are constructed by the four-momenta pi,P2, Dirac matrices and bispinors. 
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The final result for a parametrization of the matrix element of the EM current of the octet l/2 + 
baryons takes the form 

(p 2 \J tl {Q)\ Pl ) = ^u(p 2 ){^F 1B (t) + J—a^fa -p 1 ) u F 2B (t)}u{p 1 ) 1 (2.10) 

where F\ B {t) and F 2B (t) are Dirac and Pauli FFs, respectively and m B is the baryon mass. 
From the practical point of view it is more suitable to describe the EM structure of the octet 
1 /2+ baryons by means of the Sachs electric Geb (t) and magnetic Gmb (t) FFs, defined by the 
following expressions 

G E B(t) = F 1B (t) + -^F 2B (t); G MB (t)=F 1B (t) + F 2B {t). (2.11) 

There is a special coordinate system (the Breit reference frame), in which GEB{t) and Gmb(() 
describe a distribution of the charge and magnetic moment of the baryon. Hence they are called 
the electric and magnetic FFs to be normalized to the charge and magnetic moment of the baryon, 
respectively, for t = 0. 

Similarly to kaons, one can consider instead of the EM current of every member of the octet 
l/2 + baryons the EM current of the corresponding isomultiplets and to look for their splitting 
into isoscalar and isovector parts. As a result one finds the following decomposition of the 
nucleon and A—, E— and S— hyperon electric and magnetic FFs into isoscalar and isovector 
parts of the Dirac and Pauli FFs 

G Ep (t) = [FlAt)+F?At)] + ^[Fi N (t) + FZ N (t)} (2.12) 



Gu P (t) = [F s 1N (t) + F? N (t)] + [F 2 n (t) + FZ N (t) 

t 

4m;_ 

G M n(t) = [Fl N (t)-F? N (t)] + [Fi N (t)-FZ N (t)} 



G E „W = [Fl N (t)-F7 N (t)] + ^[Fi N (t)-F^ N (t)} (2.13) 



GeaW = F s lK (t) + ^F s 2K (t) 

Gma(*) - F! A (t) + F* A (t) (2.14) 



G ES+ (t) = (t) + F^ (t)} + i - 3 - [F| s (t) + F^ (t)} 

Gme+ (*) = (*) + ^ (*)] + (*) + F^ (t)} 

Get, (t) = F^(t) + ^-F^(t) 

G ME o(t) - F^(t) + F| s (t) (2.15) 

Gee-(*) = [F^-FM^ + ^iF^-FMt)} 

Gut- (*) = (t) - (t)] + [F| s (t) - FJ E (t)} 
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Gee" (t) = [F! 3 (t)+F^(t)] + ^-[F! s (t)+F^t)] 

G M ~°(t) = [F s 13 (t)+F^(t)} + [F^(t) + F^(t)} (2.16) 

t 

Ami 

Cue- (t) = [F? E (t) - F? E (t)} + [F| - F^ (t)} . 



The experimental information on GEB(t), Gmb^) in i < region can be easily determined 
from parameters of the straight-line (so-called Rosenbluth) plot of 

da(e-B ->e-B) f a 2 cos 2 (tf/2) 1 

_ ' \4i? 2 sin 4 (??/2)[l + (2E/m B )sm 2 (-d/2)] J " 
= + 5 (i) tan 2 ((9/2), (2.17) 

where 



GebW ~ 4^G 2 MB (t) 
1 *_ " ~~ ' 

t 



B W = - 2 ^ G mbW (2-18) 



4m 



2/ 



in the laboratory system versus tan 2 ($/2) at fixed t and for nuclei again diffraction minima 
appear. 

Till now all existing data on Geb (t), GMs(*)i n * > region are obtained from cr tot (e + e~ o 
BB) under the assumption that \GEB{t)\ = \GMB{t)\- 

The covariants i? M (pi,p 2 ) for EM FFs of the nonet of vector mesons, p+, p°, p~, if*+, 
if* , A"* , if* - , ui, $ and also of the deuteron are constructed by the four-momenta p\, p 2 
and polarization vectors. Then a parametrization of the matrix element of the EM current of 
vector-particles takes the following form 

( P 2\J* M \pi) = ■ZK+Fzm.ie - g )-ff (g. g )]- F 3 (t) (C ' ^ ' g) d„, (2.19) 

where £ and £' are polarization vectors for incoming and outgoing particles of four-momenta p\ 
and p 2 , respectively 

C'-P2 = 0; e- P i = 0; C' 2 = -i; C 2 = -i; 

d M = (P2+Pl) M ; Qn = (P2-Pl)p- 

Practically, it is convenient to describe the EM structure of vector particles by an analogue of 
the Sachs FFs of nucleons 



G c (t) = Fl (t) - -^G Q (t); G M (t) = F 2 (t); 
bm 



H 

t 



G Q {t) = Fl (t)-F 2 (t) + (l-— 2 -)F 3 (t) (2.20) 
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the names of which, the charge Gc{t), the magnetic Gm(£) and the quadrupole Gg(t) FFs, are 
derived from the fact that their static values correspond to the charge, magnetic and quadrupole 
moment of the vector particles. 

One can determine all Gc(t), Gm(0> Gq(t) FFs from da/dfl of the e~V — > e~V process 
provided that polarized particles are used in the corresponding experiments. Otherwise only the 
elastic structure functions A(t) and B(t) can be drawn out from 



da 

dn 



(2.21) 



where 



A(t) = - 
B{t) = 



4 (1 ~4 )Gyt)+G ^ )+ 

3m^> Am 2 



t 



fG 2 M (t)- 



On the other hand, from 

<Jtot{e + e- VV) = 



3t 



3|G c (t)| 2 + 



\GM(t)\ 2 + -^- T \G Q (t)\ 



6my 



(2.22) 



(2.23) 



one can see immediately that it is not a single task to obtain any experimental information on the 
corresponding EM FFs in t > region. 

For strongly interacting particles h with S > 1 a situation is even more complicated and 
generally it is not solved up to now. 



2.2 Properties of electromagnetic form factors of hadrons 

Summarizing our knowledge about the experimental behavior of EM FFs we come to a conclu- 
sion, that all of them have a similar behavior in the shape. But they differ in the asymptotic 
behavior, normalization, number of bumps corresponding to vector-meson resonances and also 
in the shape and height of those bumps. 

A behavior of EM FFs is a matter of predictions of a strong interaction dynamical theory. 
However, there is no such theory able to predict a correct behavior of \Fh (t) | for — oo < t < +oo 
up to now and only partial successes were achieved in this direction. 

The great discovery in the elementary particle physics was a revelation of the quark-gluon 
structure of hadrons and its direct relation [3, 4] to the asymptotic behavior of EM FFs to be 
determined by a number of constituent quarks n q of the hadron h as follows 

*fc(*)|t|-oo~* 1 " ne , (2-24) 

which is in a qualitative agreement with existing experimental data. 

On the other hand, it is well known that on the role of a true dynamical theory of strong 
interactions QCD, the gauge-invariant local quantum field theory of interactions of quarks and 
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gluons, is pretending. But, as a consequence of the asymptotic freedom of QCD, in the frame- 
work of the perturbation theory, the latter is able to reproduce [5-7] just the asymptotic behavior 
(2.24) up to logarithmic corrections. 

Not even the nonperturbative QCD sum rules [8] by means of which a prediction [9, 10] of a 
behavior of EM FFs in a restricted t < region is achieved, solve the problem of a reconstruction 
of EM FFs in the framework of QCD completely. 

For a completeness we mention also the chiral perturbation approach [11], in the framework 
of which a correct behavior of EM FFs of hadrons around the point t = is predicted. This 
is very important to be mentioned as the chiral perturbation approach is equivalent to QCD at 
low energies where the running coupling constant a s (t) takes large values and the PQCD is 
nonapplicable. 

Summarizing, QCD (not even its equivalent form) gives no quantitative predictions in the 
most important part of the time-like (4m^ < t < 4GeV 2 ) region, where EM FFs are already 
complex functions of t and the e + e~ annihilation experiments exhibit a nontrivial behavior of 
measured cross-sections caused by a creation of various unstable vector-meson states. 

Therefore for the present an appropriate phenomenological approach based on the synthesis 
of the experimental fact of a creation of vector-mesons in e + e~- annihilation processes into 
hadrons, the asymptotic behavior (2.24) and the well-established analytic properties, leading to 
the U&A model of EM structure of strongly interacting particles, is still the most successful way 
in a global theoretical reconstruction of EM FFs of hadrons. 

The vector-meson creation in e + e~- annihilation processes into hadrons is taken into account 
by means of the vector-meson-dominance (VMD) model given for isoscalar and isovector parts 
of EM FFs by the relation 



where f V hh an ^ fv are tne vector-meson-hadron and the universal vector-meson coupling con- 
stants, respectively, and my is the vector- meson mass. 

The analyticity consists in a hypothesis that all EM FFs are analytic functions in the whole 
complex t- plane besides infinite number of branch points on the positive real axis corresponding 
to normal and anomalous thresholds. 



There are experimentally confirmed neutral vector- mesons [12] with quantum numbers to be 
identical with photon. They have the isospin either or 1 . 

On the other hand the EM current of hadrons is by a rotation in the isospin-space trans- 
formed like the sum of isotopic scalar and the third component of isotopic vector. The latter 
transformation properties reflect well known fact of the non-conservation of the isospin in the 
EM interactions and lead automatically to a phenomenon, observed experimentally, that at the 
absorption and creation of virtual photon by hadron the isospin value can be changed by 1. 
Therefore the photon can be considered to be a superposition of states with isospin value (the 
isoscalar photon) and isospin value 1 (the isovector photon). Then neutral vector-mesons with 
quantum numbers J PC = 1 differ from the photon only by the mass and there is no obstacle 




(2.25) 



2.3 Vector-meson-dominance model for form factors of hadrons 



12 



Analyticity in a phenomenology of electro-weak structure of hadrons 




t 

Figure 2.1. Decay of neutral vector meson into lepton pair. 

h 






V 



h 

fvhh 



fv 



r 



Figure 2.2. Approximation of hadron FF by a sum of VMD terms. 

to assume that there are between virtual photons and vector-mesons transitions with a definite 
probability determined by the coupling constant 



9-yV 



'fv 



(2.26) 



where e is the electric charge , my is the mass of vector-meson and fy is the so-called universal 
vector-meson coupling constant. 

The transition of virtual photons to neutral vector mesons is confirmed practically by the 
neutral vector-meson decay into lepton-antilepton pair. Really, if we assume the transition 7 O 
V°, then the V° — > £ + i~ decay can be explained as a result of a change of V° into virtual photon 
7* which subsequently is creating the lepton-antilepton pair (Fig. 2.1). 

It seems to be natural a generalization of this mechanism to any process of an interaction of 
photon with hadron, in which first the photon is changed to vector-meson and then the latter is 
interacting with the hadron like in other hadron collisions by strong interactions. 

In conformity with the idea of VMD model any EM FF of hadron in the first approximation 
can be represented by a sum of Feynman diagrams with and exchange of vector-mesons V° 
(Fig. 2.2). 

By an application of standard methods of quantum field theory to an explicit calculation of 
contributions of this sum of Feynman diagrams one obtains 



{9-yV-fvhh) 



(2.27) 



the VMD parametrization (in the zero width, i.e. Ty = approximation) of the hadron EM FF, 
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which is, taking into account also the relation (2.26), normalized for t — in the form 

F h (0) = ^(W/v) (2.28) 
v 

and it possess the asymptotic behavior 

to be the same for EM FFs of all strongly interacting particles. 

The relation (2.26) can be derived by taking, first the expression (2.27) for the pion FF, the 
isovector part of the kaon FF and the isovector part of the Dirac nucleon FF to be saturated 
by only the lowest vector-meson, the p-meson. As a result one can write the following three 
independent expressions 



eF^t) = g7p / p "" (2.30) 



ml 



eF v K (t) = 9 ^ 2 fpK t K (2.31) 



in which the normalization condition at t = gives 

(2.33) 



m 2 p 



2 ml 



9_2pJ%n (235) 



e 

2 = ~m 2 p 



By a multiplication of the last two equations by 2, one can obtain the following universal inter- 
action of the p-meson to be expressed by the relations 

fpTTTT = 2.fpKK = ^fpNN = /p' (2.36) 

where the coupling constant f p is named to be the universal coupling constant of the p-meson 
and the coupling constant of the p-meson with the virtual photon g Jp can be written as follows 

m 2 

9 1P = e-^. (2.37) 
Jp 

Generalization of the previous relation to any other vector meson V with quantum numbers 
of the photon leads to the form of (2.26). 
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2.4 Asymptotic conditions for form factors represented by VMD model 

As we have noticed above, from one side the quark model of hadrons predicts (2.24) the asymp- 
totic behavior of EM FF of hadron to be dependent, in conformity with existing experimental 
data, on the number of constituent quarks of considered hadron. 

On the other hand, the VMD model for EM FFs (2.3) of strongly interacting particles gives 
the same asymptotic behavior (2.29) independently on the number of constituent quarks. 

This conflict is solved in this paragraph [13] by a derivation of the so-called asymptotic 
conditions. However, there are known in practice two seemingly different asymptotic conditions. 
Further we clearly demonstrate their equivalence. 

Generally, let us assume that the FF in (2.25) is saturated by n different vector meson pole 
terms and it has to have the asymptotic behavior 

F h (t)\ tHoo ~t- m , (2.38) 

where m < n. 

Transforming the VMD pole representation (2.25) into a common denominator one obtains 
FF in the form of a rational function with a polynomial of (m — 1) degree 

P n _i(i) = A + A l ■ t + A 2 ■ t 2 + ■ ■ ■ + A n _! • i"" 1 (2.39) 

in the numerator, where 

n 



in-2 — 



(-i)"- a E^ a E m J°i 



<=i 3 = 1 



An-3 = (-1)" 3 Yl E'"-"' 

«1.«2 = 1 j = l 

n n 
An-4 = (-IT- 4 E m l m l m lY, m hl 

i l< i 2< i 3> i rJ^3 



*n— (ra— 1) 



^_^n—m+l 
^ -yyi—m 



E 



■ rn 



H >*2>"-*m-2 = A 
H<*2- ■-< i m-2' i r^j 



. E "'""■> 

3 = 1 



(2.40) 



E 



m 2 m 2 2 



,m L-2 TO L 



•l>*2>"-*m-l =1 
*l<*2'"<*m — l^r^S 



-iE m '^ 
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A 2 = (- 1 ) 2 m n m l"' m l3E m i fl J 



*l>*2'"'*n-3 — 1 j — 1 

il<i2--<i n -3> i r¥ : 3 

71 n 

M = (-1) m n m l'" m l3 m LE m J fl J 

*li«2>"-*n-2 = 1 j — 1 

*l<*2"-<*n — 2>*r^ 

n n 

i l > i 2r--i n — 1=1 j = l 

il<*2""<*n-l>*r7 £ 3 

anddj = {fjhh/fj)- 

In order to achieve the assumed asymptotic behavior (2.38) one requires in (2.39) the first 
(to — 1) coefficients from the highest powers of t to be zero and as a result the following first 
system of linear homogeneous algebraic equations for the coupling constant ratios is obtained 

n 

E m ) a o = 

3=1 

E m i z2 m2 j a 3 = 

E m?m? ZX°j = ( 2 - 41 > 

n,'2=i ■_ 1 

E m l m l m l22 m2a i = 

t l,»2,» 3 = l ._. 
•l< l 2<*3.»r^J J 



n ^ 

E ™\ m \ ■ ■ ■ m L-2 z2 m2 3 a o = °- 

«l.«2.'--»m-2-l - = 1 

il<i2-'-<*m-2.*r^3 J 

As one can see from (2.41) with increased to the coefficients become sums of more and more 
complicated products of squared vector-meson masses. 

For a derivation of the second system we employ the assumed analytic properties of EM FFs 
of hadrons, consisting of infinite number of branch points on the positive real axis, i.e. cuts. The 
first cut extends from the lowest branch point to to +oo. Then one can apply the Cauchy theorem 
to FF in t- plane 

±-.j> F h {t)dt = (2.42) 

where the closed integration path consists of the circle Cr of the radius R — > oo and the path 
avoiding the cut on the positive real axis. As a result (2.42) can be rewritten into a sum of the 
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following four integrals 

-M f F h (t)dt+ [° F h (t-ie)dt+ f F h (t)dt + f F h (t + ie)dt \ =0 (2.43) 

2m [Jc R J+oo Jc r /2 J t J 

where e <C 1 and C r /2 is the half-circle joining the upper boundary of the cut with the lower- 
boundary of the cut around the lowest branch point t . The contribution of the first integral in 
(2.43) is zero as Fh (t) for R — > oo is vanishing. One can prove also that the third integral in 
(2.43) for r — > is zero. As a result one gets 

1 r°° 

- [F h (t + ie)-F h (t-ic)]dt = 0. (2.44) 

J t 

Then, taking into account the reality condition of FF 

FZ(t)=F h (t*) (2.45) 

following from the general Schwarz reflection principle in the theory of analytic functions, one 
arrives at the integral superconvergent sum rule 

i r 00 

- / ImF h (t)dt = (2.46) 

n J to 

for the imaginary part of the FF under consideration. 

Repeating the same procedure for the functions tFh(t), t 2 Fh(t), • • • , t m ^ 2 Fh(t) which pos- 
sess the same analytic properties in the complex i-plane as F h (t), one gets another (m — 2) 
superconvergent sum rules 

i r 00 

- / t-ImF h (t)dt = 

* J to 

1 f 00 

- / t 2 • ImF h (t)dt = (2.47) 
n J to 



-I 





OO 



t m - 2 ■ ImF h (t)dt = 0. 

'u, 

Now, approximating the FF imaginary part by 6- function in the following form 

n 

ImF(t) = 7T aiS(t - m\)m\ (2.48) 

i 

and substituting it into (2.46) and (2.47) one obtains the second system of (m — 1) linear homo- 
geneous algebraic equations for coupling constant ratios a, = (fihh/fi) 



^2 m t a * = ° 

n 
n 

^mfa l= (2.49) 



i=l 
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E2(m-2) „ 

i=i 

n 

E2(m-1) „ 

i=l 

where coefficients are simply even powers of the vector-meson masses. 

Further we demonstrate explicitly that both systems of the algebraic equations, (2.41) and 
(2.49), are equivalent, despite the fact that they have been derived starting from different proper- 
ties of the EM FF, and thus they appear to be different. 

We start with the equations (2.41). From a direct comparison of systems (2.41) and (2.49) 
one can see immediately the identity of the first equations in them. 

The second equation in (2.41) can be written explicitly as follows 

(m\ + m\ H h rn^)m\a\ + (m\ + m\ -\ V m 2 n )m\a 2 H h 

+ (m 2 1 +ml + ---+m 2 n _ 1 )mla n = 0. (2.50) 

Adding and subtracting m\ a\ to the first term of the sum, m\a 2 to the second term of the sum 
• • • etc. and finally m\a n to the last term of the sum, the equation (2.50) can be modified into 
the form 

n n n 

mj Y m]aj - ^ m ] a 3 = °> ( 2 - 51 ) 
i=i j=i j=i 

from where one can see immediately that the second equation in (2.49) is fulfilled 

n 

Y mjdj = (2.52) 

as X^/=i m ] a j = is just the first equation in (2.41) and (2.49) as well. 
The third equation in (2.41) can be written explicitly as follows 

/ 2 2 22, , 22, 22, 22, 

(m 2 m 3 + m 2 m 4 + ■ ■ ■ + m 2 m n + m 3 m 4 + m 3 m 5 + ■ ■ ■ 
+ m\m 2 n + ■ ■ ■ + m n _ 1 m„)m 1 ai + 

+ (m\m\ + m\m\ + • • • + m\m 2 n + m\m\ + m\m\ + • • • 
+ m l m l + ■ ■ ■ + m l-i m l) m i a 2 + 

+ (m\m\ + m\m\ + • • • + m 2 m 2 n + m\m\ + m\m\ + ■ ■ ■ 
+ mlml-\ h m 2 n _ 1 m 2 n )mla 3 + 
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+ (m\m\ + m\m\ + • • • + m\m 2 n _ 2 + m\m 2 n + m\m\ H 

+ m 2 m n-2 + m 2 m l 1- 

+ m n-2 m n) m n-l a n-l + (2-53) 

+ (m\m\ + m 1 m 3 + ■ ■ ■ + m\m 2 n _ 2 + m 1 m n _ 1 + m\m\ + ■ ■ ■ 

+ m 2 m l-2 + m 2 m n-l + 

+ •••+"4-2 m n-l)™« a « =0. 

Now adding and subtracting all missing terms in (2.53) from ^"i,< 2 =i TO f 1 TO i 2 S?=i TO ?°j 

»1<»2 

which in the substraction form can be written explicitly as follows 

- (ml +ml + m1-\ m 2 n )m\ai - 

- (m\ +ml + ml~\ m 2 n )m 2 a 2 - 

- (m\ + ml + ml + ■ ■ ■ m 2 n )m 3 a 3 — (2.54) 



- (mf+ml-\ hm^_ 2 + TO «) TO «-i a «-i - 

/ 2 i 2 \ 2 i 2 \ 4 

- (m 1 +m 2 -\ Vm n _ 2 +m„_ 1 )m„a„ 

and again subtracting and adding m\a\ in the first line of (2.54), m\ a,2 in the second line of 
(2. 54). ..etc., and finally m^a n in the last line of (2.54), one can rewrite (2.53) into the form 

n n n n n 

Y m ii m X Y m2 j a 3 - Y m * Y m j a i + Y m6 j a 3 = °- ( 2 - 55 ) 

*l> i 2 = 1 j — 1 i—1 j — 1 j — 1 

From this expression, taking into account the first two equations in (2.49), the third equation of 
(2.49) 

n 

rrfaj = (2.56) 

3 

follows. 

The fourth equation in (2.41) takes the following explicit form 

/ 2 2 2 i 222, 222, , 222, 

(m 2 m 3 m 4 + • • • + m 2 m 3 m n + m 2 m 4 m 5 + ■ ■ ■ + m 2 m 4 m n + • • • + 

+ ™™-2 m n-l™™) TO l a l + 

+ (m\m\m\ + • • • + m\m\m 2 n + m\m\m\ + • • • + m\m\m 2 n + • • • + 

+ ™™-2 m n-l TO ™) TO 2 a 2 + (2.57) 
+ (m\m\m\ + • • • + m\m\m 2 n + m\m\m\ + • • • + m\m\m 2 n + • • • + 
+ m 2 n _ 2 m 2 n _ l m 2 n )m 2 3 a 3 + 
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+ {ra\m^m\ + • • • + m\m%m 2 n + ra\ra\m\ + • • • + mfm^m^ + • • • + 
+ ™™-3 m n-2 TO ™) TO ™-i a «-i + 

+ (m 1 m 2 m 3 + • • • + ri{\m\m 1 n _ l + m\m^m\ + • • • + m 1 m 3 m„_ 1 + • • • + 

+ m «-3 m n-2 m ll)™2 fl « = 0. 

First, adding and subtracting all missing terms in (2.57) from 



E 



il>*2'*3 — 1 j — l 

*1 <^2 <^3 



the equation (2.57) takes the form 



E2 2 2 \ ^ 2 



*1,*2^3 = 1 7 — 1 



i 1 <i 2 ,i r ^j 



Second, subtracting and adding of all the missing terms in (2.58) from Y^%^ =1 m^m 



1 < ; 2 



2 ™2 

«2 



Sj=i one gets the equation 

n n 

E m l 1 m l 2 m i 3 ^2 m2 j a i ~ ( 2 - 59 ) 



»1<»2<»3 



- E m ii m *2 S TO J a i + E m i E TO J °J = °- 

«1,»2 = 1 j = l » = 1 j = l 

»1<»2 

Finally, additions and substractions of all missing terms in (2.59) from m f Sj=i TO j a j 

lead to the definitive form of the fourth equation in (2.41) 

n n 

i l>*2' i 3 = 1 j — l 

il <i 2 <"i3 



i l.*2 = 1 j — l i— 1 j — l j — l 

•1 <»2 

From here, taking into account the first three equations in (2.49), the fourth equation in (2.49) 



^mfa J= (2.61) 
j=i 
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follows. 

It is now easy to give a straightforward generalization of the above procedures 

i) the q-th equation in (2.41) can be decomposed into q-terms (see (2.51), (2.55) and (2.60)) 
consisting of the product of two parts, where the first part is just the sum of decreasing 
numbers of products of different vector- meson masses squared, starting from (q-1) coeffi- 
cients and ending with the constant 1. The second term takes the form Y^j=i m< j a j with 
increasing even power a starting from a = 2 up to 2q; 

ii) there is an alternating sign in front of every term in that decomposition, while the first term 
is always positive. 

Now, in order to carry out a general proof of the equivalence of the two systems of algebraic 
equations under consideration, let us assume an equivalence of (to — 2) equations in (2.41) and 
(2.49). Then, taking into account a generalization of our procedure defined by rules i) and ii) 
above, one can decompose the (to — l)-equation in (2.41) into the following form 



E m l m i--- m L-2 E m2 i a 3 - 

*li*2i*3'" - > i m — 2 = 1 j 
i l< i 2< i 3<---< i m-2 

n 

E <<"-<,-.E m &+ ,2 - 62) 

*1 J *2> < 3'" - >*m-3 = 1 J — 1 

*l<*2<*3<'"<*m-3 

n n 

+ E m n TO *2 ' ' ' m L-4 E TO i a J + • • • + 



«1,*2>«3.-'- .'m-4 = 1 

H <;.><(.,<••••<(,„ .i 



+ (-i) m - 3 E^ 2 E-? m " 2) ^ + (-i) ro - 2 E-? m_1) ^ = o. 

»=1 j=l j = l 

from where one can see immediately that the (m — 1) equation in (2.49) is satisfied 



E 



mf m 1] a 3 - = (2.63) 



as Y^j=\ m2 j a j = 0' Ej=i TO ^ a i = 0' ' ' ' ' Sj=i m 2 ^™ 1 ^ aj = are just the first (m — 2) 
equations in (2.49) assumed to be valid. 

At the end we would like to draw an attention to the proof of the equivalence of the systems 
of algebraic equations (2.41) and (2.49) from the other point of view. 

If the sums £\ =1 mpj, £ j=1 m)a 3 , £ j=1 mfa, • • • , £ J=1 "*/ a j> Ej=i m j 

aj, X^/=i ^ considered to be independent variables, then the first equation in (2.41) 

together with the modified forms (2.51), (2.55), (2.60),.. ,(2.62) form a system of (to - 1) ho- 
mogeneous algebraic equations for these variables and the equations (2.49) are just its trivial 
solutions. 
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2.5 General solution of asymptotic conditions 



In the previous paragraph we have derived two different systems of (m — 1) linear homogenous 
algebraic equations for coupling constants ratios, starting from different properties of EM FF 
Fh (t) of strongly interacting particles. 

In this paragraph, with regard to the proof of equivalence of both systems, we shall be inter- 
ested in (2.49) (the simpler one of them), though it is derived in our opinion by incorrect way by 
means of the superconvergent sum rules for the imaginary part of the EM FF to be multiplied by 
the powers of the momentum transfer squared. The coefficients in this system are even powers 
of the vector-meson masses. We find a general solution [14] of it, finally producing the VMD 
representation of Fh(t) with the required asymptotic behavior. 

First, we look for a general solution of the asymptotic conditions to be combined with the FF 
norm when FF is saturated by more vector-meson resonances than the power determining the FF 
asymptotics. 

If we assume that EM FF of any strongly interacting particle is well approximated by a finite 
number n of vector-meson exchange tree Feynman diagrams (see Fig. 2.2), one finds the VMD 
pole parametrization (2.25) and its asymptotics (2.38) is required to be determined by the power 
m. The normalization of (2.25) at t = is 



F h (0) = F . 



(2.64) 



The requirement for the conditions (2.64) and (2.38) to be satisfied by (2.25) (including also 
the results of ref. [13]) leads to the following system of m linear algebraic equations 



(2.65) 



^m 4 2r a 4 = 0, r = 1,2, ...,m - 1 



for n coupling constant ratios cij = (fihh/fi)- Therefore, a solution of (2.65) will be looked for 
m unknowns a\,...,a m and a m+ i,...,a n will be considered as free parameters of the model. Then, 
the system (2.65) can be rewritten in the matrix form 



Ma = b, 



(2.66) 



with the m x m Vandermonde matrix M 



M 



/ 1 



i 2(m-l) 



in. 



2(m-l) 



1 \ 



ml 



(2.67) 



2(m-l) i 
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and the column vectors 

a 3 



\a m J 



( F ~ ELm+1 a k \ 

En 9 
k=m+l m k a k 



En 
k= 



V-EL 



m+l m k a k 
2(m-l) 



(2.68) 



ak) 



The Vandermonde determinant of the matrix (2.67) is different from zero 



detM = ) J (mf -mfj. 



(2.69) 



3<l 



This has been proved explicitly by reducing the matrix (2.67) to the triangular form and then 
taking into account the fact that the determinant of a triangular matrix is the product of its main 
diagonal elements. 

As a consequence of (2.69) a nontrivial solution of (2.66) exists. To find the latter we use 
Cramer's Rule despite the fact that computationally Cramer's Rule for m > 3 offers no advan- 
tages over the Gaussian elimination method. However, in our case (as one can see further) all 
calculations are for the most part reduced to a calculation of the Vandermonde type determinants, 
and there is no problem to come to the explicit solutions. 

So, the corresponding solutions of (2.66) for i = 1, m are 



_ detMi 

where the matrix Mi takes the following form 



(2.70) 



1... 
m\... 



1 



Fa 



"0 — Z^fe=m+1 ak 
- Efe=m+1 m W 



1 \ 



o-EL 



k— m+1 



m k a k 



H+l 



2(m-2) 

nil 

2(m-l) 

\rrii 



in 



2(m-2) 



En 
k— m+1 

En 
k= 



2(m-2) 



m 



m+1 



2(m-l) 



2(m- 

t+1 
2(m- 

i+1 



-2) 
-1) 



m 



2(m-2) 



m 



2(m-l) 



(2.71) 



Since any determinant is an additive function of each column, for each scalar C we have 

det(A 1 ,...,CA i ,...A n ) = Cdet(A l ,...A l ,...A n ) 

and 

det(Ax, .., x fcJ 4fc, A„) = x fe det(Ai, A fc , A i+1 , ...,A n ). 



Electromagnetic form factors of strongly interacting particles and their properties 



23 



As a result, for a determinant of the matrix Mj one can write the decomposition 



detMi 



1 


1 




• ^0 


m\ 


ml 




. 


4 

m\ 






. 


2(m-2) 
1 


2(m- 


2) 


. 


2(m-l) 
1 


2(m- 

m 2 v 


1) 


. 



4 

2(m-2) 
2(m-l) 



(2.72) 



E ak 

k—m+1 



m 



mt, 



2(m-2) 2(m-2) 



"I) 



rri 



2 

2(m-l) 



2(m-2) 
k 

2(m-l) 
k 



m 



2(m-2) 
m 

2(m-l) 



from where, if in the first determinant the Laplace expansion by the entries of the column i is 
used, the explicit form is obtained 



detMi = F a (-l) 1+l \\m 2 ]"[ {mj - m?) - 

3 = 1 3,1=1 

m n m 

- n (rf-rf) e ^iiK 2 -^)- 



(2.73) 



j,l=i 

3<l,3,l^i 



k=m+l j = l 



Now, substituting (2.69) and (2.73) into (2.70), one gets the solutions of (2.66) as follows 
^l+i rr m . ™2 rr m . . . i^2_^2 

en = 



^o(-i) i+< nr=im?n m w=i K-™ 2 ) 

j¥* 3<l,3,l¥=i 



UTl=l( m l ~ m j) 
3<l 



(2.74) 



3<l,3,l¥=i 



lLv=iK-m 2 ) 

i<i 

In order to find, by means of (2.74), an explicit form of Fh(t) to be automatically normalized 
with the required asymptotic behavior, let us separate the sum in (2.25) into two parts with the 
subsequent transformation of the first one into a common denominator as follows 



F h (t) 



E 



E 

k—m-\-l 



m\a k 



(2.75) 



E^iIL^iK-tKoi 



nr=iK-*) 



m 2 k a k 



k=m+l m k 



+ E 
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Then (2.74) together with (2.75) gives 

E™i(-i) 1+i -? nr=i ™ 2 nr=iK - *) m«=i k 

F, (t) - Fn ^ ^ i<h2Ml 

h[) nr=i( m ?-*)i%i("»?-m?) 

E^iC- 1 )* -1 "*? nr=iK 2 - *) m,*=i k - ™ 2 ) el™ +1 ^ nr=iK 2 - o 
nr=iK-*)n^=iK-m|) 



+ 



+ £ ^ (2.76) 

t , 1 m fe ~ * 

k=m+l K 



The first term in (2.76) can be rearranged into the form 

rr m2 E™!(-i) 1+i nr=iK 2 -i)n ro ^=i K-m 2 ) 

1 lj=l m 7 jfr 3<U,W 



) 1+i nr=iK 2 -i)n ro ^=i 
°nr=iK-*) muK-m 2 ) 



in which one can prove explicitly the identity 

m m m m 

EH) i+i nK-t) n (w-m 2 )^ n ^ 

»=i j=i j,i=i j,;=i 

j¥* j<l,j4^'i 3<l 

leading to remarkable simplification of the term under consideration as follows 

Il m -i ™ 2 
llj=i(n»i -*) 

One could prove (2.78) by rewriting its left-hand side into the following form 

m 

£M) 1+l x 

(m 2 - t) . . . (m 2 _j - i) (m 2 +1 - t) . . . [m 2 m - t) 

m\{m\-t)... m 2 _ x {m 2 _ x -t) m 2 +l {m 2 +l - t) . . . m 2 m {m 2 m -t) 

m\{ml-t)... mf^im^-t) mf +1 (m 2 +1 -t)... m i m {m 2 m -t) 



m f m '"\m 2 -t)... mfT*\mU-t) m^K+i ~ *) • • • m^\m 2 m -t) 
m \^- 2 \m 2 -t)... mfT 2 \mU-t) m^K-i - t) . . . (m 2 m - t) 

(2.80) 

and then by using various basic properties of the determinants decomposing it into the sum of 
large number of various determinants of the same order with their subsequent explicit calcula- 
tions. Since this procedure seems to be, from the calculational point of view, not simple, with 
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the aim of a proving (2.78) let us define the new matrix 

1... 



D(t) = 



/ 1 

(ml - t) 
(ml - tf 
(ml -if 



(m\ 
(ml - 
(m\- 



-t). 
tf. 
tf. 



(ml - t) 
(mf - tf 
(m 2 -if 



1 \ 

« - 
{ml - 1? 
(m 2 m - tf 



(m 2 
\(ml 



-t" 

ty 



(ml 
(ml 



ty 
ty 



{m\ 

K 2 



ty 
ty 



{ml 



ty 
ty 



V 



(2.81) 



Denoting (to 2 — t)=Xi one gets the Vandermonde matrix 





/ 1 




1 


■ 1 \ 


D(t) = 


Xl 

x l 

r 3 


#2 • • • 
a?2 . . . 
a; 2 . . . 


x 2 . 
xl • 


x 2 
x 3 




xT~ 2 

Ur 1 


T m-2 
2 

x 2 


x?- 2 ■ 
x?- 1 ■ 


^ra-2 

-ro-1 
x m / 



(2.82) 



the determinant of which is equal just to the right-hand side of (2.78) 

mm m 

rfeiD(t) = Y[ (xi - xj) = ~[[(m 2 -t-m 2 +t) = ]J (mf - to 2 ). 



(2.83) 



j,l=i 

3<l 



3,1=1 
3<l 



3,1 = 1 
3<l 



On the other hand, if in the determinant of the matrix (2.81) the Laplace expansion by the 
entries of the first row with a subsequent pulling out of common factors in all columns of the 
subdeterminants is carried out, one gets the expression 



detD(t) = ^(-l) 1+j nK 2 -^)x 



(2.84) 



i=\ 


3 = 1 














3^1 












1 




1 




1 




i 


(ml - t) ... 


(ml 


-1 


-t) 




-t)... 


{ml - 1) 


(ml - t) 2 ... 


(m 2 _ 


1 


-ty 




-tf... 


{m 2 m - tf 


(m\-tf ... 


(m 2 _ 


1 


-ty 




-tf... 


(ml - tf 


(ml - t) m ~ 3 . . . 






t yn-3 


K 2 +i - 


t) m - 3 . . 


(ml - t) m - 3 


(ml - t) m - 2 . . . 






fyn-2 


K 2 +i - 


ty n - 2 . . 


(m 2 m -tr- 2 



Then calculating explicitly the determinant in (2.84) by using again the denotation Xk = (m\—t) 
for k = 1, i — l,i + 1, ..to, one finally obtains 



detDftl^i-ll'+'HK 2 -;) I] {m 2 -m 2 3 ) 



(2.85) 



3 = 1 
3^1 



3,1=1 
3<l,3,l^i 



26 



Analyticity in a phenomenology of electro-weak structure of hadrons 



just the left-hand side of (2.78) and in this way the identity under consideration is clearly proved. 

The second and third term in (2.76), transforming them to a common denominator, can be 
unified into one following term: 

9 TT m I ™2 jA Tf m (™2 ™2\ 



V 



m 2 fe nr=iK 2 -*)n^=iK-^ 2 ) 



(2.86) 



K-t)E™i(-i)^fnr=i(m; 



+ 



(m 



K - t) nr=iK - *) n™/=iK - m|) 



>a fe 



the numerator of which is exactly the Laplace expansion by the entries of the first row of the 
determinant of the matrix of the (m + 1) order 



N(t) = 



K - 1) 

{ml- if 
{ml-tf 



{m 2 k 



, - t) m 

v k - v 



\ 





-t) . 


(m 2 m 


-t) 


{m\ 


-t? ■ 


■ (m 2 m 


-tf 


{ml 


-tf . 


■ (m 2 m 


CO 


{ml - 


t)" 1 - 1 . 


■ K- 


j\m— 1 


{m\- 


-t) m . 


• ("4 





(2.87) 



If we define the new matrix of the (to + 1) order 



/ 



R(t) = 



( m l 


~t) 


{ml 


-t) . 


(m 2 m 


~t) 


( m l 


-tf 


{ml 


-tf . 


■ (m 2 m 


-tf 


(ml 


-tf 


{ml 


-tf . 


■ (m 2 m 


-tf 


(m 2 k - 


t yn-l 


{ml — 


t)™- 1 . 


■ (m 2 m ~ 


^.yn— 1 


V (ml - 


-t) m 


(ml- 


-t) m . 


■ (m 2 m - 


-t) m J 



(2.88) 



then for the determinant of both matrices, (2.87) and (2.88), the equation 
deiN(t) - t ■ detR{t) = detS{t) = 



(2.89) 



is fulfilled under the assumption that defS(t) is obtained by multiplication of the first row of 
detH{t) by t, and the substraction of the resultant determinant from N(t) is carried out explic- 
itly. 

There is valid also a relation 



detN(O) = 



(2.90) 



as in cfeiN(O) (like in detS{t)) the first two rows are identical. 

Now, in order to arrange the numerator of (2.86) conveniently, we write cfeiN(t) in the form 



deiN(t) = t ■ detR{0) - deiN(O), 



(2.91) 
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taking into account (2.89), (2.90) and the identity 

deiR(t) ee detR(O). 

In (2.91) we apply the Laplace expansion by entries of the first row to cfeiR(O) 
and deiN(O), separately. As a result, one gets 



(2.92) 



deiN(t) = t ■ 



m? 



m\ 



ml 



2(m-l) 
V 

m 2ra 



2(m-l) 



m 



2m 



2(m-l) 

m m 
m 2 ™ 



+ 



+ 



m l 



2 

mf 
mf 


m| 


m^ 
mi 


2(m-l) 

m? m 


2(m-l) 
TO 2 V 


2(m-l) 
™2m 



E 

»=i 



(-i)x 2 



or calculating explicitly the corresponding subdeterminants 

m m 

de*N(t) = (t - m|) JJ JJ (mf - to?) + 



3<l 



(2.93) 



4 


2 

mj . . . 

4 

mf . . . 


m|_! 

4 


m 1+l 
4 

"¥+1 




2(m-l) 

m| m 


2(m-l) 

m? m . . . 


2(m-l) 
2m 


2(m-l) 

2m 
TO i+l 


2(m-l) 

• m m 

m 2m 
ln m 



m l 


2 

mf . . . 


m|_! 






m t 


4 

mf . . . 


4 






2(m-l) 

m fe 
m 2m 


2(m-l) 

m x v ' . . 
mf m . . . 


2(m-l) 


2(m-l) 

2m 
m i+l 


2(m-l) 

• m m 

" l m 



(2.94) 



i=i 



j.;=i 



Substituting the latter into (2.86) one obtains 



v f nr=t 
fc ir + it"nr=i(-: 



+ 



nm 



n?^i(™j -™fc)' 
.7'7 ^ » 



m?) 



•Ofc (2.95) 



and combining this result with (2.79), one gets the form factor Fh{t) to be saturated by n-vector 
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mesons (n > m) in the form suitable for the unitarization 



F h (t)=F 



n; 



+ 



E 

fe— m+l 



E 



n 



(m 



*) nr=iK 2 



*) 11^=1 ("» 



t) 



(2.96) 



for which the asymptotic behavior (2.38) and for t = the normalization (2.64) are fulfilled 
automatically. 

The asymptotic behavior in (2.96) is transparent. However, for the normalization (2.64) the 
following identity 



2J 



nrn / o 
j=i{mj-m 



2 _ „2\ 



(2.97) 



has to be valid in the second term of (2.96) generally. 

For to = 2, 3, 4, 5 it can be proved explicitly. And for an arbitrary finite to it follows directly 
from (2.86), the numerator of which is exactly the Laplace expansion by the entries of the first 
row of the determinant of the matrix (2.87). Then, just relation (2.90) causes the term (2.86) 
and also (2.95) at t = for arbitrary nonzero values of ak to be zero. Hence, every term in the 
wave-brackets of (2.95) for t — has to be zero and this is true if and only if identity (2.97) is 
fulfilled. 

Now we consider the case of equations (2.65) for n = m. Then it can also be rewritten into 
the matrix form (2.66) with the m x to Vandermonde matrix M (2.67) and the same column 
vector a, but with the b vector of the following form 

o 




So, the corresponding solutions are again looked for in the form 

detMi 



(2.98) 



(l; = 



detM 



with the matrix Mi 



M i = 


/ 1 

2 

m\ 




1 

2 

m i-l 

n*i-i 


Fo 




1 

m t+l 




2(m- 
TO 1 V 

2(m- 


2) 
1) 


2(m-2) 

• m i-l 

2(m-l) 

• "Vi 






2(m-2) 
m i+l 

2(m-l) 
m i+l 



to; 
mi 



rn 



m 



2(m-2) 
m 

2(m-l) 



(2.99) 
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and its determinant to be 



detMi = F (-1) 1+1 f[m? [j {mf-m*). 



(2.100) 



The solutions 



= Fn 



(-l) 1+, Ilp^nVi .K-m|) 



ILy=i(m?-m?) 
nr=imf(-l) 1+i 



= F n 



n?=i(mj-m?)(-l)«-i 



(2.101) 



are then completely expressed through only the masses of m vector-mesons by means of which 
the considered FF is saturated. 

Substituting the solutions (2.101) into the VMD parametrization (2.25) of the EM FF one 
comes to the following representation 



Fh(t) 



n 



'nr=i( 



■im 



t) 



(2.102) 



dependent only on the considered vector-meson masses and the required asymptotic behavior 
(2.38) is transparent to be fulfilled automatically. 

The third case with the (to — 1) linear homogeneous algebraic equations for the n 
(n > m) coupling constant ratios without any normalization of FF appears naturally in the 
determination of the so-called strange FF behaviors of a strongly interacting particles with the 
spin s > from the isoscalar parts of the corresponding EM FFs, as we shall see later on. 

Then, we have only the equations 



£ 

i=l 



2r 

m 4 Oi 



0. 



1,2, ...to- 1 



(2.103) 



which can be rewritten in the matrix form (2.66) with the (to — 1) x (to — 1) matrix M 



/ m\ 



M 



4 



and the column vectors 



a = 



o 2 
a 3 



\a m -ij 



ml 

4 

m% 



\ 



m 



m—l 



2(m-l) 2(m-l) 



TIL, 



m 



2(m-l) 

TCI— 1 



b = 



/ ~ 2^k= m m k a k 

En 4 
k=m m k a k 

En g 
k=m m >k 



•sr^n 2(m— 1) 

\-E fe=m "v a */ 



(2.104) 



(2.105) 
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The determinant of the matrix M 



m—l m— 1 

detM = J| m i II ( m i - (2.106) 



is different from zero, and thus, a nontrivial solution of (2.103) exists 

detMi 

ai = 

detM 

where the matrix Mi takes the form 
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(2.107) 



Then employing the basic properties of the determinants one gets 
detMi = 



ak 
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(2.108) 



As a result, 



En 2 t — ttm — 



o t = 



nm-i / 2 



U-iy-'UT^i^-ml) 

3^1 



(2.109) 
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or finally, 



E 

k—m 



i=i K - m i) 



< nr=i K-m?) 



afe, i = 1, 2..., m — 1. 



Now substituting (2.1 10) into 



Em— 1 rrm-1 / 
i=l 11.7=1 \ 



+ 



f 



k=m k 

and transforming both terms into a common denominator one gets the relation 



Fh{t) = E 



n^K-t)n^K-^) 



,-, ; ,vnr=7K-i)n,T=iK-m^) 

3<l 



+ 



,2 j\ ^ m_1 / T\iTf m ~ 1 f^2 4-\ TT m— 1 C™2 „,2-i rrm-1 
lj=l t m j 



(2.110) 



(2.111) 



(2.112) 



+ 



t) itv=i ( m ? - m j) n^i 1 ^ 

j<l,j,ljti jjti 



m k) 



nr^K-^n^iK-m?) 



r a fc, 



in which the numerator of the first term under the sum is just the Laplace expansion by the entries 
of the first row of the determinant of the matrix T(t) of the m order 

/ 1 1 ... 1 \ 

(m 2 k -t) (mj-t) ... (ml^-t) 
(ml-tf {m\- if ... {m^-tf 



T(t) = 



\K-tr- 1 (m?-*)" 1 - 1 ... (m^i-tr-v 



One can immediately prove that 

deiT(t) = detT(O). 
Then calculating detT(0) explicitly 



(2.113) 



(2.114) 



detT(0) = Yl ( 

i=i 



2 2 

- m fe 



) n(m?-m?) 



(2.115) 



and substituting the result into (2.1 12) instead of the numerator of the term in the wave-brackets, 
one finally obtains the parametrization 



F h (t) = 



mi 



-a k 



for which the asymptotic behavior (2.38) is fulfilled automatically. 



(2.116) 
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2.6 Analytic properties of electromagnetic form factors 

In principle there are two sources on the analytic properties of FFs of strongly interacting parti- 
cles. 

The first one resides in the exact proof of the analytic properties of FFs starting form the first 
principles of the local QFT. In this way the analytic properties of the pion FF were proven [15], 
though at present days in connection with quark-gluon structure of hadrons this pretentious proof 
is possible to accept in such approximation, in which the considered hadron can be brought 
into compatibility with the local quantum field. Unfortunately, there are no exact proofs of the 
analytic properties of FFs of any other hadrons, though they are utilized practically for many 
years. 

There is a general belief that all EM FFs are analytic functions in the complex plane of the 
momentum transfer squared t besides a cut from the lowest branch point to on the real axis to 
+oo. The positions of the branch points are found by an investigation of the analytic properties 
of Feynman diagrams [16-18] representing separate terms of a formal series of EM FFs obtained 
in the framework of a quantum field perturbation theory. 

As a consequence of the hadron EM current to be Hermitian, all EM FFs are real on the real 
axis for t < to- Then by an application of the Schwarz reflection principle to EM FFs one finds 
the so-called reality condition 

F* h {t) = F h {t*) (2.117) 

reflecting the reality of the EM FFs on the real axis below to and the relation of values of EM 
FFs on the upper and lower boundary of the cut 

F*(t + ie)=F h (t* -te), e«l (2.118) 

automatically. 

The discontinuity across the cut is given by the unitarity condition 

Y.{{hh\M0)\0) - (0\J,(0)\hh)*} = J2( h ~ h \T + \n)(n\J^0)\0), (2.119) 

n 

where the sum in (2.119) is carried out over a complete set of intermediate states allowed by 
various conservation laws and T + means Hermitian conjugate amplitudes. 

Moreover, just from the unitarity condition (2.1 19) it follows that there is an infinite number 
of branch points on the positive real axis between t = t and +oo, which always correspond to 
a new allowed intermediate state n in (2.1 19). In order to fulfil the reality condition (2.1 17), the 
cuts associated with these branch points are chosen to be extended to +oo along the real axis. 

Further we demonstrate the main principles of investigation of the analytic properties of 
Feynman diagrams. 

Individual terms of the formal series of the perturbation theory possess the analytic proper- 
ties which are in no contradiction with first principles of the local QFT. An arbitrary Feynman 
diagram can be represented by the integral 

I £ (t) = f d*k 1 ...ct l k l =r i . 2 N 2 e > 0, e « 1, (2.120) 
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where qj is four-momentum of j propagator at the diagram and mj is the corresponding 
mass. N at the numerator is representing a spin structure of Feynman diagram and it does not 
contribute to an appearance of any singularity. The variables kj (j — I,..., I) are four-momenta 
related to I independent loops to be chosen arbitrary, nevertheless with regard to the conservation 
law of four-momenta. 

At the investigation of the analytic properties of Feynman diagrams appears to be the most 
suitable the so-called a-approach. By denotation qj + mj = Uj and an application of the 
Feynman relation 



Now, to find singularities of an individual term of the perturbation expansion means to scru- 
tinize the analytic properties of integrals of the type I e (t) in (2.122). 

The problem of a search of singularities of the general integral (2.122) has been reduced by 
Landau [16] and Cutkosky [19] into a series of necessary conditions (so-called Landau-Cutkosky 
equations) to be formulated in the following way. The integral I s (t) possesses a singularity: 

• if either ctj = or q? = —m? for every j inner line; 

• Ym=i a iQi = f° r ever y independent loop. 

If one is restricted to a search of singularities only on the physical sheet of the Riemann surface, 
then to the above-mentioned conditions another one is joined in 

• otj have to be real positive numbers. 

From the first condition it follows, that at the diagram generating a singularity either every 
inner four-momentum is on the mass shell or the corresponding aj is equal zero. In the last case 
the four-momentum qj does not appear in the second condition, so a presence of the inner line 
at the diagram does not influence the created singularity. Or in other words, the same singularity 
can be found from the Feynman diagram, in which the inner line belonging to a.j at the primary 
diagram, is substituted by point. 

Diagrams, obtained from the primary diagram by removing one or more inner lines (the 
corresponding vertices are joined into one point) are called to be reduced diagrams. 

Further there are described two types of Feynman diagrams, a triangle presented in Fig. 2.3 
and a two-point diagram presented in Fig. 2.4. Both of them are creating branch points, the 
triangle the anomalous and the two-point the normal (physical) threshold. The two-point diagram 
can in another time appear to be joined with some triangle or more-component diagram. 




(2.121) 



the integral (2.120) can be transformed into the following form 




(2.122) 
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Figure 2.3. Triangle diagram as reduced from primary diagram 




Figure 2.4. Two-point diagram as reduced from primary diagram 



At the determination of the singularity one can start from the two-point diagram (see Fig. 2.4), 
from which it follows 



92 = ~ m h <?3 = - m h a 2l2 + a 3 q 3 = 0. (2.123) 

By a multiplication of the last equation by q 2 and q 3 one after the other one obtains the system 
of algebraic equations 



-ot 2 m\ + a 3 q 2 q 3 = 
"29293 - ct 3 ml = 0, (2.124) 

which possesses nontrivial solution for aj at that time and only at that time when the deter- 
minant of the system is equal zero, i.e. 

det (-ml <12QA =0 (2-125) 
V<?2<73 -mi) 

orra\m\ - (q 2 q 3 ) 2 = 0. 

By utilization of the conservation law of the four-momenta Q = p 2 — pi =93 — 92 and 
t = —Q 2 one gets the expression 



9293 = - [*- (m| + m§)], 
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which together with the relation for the determinant gives the following quadratic equation 

t 2 - 2(m 2 2 + m\)t + (ml - m\f = 0. (2.126) 
Its solution gives 

t=(m 2 ±m 3 ) 2 , (2.127) 

in which only the root with the positive sign fulfills the third condition. As a result the two- 
point diagram in Fig. 2.4 creates the branch point on the physical sheet of the Riemann surface 
at t = (m 2 + m 3 ) 2 to be known as the normal threshold. 

Further we are interested in the triangle diagram in Fig. 2.3, where 

2 2 2 2 2 2 

9l = - m l> 9 2 = ~ m 2: 9 3 = - m 3 

and 

a 1 q 1 + a 2 q 2 + a 3 q 3 = Q. (2.128) 

From the last equation by a subsequent multiplication by q\, q 2 a q 3 one obtains the system 
of three algebraic equations 

- a x m\ + a 2 q x q 2 + a 3 q ± q 3 = 

ai9i 92 - a 2 ml + a 3 q 2 q 3 = (2.129) 
ai<7i<?3 + a 2 <?2<?3 _ a 3 m i = °> 

possessing a nontrivial solution only in the case of the zero determinant 

(-m\ qiq 2 qiq 3 \ 
det qi q 2 -ml q 2 q 3 = 0, (2.130) 
\9i93 9293 -ml) 

or 

-m\m 2 2 m 2 3 + 2(q l q 2 )(q 1 q 3 )(q 2 q 3 ) + ml(q 2 q 3 ) 2 + m^q^f + m^q^f = 0. 

By utilization of the four-momenta conservation laws at the corresponding vertices at the 
triangle diagram 

P2 - Pi = 93 - 92 

Pi -9i -93 (2.131) 
Pi = 9i - 92 
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and the relation p\ = p\ = —M 2 one obtains for scalar products of four-momenta appearing 
at the determinant 

(9293) = ^[t- (ml +m§)] 
(q^) = \[M 2 -(ml + ml)] 
(«i<fc) = \[M 2 - (m? + ml)}, 

on the base of which the condition (2.130) leads to the quadratic equation 

m\t 2 + |[M 2 - (m\ + m 2 2 )][M 2 - {m\ + ml)] - 2m 2 l (m 2 2 + m 2 3 )}t + 

+ ml(ml - mj) 2 + m%[M 2 - (ml + m 2 )} 2 + ml[M 2 - (ml + ml)] 2 - (2.132) 
- (ml + mj)[M 2 - {ml + m 2 2 )][M 2 - (ml + ml)] = 0. 

Its solution takes the form 

t = (ml + m 2 ) - -L [M 2 - (ml + m 2 2 )][M 2 - (m? + m 2 )] 



± 2mJ \/ 4to ? to 2 - [M 2 - ( m l + m 2 2 )] 2 ^mlm 2 - [M 2 - (ml + m 2 3 )] 2 , (2.133) 

determining singularities of EM FFs from the triangle diagram to be called the anomalous 
thresholds. The latter, unlike the normal thresholds, do not correspond to some physical pro- 
cesses. 

The anomalous thresholds are found on the physical sheet of the Riemann surface only in 
the case of the real positive solutions aj of the system of equations (2.129). This condition is 
fulfilled by the solution (2.133) with the positive sign. 

The most effective approach for a determination of the position of the anomalous threshold 
of the triangle diagram seems to be a geometrical way of a solution of the Landau-Cutkosky 
equations. It consists in the following. If Q is expressed as i\fi, then from the four-momenta 
conservation law (2.131) follows 

Pi +iVt = p 2 , 

to be pictured graphically in Fig. 2.5 and called a dual diagram to Fig. 2.3. 

Every inner triangle on Fig. 2.5 represents the four-momenta conservation law at the vertices 
of the diagram on Fig. 2.3. From its form one can reveal wether the singularity is placed on the 
physical sheet or on one of the un-physical sheets of the Riemann surface. 
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The Landau-Cutkosky equation (2.128), expressing the linear dependence among qi, q 2 , 53,, 
requires all three vectors to be in one plane. The lengths of qj are equal to the corresponding 
masses rrij and the lengths of pi a p 2 to the mass M. From the latter there is unambiguous 
boundary on the form of the outer triangle in Fig. 2.5. 

The anomalous threshold is determined by an expression for t in Fig. 2.5 to be found by the 
methods of elementary geometry. 

From the form of the dual diagram one can conclude, whether the anomalous threshold is 
located on the physical sheet or on the one of the un-physical sheets of the Riemann surface. 

The condition of a positivity of the parameters aj, securing for the singularity to be on the 
physical sheet, leads the point to be inside of the outer triangle in Fig. 2.5. If the point 
is outside of the outer triangle, the corresponding anomalous singularity is placed on some un- 
physical sheet. 

So, the position of the anomalous threshold of EM FF on the physical sheet of the Riemann 
surface, generating by the triangle diagram, is numerically determined by (2.133) with the posi- 
tive sign. 



2.7 Unitary and Analytic model of electromagnetic structure of hadrons 

In the previous paragraphs we have solved the conflict between the uniform asymptotic behavior 
of the VMD model (2.29) and the asymptotic behavior (2.24) of EM FFs as predicted by the 
quark model of strongly interacting particles by finding the three VMD parametrizations (2.96), 
(2.102) and (2.116) fulfilling the quark model asymptotic behavior (2.24) automatically. Every 
of them are compound of the products of the resonant terms 



ml 



t 



(2.134) 



which further will be unitarized by an incorporation of the well known analytic properties of 
EM FFs consisting of an infinite number of branch points on the positive real axis. They are 
branch points corresponding to normal (given by the unitarity condition) and anomalous (given 
by allowed triangle diagrams) thresholds generating many-sheeted Riemann surface. The first 
sheet is called physical, all other sheets of the Riemann surface are unphysical. 

Practically, further we are restricting ourselves to a two-square-root-cut approximation of 
the analytic properties generating the four-sheeted Riemann surface. Then any resonant stay is 
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always associated with a pair of complex conjugate poles on unphysical sheets to be generated 
by the branch points on the positive real axis of the t-plane. 
In order to transform (2.134) into one analytic function 

i) with two square-root branch points on the positive real axis, 

ii-j with two pairs of complex conjugate poles on unphysical sheets corresponding to the res- 
onance r, 

one proceeds as follows 

- first, the nonlinear transformation 

ai35) 

with to - the square-root branch point corresponding to the lowest threshold and tj„ - an 
effective square-root branch point simulating contributions of all other relevant thresholds 
given by the unitarity condition is applied, which automatically generates the relations 

2 . 4(tj„ — to) /0 n--\ 

mt = t + . jr — — (2.136) 



[l/W r0 - T^o] 2 



and 



= to+ M/ 4 T ^" §^ (2.137) 
[1/W N - W N \ 2 

- then relations between W r o and W* Q are utilized 

- and finally, the instability of the resonance is introduced by its non-zero width T r ^ 0. 
The application of the (2.135)-(2.137) to (2.134) leads to the following factorized form 

ml _ ( 1 - W 2 \ 2 (W N - W r0 )(W N + W r0 )(W N - l/W ra )(W N + l/Wta) 



W%) (W- W r0 )(W + W r0 )(W - l/W r0 )(W + 1/W r0 ) 



(2.138) 



with asymptotic term (the first term, completely determining the asymptotic behavior of (2. 134)) 
and on the so-called finite-energy term (for \t\ -4- oo it turns out to be a real constant) giving a 
resonant behavior around t = m 2 . 
One can prove 

a-, if m 2 r - T 2 /4 < t in => W rQ = -W* 

(2.139) 

b) if ml - T 2 /4 > t in => W ra - l/W*o 
which lead the eq. (2.138) in the case a-j to the expression 

ml _ ( i - w 2 \ 2 (w N - w r0 )(w N - w r * )(w N - i/w r0 )(w N - i/w; ) 



m 2 -t \l-W 2 ) (W-W r0 )(W-W? )(W-l/W r0 )(W-l/W? ) 
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and in the case b) to the following expression 

m 2 r f i - w 2 \ 2 (w N - w r0 )(w N - w; )(w N + WtoHWn + w; ) 
m*-t \i-w%) (w- w r0 )(w - w* )(w + w r0 )(w + w; ) • ' 

Lastly, introducing the non-zero width of the resonance by a substitution 

ml -> (m r - T r /2) 2 (2.142) 
i.e. simply one has to rid of "0" in sub-indices of (2.140) and (2.141), one gets: in a) case 

m 2 r (\-w 2 \ 2 (w N -w r )(w N -w;)(w N -i/w r ){w N -i/w;) 



l-W^J (W- W r )(W - W*){W - l/W r )(W - 1/W*) 

\-w 2 \ 2 

' L(W r ) (2.143) 



1-W 2 N 



and in b) case 

i - w 2 \ 2 {w N - w r )(w N - w;){w N + w r )(w N + w*) 
i-w 2 ) (w - W r ){W -w*){w + w r ){w + w*) 

i-w 2 \ 2 

iTwf) H <W (2-144) 

where no more equality can be used in (2.143) and (2.144) between the pole-term and the trans- 
formed expressions. 

The latter are then analytic functions defined on four-sheeted Riemann surface, what can be 
seen explicitly by the inverse transformation to (2.135) 



Wit) = i \ K ' V ' = \ ■ (2.145) 

// \ 1/2 / \ 1/2 // \ 1/2 / \ 1/2 

vOvO + (^) + v(^) -(^) 

These expressions have two pairs of complex conjugate poles on 
a-j case (2.143) - the second and fourth sheets 
b) case (2.144) - the third and fourth sheets, 

describing always one resonance under consideration, and the asymptotic behavior <~ 1 /t, to be 

completely given by the asymptotic term ( ) > more specifically, by its power "2". 

Here we would like also to note, that expressions (2.143) or (2.144) (it depends on the fact 
if the resonance is under the threshold t in or above it) are more sophisticated analogue of the 
Breit-Wigner form and they can be applied to determine resonance parameters m r , T r from 
experimental data on EM FFs in the region of the resonance under consideration. 
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3 Electro-weak structure of nonet of pseudoscalar mesons 

There is only one scalar function of the momentum transfer squared t completely describing the 
electro-weak structure of every member of the nonet of pseudoscalar mesons. 

The EM FFs are identically equal to zero for ir°, rj and rf for all t from —00 to +00, if one 
takes into account the relation (2.6). Then practically one has to consider only three EM FFs for 
psedoscalar meson nonet corresponding to the charged pions, charged kaons and neutral kaons, 
respectively. 

The pseudoscalar mesons are the bound states of qq pairs and as a result of (2.24) the asymp- 
totic behavior of their FFs is 

*H*)|t|-oo ~ \ (3-D 

to be identical with the VMD model (2.25) asymptotics. 

Consequently, the pseudoscalar meson FF in the framework of the U&cA model can be rep- 
resented by a sum of i-terms (2.143) of resonances below the threshold t in and j-terms (2.144) 
of resonances above the threshold t in as follows 

/ 1 - W 2 \ 2 

MW{t)\ = f^-^) x (3.2) 

Wi){W N - W*)(W N - l/Wi)(W N - 1/W*) , 




E 



Wi){W -W*){W- l/Wi){W - 1/W*) 
{W N - Wj)(W N - W*)(W N + Wj)(W N + W*) 



-(fipp/fi)+ 



{W - W.j)(W - W*){W + Wj)(W + W*) 



(fjpp/fj) 



which takes into account just n = i + j vector-meson resonances with quantum numbers of the 
photon. It is analytic in the whole complex t-plane except for two cuts on the positive real axis. 



3.1 Electromagnetic form factor of charged pions 

The 7r- meson is the lightest hadron with spin S = and isospin 1 = 1, i.e. three charge 
(positive, negative and neutral) states of the pion do exist. We consider the pion to have an 
internal (commonly quark-gluon) structure that for the time being we do not know to reproduce 
theoretically but that can be parametrized in some rather general way. 

The virtual photon-pion vertex with the internal structure of the pion is graphically presented 
in Fig. 3.1. 

The blob means various virtual processes caused by the strong interactions and theoretically 
it is expressed by a matrix element (p2|^ M (0)|pi) of the electromagnetic current j[f M (x) of 
the pion, that can be decomposed into a sum of products of invariant coefficients and linearly 
independent covariants constructed now only from the four-momenta pi , p 2 as the spin of the 
pion is zero. 

The corresponding coefficients, so-called invariant FFs, are functions of only one invariant 
variable. 
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Figure 3.1. The virtual photon-pion vertex 



Really, from the conservation of four-momenta at the virtual photon-pion vertex 

Pi+q = P2 (3.3) 

and by a multiplication of both sides of the latter by one after the other four-momenta q, p\ and 
P2, one gets a system of the following three algebraic equations 

(qpi) + q 2 = [qp2] 

p\ + {qpi) = fan) (3.4) 
(P1P2) + (qp2) = p\ 

for four unknown invariant variables, (qpi), {qpi), q 2 , and (P1P2), provided p\ — p\ = — toJ, 
where is the pion mass. 

By a solution of (3.4) one can express three invariant variables by means of an arbitrary value 
of the fourth one. The latter is usually chosen to be just the four-momentum transfer squared 

t = q 2 = -Q 2 . 

A decomposition of the matrix element of the pion electromagnetic current looks as follows 

( P2 |jf»| Pl ) = e[Ff (t)(pi +p 2 ) M + G* (t)( P2 (3.5) 

where instead of two independent four- momenta pi, p 2 we have used their suitable combinations, 
(Pi + P2)p. and (p 2 -pi)n = gy, and F®(t), G® (t) are invariant FFs. 

The gauge invariance of the EM interactions implies the current conservation 

d^ M (x) = (3.6) 

from where the following condition 

qMrf M (o)\ Pl ) = (3.7) 

has to be fulfilled. Since q^(pi+P2)fi =p\—p\ = 0, the condition (3.7) will be valid if in (3.5) 
(t) — 0. Then the matrix element of the pion EM current takes finally the well-known form 

(jp 2 \J* M (0)\ Pl ) = eff (t)(pi + P2 )^ (3.8) 
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where (t) is the so-called pion EM FF. 

On the other hand, making use of transformation properties of the EM current operator 
J^ M (x) and the one-particle state vectors of pions, \ir~) and \tt°), with regard to the 

all three discrete C, P, T transformations simultaneously, one can derive relations as follows 

(7r-|jf"(0)|7r-) = -<^ + |J^ M (0)|7r+) (3.9) 

and 

(7rVf»k > = -<7r°|jf"(0)|7r°>. (3.10) 

Then a substitution of the parametrization (3.8) into (3.9) and (3.10) provides for the EM FFs of 
charged and neutral pions relations 

FfL(t) = -F£(t) (3.11) 

and 

F${t) = 0, (3.12) 

respectively, for all values of the four-momentum transfer squared — oo < t < +oo. 

Since the charge e appears as a prefactor in (3.8), the pion EM FF at t = is normalized to 

one 

i^(0) = 1. (3.13) 

In the framework of the axiomatic quantum field theory one can prove [15] that F 7T ± (t) is 
an analytic function in the whole complex t- plane, besides a cut from 4m^ to +oo. The same 
analytic properties can be derived within the framework of QCD [20]. Further, as a consequence 
of the pion EM current J^{x) to be Hermitian, the pion EM FF is real on the real axis for 
t < 4ra^. Then by an application of the Schwarz reflection principle to the pion EM FF one 
finds the so-called reality condition 

[F? (t)}* = F*{t*) (3.14) 

reflecting the reality of the pion EM FF on the real axis below 4m^ and the relation of values of 
the pion EM FF on the upper and lower boundary of the cut 

[F*{t + ie)Y =F?(t-ie), e«l (3.15) 

automatically. 

The discontinuity across the cut is given by the unitarity condition 



Yi {(^-|^ M (0)|0) - (0|J,f M (0)|7r+O*} - 52(* + *-\T+\n)(n\Jj? M (0)|0), (3.16) 
where 



(tt+tt- 1^(0)10) = e(k w+ K-)^{t) (3.17) 
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is a matrix element of the pion EM current defining the pion EM FF in the time-like (t > 4m 2 ) 
region. The sum in (3.16) is carried out over a complete set of intermediate states allowed by 
various conservation laws and T+ means Hermitian conjugate amplitudes. 

Moreover, the unitarity condition (3.16) tells us that there is an infinite number of branch 
points on the positive real axis between t = 4m 2 and +oo, which always correspond to a new 
allowed intermediate state n in (3.16). In order to fulfil the reality condition (3.14), the cuts 
associated with these branch points are chosen to extend to +oo along the real axis. 

We would like to note that one can specify [16-18] the same singularities of the pion EM FF 
by means of an investigation of the analytic properties of the corresponding Feynman diagrams 
of a formal pion EM FF perturbation series. 

Consequently, the first branch point of the pion EM FF is at t = 4m 2 , the second one at 
t = 16m 2 , then at t = (m,o + m^) 2 , t = Am 2 K , etc. 

If we restrict ourselves in (3.16) only to \n) = \ir + ir~), then the elastic unitarity condition 
of the pion EM FF is obtained 

1 {F?(t + ie) [F*(t + ie)}*} = [A\{t + ie)]*F?(t + ie) (3.18) 

where A\(t + ie) is the P- wave isovector tttt- scattering amplitude. 

If the complete set of intermediate states in (3.16) is taken into account then the pion EM FF 
unitarity condition can be written in the form 

i {F?(t + ie) [F*(t + ze)}*} = [A\{t + ie)]*F*{t + ie) + a(t + ie) (3.19) 

where a(t + ie) represents all higher contributions. 

Now we prove that the lowest singularity of the pion EM FF at t = Am 2 is a square-root type 
branch point. 

The idea consists in the analytic continuation of the pion EM FF through the upper and lower 
boundary of the cut between t = Am 2 and t = 16m 2 on the next sheet of the Riemann surface 
by means of the elastic unitarity condition (3.18) and in a subsequent comparison of obtained 
functional expressions. 

By using the reality condition of the pion EM FF (3.14) and similar condition [21] 

[A{(t + ie)\* = -A\(t-ie) (3.20) 

for the P- wave isovector tttt - scattering amplitude, one can rewrite (3.18) into the form as 
follows 

g^'- i+t^A) - (3 - 21) 

If we denote by [F^(t — ie)] + and [F^(t + ie)]~ the analytic continuation of the pion EM 
FF on the Riemann sheets achieved through the upper and lower boundary of the elastic cut, 
respectively, then as a consequence of a continuity of the pion EM FF we have 

\F?(t-ie)]+ = F?(t + U) (3.22) 
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[F^(t + te)}- = F?{t-ie). (3.23) 
Then by a substitution of (3.22) into (3.21) one gets 

- i^gfe <3 - 24) 

On the other hand, by a complex conjugation of (3.18) and the reality conditions (3.15) and 
(3.20), the unitarity condition of the pion EM FF on the lower boundary of the elastic cut is 
obtained 

i {F?(t - ie) [F*(t ze)}*} = {A\{t - ie)}* F* \t re), 
from which the relation follows 

F?(t - ie){l + 2iA\(t + ie)} = F^(t + ie) 
or by using (3.23), finally one gets 

■g" +fc "-- iTOW a25) 

Comparing (3.25) with (3.24) we see that by the analytic continuation of the pion EM FF 
through upper and lower boundary of the elastic cut one gets the identical functional expressions. 
The latter convinces us that one has come on the same sheet of the Riemann surface, which 
further will be called the second one and denoted by II. It follows just from here that the branch 
point t = 4toJ is of a square-root type, i.e. it generates just two sheets of the Riemann surface, 
on which the pion EM FF is defined. 

Moreover, the analytic continuation of (3.24) and (3.25) leads to the same expression of the 
pion EM FF on the second Riemann sheet as follows 

I F - E M1" - Y^wy <3 ' 26) 

One can read out from (3.26), that the pion EM FF on the second Riemann sheet has all the 
branch points of the pion EM FF on the first sheet and besides also the branch points of the P- 
wave isovector tttt- scattering partial amplitude A\(t). 

It is well known [21] that the analytic properties of A\(t) consist of the right-hand unitary 
cut for 4m \ < t < oo and the left-hand dynamical cut for — oo < t < 0. From the latter and 
(3.26) it follows that the pion EM FF on the second Riemann sheet has the left hand cut too. 

The unitarity condition in the language of the Sl(t) matrix, connected with the amplitude 
A\ (t) in (3.26) by the relation 



Sl{t) = l + 2iA\(t), 



(3.27) 
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has the form 

Sl(t)[Sl(t)]* = 1 (3.28) 
from where one can parametrize 

Sj(t)=e 2i4 i<*> (3.29) 

with 8\(t) to be the P- wave isovector 7T7r-scattering phase shift. By means of a substitution of 
(3.27) into (3.28) one gets the unitarity condition for A\(t) in the form 

ImA\{t) = \ A\{t)\ 2 (3.30) 

which is exactly valid for 4m 2 < t < 16m J. However it follows just from data [22-24] on 
A\{t) that it works approximately up to lGcV 2 . 

On the other hand by a combination of the relations (3.29) and (3.27), the following parame- 
trization of A\{t) is obtained 

A\(t) = (e 2 ^ (t) - l) - e^ (t) sin^(i) (3.31) 

which is consistent with the unitarity condition (3.30) automatically. 

We note, that the partial wave amplitude A\{t) has a resonant behavior as the derivative 
d5\ (t) I dt has a pronounced maximum to be not caused by an opening of some new threshold. 
Since A\(t) is an analytic function, the rapid variation of 8\{i) can be described by a pole in 
the second Riemann sheet, which corresponds just to the well known pillO) resonance. Farther, 
starting only from general considerations we show that the latter pole has to appear also on the 
second Riemann sheet of the pion EM FF. 

Really, starting from the unitarity condition (3.30), analogically to the pion EM FF, one can 
prove that the branch point 4m 2 of the amplitude A\(t) is a square root type and also one can 
find an expression of the amplitude A\ on the second Riemann sheet to be 

By a comparison of (3.32) and (3.26) we see they have an identical denominator. So, if 
[^4}(£)] /7 has the p(770)-meson pole created by a zero of the denominator 1 + 2iA\(t), then the 
same zero has to be present in the denominator of [F^lt)] 11 , and as a consequence, the p(770)- 
meson pole has to be created on the second Riemann sheet of the pion EM FF. 

Since such a pole is shifted (due to the instability of p(770)) from the real axis into the 
complex plane, as a result of the reality condition (3.14), there must exist also the complex 
conjugate pole to the latter. 

We note, there are radial excitations [12] of the p(770)-meson. The corresponding poles are 
also placed on the unphysical sheets of the Riemann surface on which the pion EM FF is defined. 
Then all information on the analytic properties of the pion EM FF looks like it is presented in 
Fig. 3.2. 

Defining a phase of the pion EM FF for t > 4m 2 by the relation 

F®(t) = |J^(t)|e i4 'W (3.33) 
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Figure 3.2. Analytic properties of the pion EM FF 



and substituting the latter, together with (3.31), into the unitarity condition (3.18) one finds the 
identity 



6 n (t) = Sl(t) 
valid in the whole elastic region. 



(3.34) 



The next property of the pion EM FF follows from the threshold behavior of the P-wave 
isovector tttt- scattering partial wave amplitude 



Al(t)\^ ~ a\q 3 



(3.35) 



where a\ is the P- wave isovector nn- scattering length and q is the absolute value of the pion 
three momentum in the center of mass (cm.) system, which is connected with the momentum 
transfer squared by the following relation 



t = 4(q 2 +ml). 



(3.36) 



Taking into account (3.35) and the parametrization (3.31) for q — > one gets the threshold 
behavior of the phase S\ {t) q ->o ~ a\q 3 which by means of the (3.34) gives the threshold behavior 
of the pion EM FF phase 



*ir(*)|,-«, ~alq 3 . 



(3.37) 



One of the basic properties of the EM FFs of strongly interacting particles is their asymptotic 
behavior. However, until the discovery of the quark-gluon structure of hadrons the latter was 
unknown theoretically and only polynomial bounds were always assumed to be fulfilled. 

According to the quark model, the large momentum transfer behavior of the hadron EM FF 
is related to the number of constituent quarks of the considered hadron (2.24). Thus for the 
pion (n q = 2) one gets (3.1), which is in a qualitative agreement with existing data. This is a 
prediction of a dimensional counting rule and some simple assumptions which are based on the 
case of an underlying scale-invariant theory [3,4]. 

The asymptotic behavior (3.1) of the pion EM FF was proved [5-7] (up to the logarithmic 
correction) in the framework of the perturbative QCD. Here the pion EM FF takes the following 
form [5] (see Fig. 3.3) 
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Figure 3.3. Representation of the pion FF in QCD. 
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Figure 3.4. Dominating Feynman diagrams in the lowest order according to a s (Q 2 ). 



F?(Q 2 ) = [ dx [ dy$ + (y, Q 2 )T(x, y, Q 2 )$(x, Q 2 ) (3.38) 
Jo Jo 

where <j>(x, Q 2 ) is a wave function of the pion, which represents a probability to find the quark to 
carry the fractional momentum x of the total pion momentum and T(x, y, Q 2 ) is the amplitude 
of an interaction of the virtual photon with quarks and gluons. For the latter one can write the 
following perturbative expansion 



T(x, y, Q 2 ) - a s (Q 2 )T B (x, y, Q 2 ) [l + a s (Q 2 )T 2 (x, y, Q 2 ) + ....] 



(3.39) 



where a s (Q 2 ) is the effective constant of the quark-gluon interactions defined (in the lowest 
order) by the expression 



a s (Q 2 ) 
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(ll-2/3n / )lnQ 2 /A 2 



(3.40) 



with n f to be a quark number flavor and A as a QCD scale parameter. 

In the lowest order according to a s (Q 2 ) there are contributing only two Feynman diagrams 
(see Fig. 3.4) into T(x, y,Q 2 ). 

By means of a calculation of the latter one gets an explicit form for the Born approximation 
Tb of the amplitude T(x, y, Q 2 ) to be 



T B (x,y,Q 2 ) 



16irC F 



Q 2 (l-x)(l-y) 
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where Cf = 3/4 is the so-called colour factor. 

Now one has to find $(x, Q 2 ). It satisfies the differential equation (see e.g. [25,26]) 



d<I>(x,Q 2 ) 



l duV qq ^ qq (u, x)$(u, Q 2 

Jo 



(3.41) 



in which 



k = 2/(3 In 



«-(Qo) 



a s (Q 2 )' 

and the kernel takes the form [5] as follows 



= 11 - 2/3n/ 



v q ^ qq {u, x ) = i/2C F + - 1 -) + e( u - x) + ^(i + - 1 -) + e(x - «)} . 

(3.42) 

If the function $(x, Q 2 ) is known at Q\, then by using the equation (3.41) one can calculate 
it for arbitrary value Q 2 . In principle Q§ can be chosen to be arbitrary as soon as the condition 



<*a{Qo) 



< 1 



47T 

is fulfilled. 

The equation (3.41) can be solved easily if the property [5] 

f dzC^ 2 \2z l)V qq ^ qq {z, x) = l/2A% s C 3 J 2 (x) 
Jo 

is used, where 



(3.43) 



(3.44) 



A N J = C F 



n+l 



2 (n + l)(n + 2) — j 



(3.45) 



are the so-called non-singlet anomalous dimensions and C?/ 2 (x) are the Gegenbauer polynomi- 
als of the order 3/2. 

If one decomposes Q 2 ) into a series according to the polynomials C^ 3 ^ (x) as follows 



*{x, Q 2 ) = x(l -x)J2 <S>n(Q 2 )Cf/ 2 \2x - 1) 



(3.46) 



n=0 



and then substitute it into (3.41), one gets the equation for the coefficients of the previous expan- 
sion in the form 



d<I>n{Q 2 ) _ . NS 



dn 



A N n ^ n {Q 2 ). 



(3.47) 
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By means of a solution of the latter one obtains 

' a s {Q 2 ) 

<*s{Ql) 



$„(Q 2 ) = exp{KA» b )<5> n {Qi) 



<MQo) 



where 

d n = -2A» s /p . 

Now collecting all partial results we come to the expression of the pion EM FF 



from which the relation 



n=0 



a s (Q 2 ) 
as(Ql) 



4,a s (Q 2 



3' Q 2 



$o(OS)r + 2i2e (MQoWAQt,)) 



a s (Q 2 ) 

MQ 2 o) 



(3.48) 



(3.49) 



(3.50) 



(3.51) 



follows with di=0.427 for n/=4. 

If one restricts in the previous result only to the first leading term, then 

F*(Q 2 ) = M 4 3 ) a ^l\MQ 2 )\ 2 - 

In this case Farrar and Jackson have found [6] the relation 

$ (Qo) = ^fn, 



(3.52) 



where f n = 92.4 ± 0.2 MeV is the constant of the weak pion decay into [l- meson and antineu- 
trino. As a result one comes to the expression for the asymptotic behavior of the pion EM FF 



F?{Q 2 )cp- 
or by using (3.40) 



Q 2 



^ 2 fl 



(3.53) 



(3.54) 



(ll-2/3n / )g 2 lng 2 /A 

which confirms (up to the logarithmic correction) the quark counting rule prediction (3.1). 
Higher orders in (3.51) depend on the models chosen for the pion wave function. 

Till now we have considered only the invariant pion EM FF in the momentum representation 
to be sufficient for a description of observable phenomena. 

Nevertheless, there is a special Breit system in which the space-component of the pion EM 
current is equal to zero and as a consequence one can write down there the Fourier transform of 
the pion EM FF giving just the static charge distribution density as follows 



(3.55) 
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The inverse Fourier transform to (3.55) has the following form 

F* (Q 2 ) = J p{r)e-^d\. (3.56) 

If p(r) is a spherically symmetric distribution, then one can rewrite (3.56) into the spherical 
coordinates and by an integration over 9 and ip angles one gets 

F*{Q 2 )=^j~^p{r)r 2 dr. 

For the case of Qr <C 1 

{Qrf , 



sin Qr ~ Qr 
and 



F*(Q 2 ) = 4ir r 2 p{r)dr - H r 4 p(r)dr + .... 

Now, taking into account the charge distribution density normalization 

p(r)4:Trr 2 dr = 1 



10 

and a definition of the mean square charge radius 



/•OO 

(r 2 ) = / r 2 p(r)4irr 2 dr 
Jo 



one gets finally 

F 7 f(Q 2 ) = l-^<r 2 ) + ... 
from where the well known rule for a calculation of the mean square pion charge radius 

{rl)=6 dFgi± (35?) 

a( - lt=o 

is obtained. 

From the relations (3.11), (3.12) and (3.57) it is transparent to see that the mean square charge 
radius of the tt + - meson takes a positive value, the mean square charge radius of the 7r _ - meson 
takes a negative value and the mean square charge radius of the ir°- meson is identically to be 
zero. 
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3.2 Experimental information on the absolute value of charged pions form factor 

The pion EM FF is the simplest one of all others EM FFs of hadrons, the absolute value of which 
can be measured (directly or indirectly) everywhere on the real axis of the complex i-plane, on 
which the pion EM FF is defined, by using the following five different processes 

• the electroproduction of pions on nucleons (e~N — > e^nN); 

• the scattering of charged pions on atomic electrons; (n~e~ — > n~e~) 

• the inverse electroproduction processes (ir~p — > e + e~n); 

• the electron-positron annihilation into two charged pions; (e+e~ — > -k + it~) 

• J/^f decay into two charged pions (J/\t — > tt + tt~). 

Up to the present time there are more than 25 independent experiments realized in which an 
information on \F®{t)\ in the space-like (t < 0) and in the time-like (t > 0) regions, for the 
range of momenta -lOGcV 2 < t < lOGeV 2 , was obtained. 

The crucial (though not very precise) experiment was carried out in Novosibirsk [27], where 
the cross-section on the e + e~ — > ir + ir~ process was measured and an information on \F^(t) \ at 
the region of the p(770)- resonance was extracted. As a result, a creation of the p(770)- meson in 
the electron-positron annihilation into hadrons was confirmed experimentally for the first time. 
Later on this resonant region was re-measured with a higher precision in ORSAY [28] and the 
isospin violating w(782) — »■ tt + tt~ decay contribution to e + e~ — > tt + -k~ process, the so-called 
p — co interference effect, was experimentally revealed. 

Afterwards, during the next four decades many new experiments with the aim of obtaining 
of the experimental information on the pion EM FF in the space-like and time-like regions were 
proposed and realized. 

Since the pion is an unstable particle, one can not make a fixed pion target experiment of 
the elastic electron scattering on pions, like on the protons, in order to investigate the pion EM 
structure in the space-like region. Therefore, for very small values of t, the scattering of charged 
pions on atomic electrons, ir~ e~ — > ir~ er , was employed [29, 30] and for higher negative values 
of t the electro-production process, e~ N — »■ e~nN, was used [31-33]. 

In the time-like region, for < t < 4m 2 , an information on \F^{t)\ was obtained [34-36] 
from data on the inverse electro-production process -4- e~e + n. Also the process pn — > 
ir~£ + £~ is a suitable candidate [37] for the latter, provided that all corresponding nucleon EM 
FFs are known in this region. 

On the other hand, as the process e + e~ — > 7r + 7r~ is of the EM nature, one can treat the 
latter in the one-photon-exchange approximation and as a result, there are no model ingredients 
in an extraction of \F®{t)\ from the measured cross-section of the process under consideration. 
Therefore, the most reliable up to now information on the \F®{t)\ was obtained in a systematic 
investigations of the e + e~ — > 7r + 7r~ process in ORSAY [38-40], Frascati [41—43], Novosibirsk 
[44-46] and CERN [47]. 

The experimental point on \F^(t)\ at the highest value of t, t= 9.579 GeV 2 , was obtained 
[48] from the J/* — > tt + tt~ decay, which must be electromagnetic (i.e. it is realized through 
one virtual photon) provided that the G- parity is conserved by strong interactions. 
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Figure 3.5. Experimental data on \F^(t)\ 

As a result there are more than 300 experimental points on \F^(t)\ (see Fig. 3.5) for the 
range of momenta — lOGcV 2 < t < +10GeV 2 . However, they do not seem to be all mutually 
consistent and as a result too high values of the x 2 /NDF is in the fitting procedure by U&.A pion 
EM FF models found. 



3.3 Unitary and analytic model of charged pions electromagnetic structure 

Further we present a construction of the model which reflects all known properties of F^(t), 
including also the threshold behavior of ImF® (t) and describes all existing data. It represents 
just the most accomplished pion electromagnetic structure model to be constructed up to now. 
However, first very briefly all known properties of the pion EM FF are reviewed. 

1. The function F% (t) is normalized (3.13) at t = 0. 
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2. (t) can be considered as a complex function of a complex variable t and it is analytic 
function in the whole complex i-plane besides the cut from Am 2 to +00. 

3. A discontinuity across the latter cut is given by the unitarity condition (3.16), in which the 
first three allowed intermediate states are \tt + tt~) \tt + tt~tt + 'k~) and \tt°cj). 

4. To every intermediate state \n) in the unitarity condition (3.16) a branch point corresponds 
at the value of t to be equal to the squared sum of masses of the corresponding particles, 
i.e. they are at t — (2m 7T ) 2 , (Am v ) 2 , (m„o + m w ) 2 , etc.. 

5. From the Hermitivity of J^ M i.e. 



the reality of F^(i) on the real axis from —00 to Am 2 , and the reality condition (3.14), as 
well as a relation (3.15) between values of F^(t) on the upper (t + ie) and lower (t — ie) 
boundary of the cut ( t > Am 2 ) follow. 

6. By using of (3.18) one can continue analytically F^(t) through the upper and lower 
boundary of the cut at the interval Am 2 < t < 16m 2 and simultaneously one can prove in 
a such way that the first branch point at t = Am 2 is a square root type. 

7. As a by-product of the latter one obtains the expression (3.26) for the pion EM FF on 
the second sheet of the Riemann surface on which F^(t) is completely defined as an 
unambiguous function of t. 

8. One can read out from (3.26) the analytic properties of F^(t) on the second Riemann 
sheet. As the analytic properties of A\ on the first (physical) Riemann sheet consist of the 
right-hand unitary cut for Am 2 < t < +00 and the left hand dynamical cut for —00 < t < 
0, the pion EM FF on the second Riemann sheet has the left hand cut for —00 < t < too. 
In order to obtain a correct reproduction of the pion EM FF phase, a contribution of the 
latter can not be neglected as it was done in many phenomenological models in the past. 

9. As the 7T7T scattering amplitude is dominated by the p(770) resonance to be placed in the 
form of the pole on the second sheet and has identical denominator with (3.26), 
then F^(t) has also the pole on the second Riemann sheet corresponding just to the p(770) 
resonance. All excited states of p(770) are also placed on some unphysical sheets of the 
constructed Riemann surface. 

10. As a consequence of the reality condition (3.14) to every resonance pole also a complex 
conjugate pole has to exist . 

11. The identity (3.34) is practically valid for Ami < 1 < IGeV 2 . 

12. Directly from (3.33) the imaginary part of the pion EM FF takes the form 




EM 



(3.58) 



ImF*{t) = \FZ{t)\sm5«{t) 



(3.59) 
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from where by using (3.37) one gets the threshold behavior of the pion EM FF imaginary 
part 

ImF*(t\ =0 = \F*(t)\a{q 3 (3.60) 
which can be transformed into the following three conditions 



\q = * I 3 = 



13. The 7T7T scattering length a\ in (3.35) can be calculated from the pion EM FF by means of 
the expression 

d :, ImF^(t) 

14. The asymptotic behavior (3.54) of the pion EM FF obtained in the framework of the pQCD 
is equivalent (up to the logarithmic correction) to the quark counting rule prediction (3.1). 

To our knowledge there are no other properties of F E (t) to be known for the time being. 

Starting from the general parametrization (3.2) one obtains for the pion EM FF the expression 



f?[wu)] x (3 - 63) 

(Wn - w p )(w N - w;)(w N - i/w p )(w N i/wp * 
(w -w p ){w -w* p ){w -i/w p ){w -i/w;) 1 f p ' 

v (w N - w v )(w N - w;){w N + w v )(w N + w*) 
v J^ p „ {W -W V ){W -W*)(W + W V ){W + w*) 1 f v ' 

to be defined on the four-sheeted Riemann surface reflecting all above mentioned properties, 
besides the left-hand cut contribution and the threshold behavior (3.60) of the ImF^(t). 

The left hand cut contribution manifests itself mainly about the first 4mJ threshold as it 
provides a possibility to achieve a correct ratio of ImF^ (t) to ReF^ (t) at the expression 

Tm F E (t) 

^ 5 ^-r^fW) (3 ' 64) 

and in a such way as a result of (3.34) also a reproduction of the experimental data on 5\ (t). 

From literature it is well known, that the contribution of a cut of an arbitrary analytic function 
can be in the framework of the Pade approximations represented by alternating zeros and poles 
on the place of that cut. 

This method was employed in [49, 50] to demonstrate explicitly that experimental data on 
\F®(t)\ and 8\{t) can be described consistently only by the pion EM FF model including the 
correct left-hand cut contribution. 
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Really, in [49] tan S{ (t) was approximated by a rational function in q- variable with a min- 
imal number of coefficients to be determined from a fit of data on Sl(t). Then by means of a 
dispersion integral a rational function for F® (t) was obtained in which one zero and one pole 
just on the place of the cut from the second Riemann sheet were found. The latter result clearly 
demonstrates that the left-hand cut contribution can not be neglected in any analytic pion EM FF 
model. 

The same result was confirmed in [50] by a different way. Here, first, contributions of the 
isospin violating cj(783) — > 7r + 7r~ decay from the data on (Jtot{ e+e ~ — >■ 7r + 7r~) were separated. 
Then the resultant data were combined with data on 8\{t) in order to get the data on ReF®{t) 
and ImF^it). The latter were described by various Pade-type approximants in ^-variable, in 
which repeatedly one zero and one pole were identified on the place of the left-hand cut from the 
second Riemann sheet. Moreover, the found positions are almost identical with those from [49]. 

By means of the conformal mapping (2.145) the left-hand cut from the second Riemann 
sheet is mapped into the interval < W < +1. So, on the base of the results obtained in [49,50] 
the left-hand cut contribution is in our model approximated by one normalized zero and one 
normalized pole as follows 

(W-W Z )(W N -W P ) 

(W N -W Z )(W-W P ) ^■ OD> 

where W z and W p are left to be free parameters. 

The most natural incorporation of (3.65) into our model of F^(t) could be as a multiplicative 
factor to the p(770)-meson term in (3.63) as the left-hand cut from the second Riemann sheet is 
near-by to the elastic region 4ra \ < t < lGeV 2 . 

However, just for the sake of a consistent incorporation of the threshold behavior of ImF^ (t), 
given by (3.60), we incorporate (3.65) as a common factor to all p-meson resonances contributing 
to F^(t) as follows 

K[w(t)] ~ \T^wl) (w N -w,)(w-w n ) x (3 ' 66) 



f (w N - w p )(w N - w;)(w N - i/w p )(w N - i/wp t 



\ (w-w p )(w-w;)(w-i/w p )(w-i/w;) v f p 



v (Wn - w v )(w N - w*){w N + w v )(w N + w:) f m 

(W- W V ){W- W*)(W + W V )(W + W*) [ f v 



v—p'^p" 

Then all three terms in the wave-brackets give a nonzero contribution to the ImF^it), though 
p'- and //'-term (unlike the p-meson term) are pure real in the interval 4mJ < t < t in . Thus the 
left hand cut contribution in a consideration of the threshold behavior of ImF^ (t) in the USzA 
pion EM FF model is crucial. 

Now, calculating the ImF®(t) from (3.66) explicitly, one can immediately verify that - 
ImF^ {t)\ q=0 =0. A requirement of the first derivative of ImF®(t) at q = to be zero gives the 
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following condition on the coupling constant ratios 

l(W z - W p ) + 2W Z ■ W p ReW p (l + \W p \- 2 )]N p (f p7T J.f p ) + (3.67) 

v=p',p" \ VVv \ 

where 

n p = (w N - w p )(w N - w;)(w N - i/w p )(w N - i/w* p ) 

and 

N v = (W N - W V )(W N - W:)(W N + W V )(W N + W*). 

If condition (3.67) is fulfilled, then the second derivative of ImF^(t) at q = is automatically 
zero. The condition (3.67), together with the relation 

l/WA) = i ; (3-68) 

V — p.p' ,p" 

following from the normalization (3.13) and (3.66), give a system of two algebraic equations 
with three unknown coupling constant ratios. A solution of the latter can be written, e.g. in the 
form as follows 



(fp'TTTr/fp') 



N P" 

\w,,\* 



\w p ,\* 

\Wp»\ 4 



N„ 



(3.69) 



i f p" 1T-K / f p" ) 1 



\w p ,,\* 



N P" 

\w p „\* 



N„ 



, Wz.W p 
'Wz-W p 

\w p ,\* 

+ 



.Re [W p (l + \W P \- 2 ) 

N p// 
|W>| 4 



{fp-XK I fp) 



\w p ,\ 4 



N P" 

\w p „\* 



+ 



\w p „\ 



S(l + 2^-.Re\W p (l + \Wp\- 2 )])N p 



'W z -W t 



|W>I 4 



(3.70) 



■ (/p-n-Tr/ fp) 



Then the expression (3.66), together with relations (3.69) and (3.70), represents the most ac- 
complished USzA pion EM FF model, which is defined on four sheeted Riemann surface with 
complex conjugate poles (corresponding to unstable p-meson resonances) on unphysical sheets 
and reflecting all properties (including also the threshold behavior (3.60) of the ImF®{t)) briefly 
reviewed at the beginning of this paragraph. It depends on 10 physically interpretable free param- 
eters, ti n , m p> T p , (f p ,mr/fp), iTi p >, T p i , nri p " , T p n , W z and Wp. They are determined from the 
fit of existing reliable experimental points on | (t)\, containing, however, also a contribution of 
the isospin violating w(783) — > tt + tt~ decay, leading to the so-called p — uj interference effect, 
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Figure 3.6. Description of pion FF data by three resonance UhA model 



which can not be excluded by experimentalists in a measurement of a tot (e + e 
order to take into account the latter effect we have carried out the fit of data by 



\F^[W(t)]+Re 



icf> 



t - im u ,Y u 



where 



arctan 



m 



p 1 - p 



mf. 



-). In 



(3.71) 



(3.72) 



is the p — uj interference phase and R is the corresponding amplitude to be left as an additional, 
eleventh, free parameter. 

A description of all existing data by means of the most accomplished up to now UhA pion 
EM FF model (3.66) with (3.69) and (3.70), and the found values of parameters 



t 



(1.296 ±0.011)GeV 2 
0.3722 ± 0.0008 
TOp = (759.26 ± 0.04) MeV 
in pi = (1395.9 ± 54.3) MeV 
m p » = (1711.5 ±63.6) MeV 



R = 0.0123 ± 0.0032 
W p = 0.5518 ± 0.0003 
T p = (141.90 ±0.13) MeV 
T p , = (490.9 ± 118.8) MeV 
f> = (369.5 ± 112.7) MeV 



(/pW/p) - 1-0063 ± 0.0024 X7NDF = 1.58 



is presented in Fig. 3.6. A prediction of S^it) behavior by this UhA model and its comparison 
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Figure 3.7. P-wave isovector tttt phase shift. 

with data on 5{(t) is shown in Fig. 3.7. 

Just the coincidence of the predicted behavior of 8^ (t) with existing data on 5\ (t) confirms 
that the U&A pion EM FF model (3.66) with (3.69) and (3.70) fulfills the unitarity condition 
(3.18), at least from 4m 2 to aproximatelly lGeV 2 . 



3.4 Prediction of P-wave isovector tttt phase shift and inelasticity above inelastic 



threshold from e + e 



process 



The unitary and analytic model (3.66) of F® \W{t)\ has one elastic cut 4mJ < s < +oo and one 



effective cut L 



!GeV z < t < +oo. 



The elastic unitarity condition (3.18) can be utilized for the analytic continuation of F® \W(t)] 
through the elastic cut on the second Riemann sheet and as a result one obtains (3.26), the FF on 
the II. sheet to be expressed by FF on the I. sheet and the tttt partial wave amplitude A\(t) on 
the I. sheet, from where one obtains 



[K E (t)] ] 



2i[Ff{t)] 1 



(3.73) 



to be valid in the whole complex i-plane. 

Now, substituting a standard parametrization of I=J=1 tttt scattering amplitude at the physical 
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Figure 3.8. P-wave isovector -kit phase shift predicted by U&A model 



region 

A\(t + ie) = 
into (3.73), one gets 



r]{{t + ie) e 2W i< t+fe >-l 



2i 



r)l(t + ie)c 



[FZ[W{t + is)]] 1 



{FZ{W(t + is)}} 1 
from where it is straightforward to find 



1 



Im 



8\{t + ie) = -arctg- 
2 Re 



[F*[W(t+ie)]] 1 
[FV[W(t+ie)]] u 
[FElWjt+ie)]] 1 
[Ff[W(t+ie)]] n 



(3.74) 



(3.75) 



Vi(t+ie) = 



[FflWjt + ie)]] 1 



[FZ[W{t + l s)t 



. (3.76) 



Substituting the UkA model (3.66) of the pion EM FF with (3.69), (3.70), (3.26) and (3.73) 
into (3.76) one predicts [51] behavior of 8\(t + ie) and r]\{t + ie) in a perfect agreement with 
existing data (see Figs. 3.8 and 3.9) also in the inelastic region, i.e. above 1 GeV 2 . 

The latter is clear demonstration of the analyticity as one of the powerful means in the ele- 
mentary particle physics phenomenology. 
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Figure 3.9. Predicted behavior of inelasticity of the P- wave 

3.5 Excited states of the /9(770)-meson 

Here we demonstrate another example of the utilization of the analyticity as the powerful means 
in the elementary particle physics phenomenology. 

From Fig. 3.5, where the data on the pion EM FF are collected, one can see dominating role 
of the p(770)-meson. Besides the latter one can notice in the data also another explicit resonance 
around the energy s — 2.9 GeV 2 . But it does not mean that there are no more other hidden 
p-resonances in the pion EM FF to be cowered in shadow of some other resonances and the 
background. In such case it is difficult to identify the resonance by the Breit-Wiegner form. Here 
the U SzA pion EM FF model, providing one analytic function in the whole interval of definition, 
is starting to be very suitable. So, the correct approach in an investigation of the latter problem 
is then to take the expression (3.66), first, with two resonances, then with three resonances etc. 
and to look always for minimal value of \ 2 in experimental data fit. 

Such program has been practically realized. Considering only two resonances in (3.66), 
X 2 /NDF = 539/279 was achieved and p(770) with //'(1700) were identified. 

However, if three resonances were taken into account in (3.66), % 2 /NDF = 382/276 was 
found and two excited states, //(1450) and //'(1700), were identified [52]. 

Finally, by a consideration of four resonances in (3.66), x 2 /NDF = 343/273 was reached and 
in addition to p'(1450) and //'(1700), also the third excited state ^'"(2150) of the p(770)-meson 
was revealed [53] (see Fig. 3.10). 

3.6 Unitary and analytic model for kaon electromagnetic structure 

Unlike the pion EM FF, little has been done theoretically for the kaon EM FFs {t) and 
F^a (t) up to now. The main reason was, first, the existence of a rather broad unphysical region 
< t < Am 2 K unattainable experimentally, in which the dominating resonances p(770), cj(783) 
and 0(1020) of the kaon EM FFs are spread out, and secondly, shortage of reliable and compact 
experimental information outside the internal < t < Am 2 K . 

For instance, there are no data on the neutral kaon in the space-like it < 0) region, and for the 
charge kaon only 25 experimental points in the very narrow interval of momenta —0.1 145 GeV 2 
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Figure 3.10. Description of pion FF data by four resonance U&zA model 



< t < -0.0175 GeV 2 exist [54,55]. Slightly better situation is in the time-like (t > Am 2 K ) 
region, where the data by means of the e+e~ — > K + K~ and e+e~ — > K^Ks processes have 
been obtained up to t = 4.5 GeV 2 . 

The most successful description of this experimental information on the kaon EM FFs was 
achieved [53,56] in the framework of the UhA model (3.2) to be applied to the isoscalar and 
isovector parts of the kaon EM FFs 



F s K {t) = 



F v K (t) 



+ 



l-V 2 
i _ v 2 



\- (Vn - V s )(Vn - V s *Wn - 1/Vs)(Vn - 1/V s *) ff 
2^ 7Y7 t7\7t7 , t7T7t7 TTTTTTTTTTTtT^ \JsKK 



(Vi 



N 



(V - V S )(V + V S )(V - 1/V S )(V + 1/V S ) 

v^)(v N - v;){v N + v^Wn + v;,) 



//.)+ 



(v vp){v v*)(v + v r )<y + v*) 

t2 \ 2 



{f(j>'KK/f<t>') 



1-W 2 

l-W 2 N 



(3.77) 



(Wn - W V )(W N - W* V ){W N - l/W v )(W N - i/w* v ) tf 
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+ 



(w N - w p „)(w N - w;„){w N + w p n){w N + w;„) 



(w - w p n){w - w;„){w + w p „){w 

which are related to the charge and neutral kaon EM FFs by the relations 



(fp"KK/fp") 



F* + {t)=Ff ( {t)+Fl{t) 



(3.78) 



FkoW = Fk(t) ~ Fl{t) (3.79) 

as it follows directly from (2.7) and (2.8). 

The expressions (3.77) are analytic functions in the whole i -plane besides two cuts on the 
positive real axis starting, in the case of Ffc (t) from ig = 9m 2 and t\ n and in the case of the 
Ffc(t) from tjj = 4m \ and t\ n , to be real on the whole negative real axis up to positive values 
to = 9m 2 and = 4mJ, respectively, automatically normalized (2.9) if 

YsUsKK/fs) = Y^UvKKlfv) = 1/2, (3.80) 

S V 

with ImFKt) starting from 9m \ up to +oo and ImF1(t) ^ starting from 4m \ up to 
+oo, with poles corresponding to vector-mesons placed in the complex conjugate pairs on the 
unphysical sheets of the Riemann surface, on which both EM FFs F^t) and F^t) are defined 
and governing the asymptotic behavior (2.24) with n q = 2 as predicted by the quark model of 
hadrons. 

The reproduction of all existing data on \F K +(t)\ and \F K o(t)\ by means of (3.77), (3.78), 
(3.79) and the free parameters of the model (the values of parameters m p , T p , m w , T w are fixed 
at the world averaged values) 

iL = VWn - 9)/9 = 2.2326[m w ] qV n = ^ - 4)/4 = 6.6721[m w ] 

(Lkk/U) = 0-14194 (fpKR/fp) = 0.5615 

m = 7.2815^] m p > = 10.3940[m T ] 

T = 0.03733[m T ] T p , = 1.6284[m 7r ] 

(Ukk/U) = 0.4002 (fp'KK/fp') = -0-3262 (3.81) 

mp = 11.8700[m„] m p „ = 13.5650[m„] 



Tp = 1.3834[m w ] T p „ = 3.4313[ 



ni- 



(fpKR/M = --04214 (./>Wi» = --02888 

X 2 /NDF = 181/166 

are presented in Fig. 3.11. 

From the obtained results one observes that the contribution of the p'"(2150) resonance to 
the e+e~ — > KK processes is favored prior to the p" (1700) one by the existing data in the 
charge and neutral kaon EM FFs. 

In the Review of Particle Properties [12] are besides w(783) meson also its higher excitation 
states, cj(1420) and w(1600), presented. Therefore it could be interesting to investigate their 
contributions to the e+e~ -> KK processes by means of the UhA model of the kaon Em FFs. 

However, the latter problem can be considered seriously only in the case that more abundant 
and more precise experimental information on the kaon EM FFs will appear. 
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3.7 Conserved-vector-current (CVC) hypothesis and the v e e — > M M° and 
t~ -> v T M~M a weak processes 

At the end of fifties, by Gerstein and Zeldowich [57], and independently by Feynman and Gell- 
Mann [58], the conserved-vector-current (CVC) hypothesis, as a theoretical ground for an ex- 
planation of an approximate numerical equality of the muon decay constant G M and the neutron 
decay vector constant Gy, was postulated in the framework of the V-A weak interaction theory. 

The latter hypothesis, besides others, provides a relation between a matrix element of the 
vector part of the weak charged hadronic current and a corresponding matrix element of the EM 
current to be taken between two pion states. As a result of the foregoing, a probability of the 
7T + -meson beta-decay 7r + — > 7r° + e + + v e was predicted [57,58] theoretically. Its agreement 
with experimental results was presented [59] as a demonstration of the general validity of the 
CVC hypothesis in the weak interaction theory. 

Further, based on the CVC hypothesis and U&zA model of the pseudoscalar meson EM 
FF's formulated above, we predict [60], [61] a behavior of cr to t(-Ev ab ) and da/dE l ^ h of the 
weak D c e~ — > M~M° processes for the first time and the ratio of total probabilities r(r~ — > 
v T M~ M°)/T e are determined and compared with experimental values and theoretical values of 
other estimates. 

In general, the differential cross-section of the weak & e e~ — > M~M° reactions in the cm. 
system is given by the following expression 

dT? ~ (2s p + l)(2s c + 1) 64^^ ^ 11 ' ( } 

where s > 4mJ is the cm. energy squared, fc c m - = [(s — 4m^.)/4] 1 / 2 is a length of a 3- 
dimensional momentum of produced pseudoscalar mesons, p c m - = y / s/4 is a length of a 3- 
dimensional neutrino-momentum (assuming m Vo = 0) and s p and s c are spins of the antineutrino 
and electron, respectively. The matrix element M in the lowest order of a perturbation expansion 
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Figure 3.12. Reduction of the Feynman diagram to contact one 



can be calculated from the Feynman diagram presented in Fig. 3. 12a (the F^ (s) is the weak FF 
of a charged W~ -boson transition (W~)* —> M~M°), which for s <C fn\r- is reduced to a 
contact diagram shown in Fig. 3.12b. 

In the standard electro- weak theory an effective Hamiltonian describing the D c e~ —> M~M° 
process is 



Hl = -^=[v e7ll (l+ l5 )e}jW + h.c. ) 



(3.83) 



where G = 1.1663 x 10 GeV is the Fermi constant of the weak interactions and J,, is a 
pseudoscalar meson weak current. Then the matrix element corresponding to the diagrams in 
Fig. 3.12 takes the form 



M = ^fMl» (1 + 75) e{p x ) (h - k 2 f F%(a) 
and as a result, the differential cross-section (3.82) can be calculated explicitly 



dJ? 128vr 2 



(3.84) 



(3.85) 



where /3m = v/l~ ^ m \i I s ^ s tr, e velocity of produced pions and ■& is the scattering angle in the 
cm. system. 

Now considering an explicit form of dJ?=sin d dti dtp and then integrating over angles d and 
ip in (3.85), one gets the total cross-section 



^tot( S ) = ^ S -/?^K( S )| 2 



(3.86) 



of the weak v c e~ — > M~ M° process. 

In order to predict the behavior of (3.85) or (3.86) as dependent on s, one is in need of a 
knowledge of the weak pion FF F^} (s) as a function of s. Since, there are neither data on 
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the latter, nor accomplished dynamical theory of strong interactions able to predict the F$ (s) 
behavior for s > Am 2 M , not even a phenomenological approach like it is in the case of the EM 
FF's, one has to use the CVC -hypothesis [57, 58] providing a relation of the weak F$(s) just 
with the pure isovector part of EM FF's of the pseudoscalar mesons. 

In order to derive the latter relations we start with the standard expression of the pseudoscalar 
meson weak current 

J™ = V^ + (3.87) 

where and are the vector and axial-vector, respectively. 

On the other hand, the EM current of hadrons is composed of the sum 

tf = MOm + OHx. (3-88) 

where ( J^)^ is a third component of an isotopic vector current J M ( Jf, J%, J^) and (Js)O jU is an 
isoscalar current. 

In the second-half of the fifties a very predictable postulation was introduced 

V 1 = Jf - iJ%, (3.89) 

i.e. that the charged weak vector current and the isovector part Jg of the EM current are 
components of the same isotopic vector J M . Since, strong interactions are invariant with respect 
to the isotopic SU(2) group, the isotopic vector current J M obeys the relation 

d» Jf = 0, (3.90) 

which directly results in the formalism of the conserved-vector-current (CVC) hypothesis 

0„V" = 0. (3.91) 

Further, in order to derive a relation between F^{s) and Ff/ -1 (s), we start with a commu- 
tation relation 

[T h J?]=ie ijk Jjt, (3.92) 

with Tj being an isospin operator. The relation (3.92) is fulfilled automatically if J M is trans- 
formed in the isospin space like a vector. Now defining 

T_ =Ti -iT 2 

and (3.93) 

= Jf -i 4, 

one can prove the following relation 

[T_, Jg] = J^ = V (3.94) 
by using expressions (3.93). 
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If we multiply (3.94) from the right-hand side by a state vector of M + pseudoscalar meson 
|M + ) and from the left-hand side by a state vector of M° pseudoscalar meson (M° | of the same 
isomultiplet, and simultaneously we use the relations 

T_\T,T 3 ) = ^(T + T 3 )(T-T 3 + 1) \T,T 3 - 1) (3.95) 
(T,T 3 \T_ = ^(T-T 3 )(T + T 3 + 1) (T,T 3 + 1\ 

we get 

(M°\V»\M + ) = (M°| [T_, J 3 M ] |M+) = (3.96) 
= y/(T-T 3 )(T + T 3 + l) (M+\Jg\M-) - y/(T + T 3 )(T - T 3 + 1) (Af°|J^|M°). 
Then from (3.96) for pions one gets the relations as follows 

(w \V»\Tr + ) = V2(tt+\ J^\n+) - V2(ir°\ J^tt ) , (3.97) 

in which the second term is zero due to the charge conjugation. Really, if Uq is a unitary charge 
conjugation operator, and 

U C \*°) = UcJgUc 1 = -4, (3.98) 

then 

(ttVsV) = (AUc'UcJ^UcItt ) = -(A4\A = 0- (3.99) 

So, finally one gets the relation 

(7r |V> + ) = v / 2>+|J 3 >+). (3.100) 

On the other hand, if we multiply (3.88) from the right-hand side and left-hand side by a state 
vector of tt + meson, we get 

(7r+|j£|7r+) = (7r+|J 3 > + ) + (T + l^k + >- (3-lOD 

Because strong interactions are invariant under the G-transformation, that is a combination of 
the charge conjugation and an isotopic rotation for an angle 180° around the second axis in the 
isospin space, the second term in (3.101) is equal to zero. Really, considering that 

U G \ir+) = [/ G J^" 1 = -j£ (3.102) 

one has 

(tt+I J£|7r+) = (^ + |C/ G 1 (7GJ^ G 1 C/Gk + > = -(n + \Jg\n+) = (3.103) 
and as a result, 

(7r+|j£|7r+) = (n + \Jg\n + ). (3.104) 
Comparison of (3.104) with (3.100) leads to 

(7r°|\^>+) = \/2> + |j£|7r + ). (3.105) 
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Then parametrizing the matrix elements in (3.105) through corresponding FF's one gets finally 
the relation between the pure isovector part of the pion EM FF and the weak pion FF in the 
following form 

F™(s) = V2F*' I = 1 (s). (3.106) 
For kaons the expression (3.96) gives the following relation 

{K°\V^\K+) = (K+\Jg\K+) - (K°\Jg\K°). (3.107) 
By using the relation (3.88) one can write down for the difference 

(K+\J£\K+) - (K°\J£\K°) = (3.108) 
= (K+\Jg\K+) - (K°\Jg\K°) + (K+\J£\K + ) - (K°\j£\K a ). 

Now applying the relations 



T+\T,T 3 ) = v /(T-T 3 )(r + r 3 + l) \T,T 3 + 1) (3.109) 
(T,T 3 |T + = v/(T + T 3 )(T-T 3 + l) (T,T 3 -1\ 

to kaons as follows 

T+\K°) = \K+); (K+\T + = (K°\ (3.110) 
one gets 

(K+\Jg\K+) = (K + \J%T+\K a ) = (K+\T + Jg\K°) = (K°\J%\K Q ) (3.111) 
which in (3.108) leads to 

(K+\J£\K+) - (K°\J£\K°) = (K+\Jg\K+) - (K°\Jg\K°). (3.112) 
A comparison of (3. 107) and (3.112) gives finally 

(K°\V^\K+) = (K + \J£\K+) - (K°\J£\K°). (3.113) 

Then parametrizing the matrix elements in (3.112) through corresponding FF's, we obtain the 
relation between the isovector part of the kaon EM FF's and the weak kaon FF in the following 
form 

F%(8) = F K+ ( S ) F K o(s) = F%(8) + F%(s) - F s K {s) + F%(s) = 2F£( S ). (3.114) 

Substituting the relations (3.106) and (3.114) into (3.85) and (3.86) and using the unitary and 
analytic model of the pion and kaon EM structure formulated in we get the differential and total 
cross-sections 

ctofee ~ di r~' r ° ) = ^/s-'i^w^ (His 
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Figure 3.13. Cross-section comparison of electromagnetic and weak process. 



a tot (u e e- Tr-n ) = ^-fi • s\F^=\s)\ 2 



and 



da(P e e' -> if-^ ) _ G 2 



32tt 2 



/^ S |F£( S )| 2 sin 2 tf 



(3.116) 



(3.117) 



crtot(^e- - • s\F%(s)\ 2 (3.118) 

of the weak annihilation processes D c e~ — > 7r~7r° and /> e~ — > K~K° respectively, for the 
first time [60,61]. The comparison of crtotl^ce - — > 7r~7r°) with (7 to t(e + e~ — > 7r + 7r~), and 
comparison of <Jtot{v c e~ — > K~K°) with cr t ot(e + e~ — > K + K~) and crtot(e + e~ — > KgK^) 
is presented in Fig. 3.13 and Fig. 3.14, respectively. 

However, the results in Fig. 3.13 and Fig. 3.14 are not very interesting for experimentalists 
because the weak annihilation process z/ e e~ — > M~M° can be in principle measured only in an 
interaction of antineutrino beams with atomic electrons in the laboratory system. 
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Figure 3.14. Cross-section comparison of electromagnetic and weak process. 



The corresponding total cross-sections in the laboratory system can be found by using the 
following relation 

s = ml + 2m e E x f° (3.119) 
which leads to the expressions 

o™(p e e- -+ tt-tt ) = ^-(ml + 2m c E^)Pl\Ff' l =\E^)\ 2 (3.120) 

and 

a[ a b t (D c e- -+ K-K°) = ^-{m 2 e + 2m c E^ h )(3 3 K \F^ (£' ab )| 2 , (3.121) 
respectively. 

It is also interesting to have a prediction of the energy distribution of the pseudoscalar mesons 
created in the final state of the weak v c eT — > M~M° process given by the differential cross- 
section da/dEffi in the laboratory system. With this aim we use in (3.115) and (3.117) the 
explicit form of di? = — d cos i? dip and integrate the corresponding expressions over the angle 
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ip. Then we insert 

dcos?? = - — "h d£ ;iab s«2m CJ B' ab ; (3.122) 

/ rem. ZT'c.m. ™ ZT'lab^ 

sin 2 ,? = 1- ( ^ M „ ° ... Wc f 2 M > , (3.123) 

where 



^ m « i? c c - m - « P cm « ^ ^25jf*>; and fc cm ' « ^(£> ab - M 0) ) 

(3.124) 

with the threshold energy )/2m e of the antineutrino beam to be Ei°J = 

76.7 GeV for pions and E^ = 962 GeV f 
the energy distributions of pions and kaons 

plab Ei(0) \ f CZ^lab o z?lab\2 



76.7 GeV for pions and = 962 GeV for kaons. As a result, the following expressions for 



da(D c e- -> n-n°) E™> - \ I _ (iff* - 2iff b ) 2 I E>1=1 lab 

8^ ^ £' ab J\ £lab(£lab J 1 * 



2 



(3.125) 



and 



d<r(P c e- -> / ^ ab -4°^ f\ _ {E^-2E^f \ ^ lab 2 

dE%> - 4n { E^ )\ E^(E^-E%)j l K[ " > l 



(3.126) 



are obtained, respectively. 



In the last few years a considerable improvement in the detection of small exclusive hadronic 
tau-decay processes has been achieved. As a result, also experimental values on the branching 
ratio BR{t~ — >■ v T M~M°) were obtained recently, where M means a pseudoscalar meson. 

In order to predict it theoretically, we start with a general expression for a probability of 

r~(ki) -> v T (k 2 )M~ (pi)M°(p 2 ) decay process 

, r , ,4 d 3 k 2 d S Pl d 3 P2 

dT = {27T) J 2m; S{h - fc2 - ^ - *») (o^)^ (2 ^3 2i;i (27r )3 2j , 2 - < 3 - 127 ) 
where the corresponding matrix element takes the following form 

M = GFC J! 6c u(k 2 )^(l+ l5 )u(k 1 )jW (3.128) 

and 9c is the Cabibbo angle. In the expression (3.128) the VF-boson exchange mechanism is 
assumed to be realized. However, since m 2 / m 2 ^ < < 1, all effects due to the VK-boson exchange 
can be neglected. Gf means the Fermi constant of weak interactions and is a matrix element 
of the weak charge vector current responsible for the W~ — > K~K° transition. 
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Then the absolute value squared of the matrix element (3.128) takes the form as follows 

\Mf = G2FC f 6c l^ + r^J?(J?y (3.129) 

with 

t-nv = 4[fci^fc 2i , + k\ v k 2jl - gn V {ki ■ k 2 ] (3.130) 

Tpv = titpvapSakip - 4m T [s Al fc 2 „ + s„fc 2M - g^(s ■ fe 2 )], 

where s is the r-lepton spin four- vector, s ■ k\ — and the bar (3.129) means an average over 
the spin of the r-lepton and a summation over the spin states of the neutrino v T . 
It is straightforward to rewrite (3.127) into the following three equivalent forms 

iT _ IMP^ 
64tt 6 m T 

1287r d m% 
\M\ 2 

dT = ^^ m rdxdy 7 (3.133) 

where E is the energy of the M~ pseudoscalar meson at the r-lepton rest reference frame, 
k 2 = (fci — k 2 ) 2 = m T (m T — 2E V ) is the effective mass squared of the M~M° system, 
x = k 2 /m 2 and y = 2E/m T . 

If we define Q — Am 2 M /m T , then 

R<x<l; and VR<y<l. (3.134) 

The matrix element of the weak vector current of the W~ — > M~M° transition can be 
parametrized in the following way 

Jj = (pi- P 2)^(k% (3.135) 

where Fj^ (k 2 ) is the weak FF of the corresponding pseudoscalar mesons. By using (3.135), 
(3.129)-(3.130) and an integration over the energy of the M~ pseudoscalar meson in (3.130) one 
gets the relation as follows 

dT = G llZZl dk2{ml - fc2)2(m ' + 2fc2)(1 ^) 3/2 ! F m (k 2 )\ 2 - 0.136) 

Now, if instead of the weak pseudoscalar meson FF in (3.136) the EM FFs by (3.106) and (3.114) 
are substituted, one gets the decay widths 

r(r" -+ v T n-n°) = jT dfc 2 (^-fc 2 ) 2 K+2fc 2 )(l-^) 3 / 2 |F 7 f=^(fc 2 )| 2 

(3.137) 
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r(r" -+ u T K-K°) = SZ * 2 (™?-^) 2 (-? + 2fc 2 )(l-^) 3/2 |^(fc 2 )| 2 , 

(3.138) 

respectively. By using in (3.137) the pion EM FF U&zA model (3.66) with the corresponding 
parameter values and in (3.138) the isovector part of kaon EM FF USzA model (3.77) with 
the corresponding parameter values (3.81), and simultaneously Gf = 0.2278 x 10~ 6 [m~ 2 ], 
sin 8c = 0.22 we finally get the decay widths 

r(r" -> v T ir-ir°) = 0.3898 x lO" 11 ^] (3.139) 

and 

r(r" -> v T R-K°) = 0.2574 x lO -13 ^], (3.140) 
respectively. Dividing the latter values by the total width of the r-lepton to be 

r*°* = 0.159 x lO" 10 ^] (3.141) 
one gets the corresponding branching ratios 

BR(t- i/ r 7T-7r°) = 24.52% 

and 

BR(t- -> v T K-R n ) = 0.16% 
respectively. They can be compared with the recent experimental values 

BR exp { T - -> v T ir-ir a ) = (25.30 ± 0.20)% 

and 

BR exp { T - -> v T K~K ) = (0.26 ±0.09)% 

obtained by ALEPH Collaboration [62] in CERN. There is very good agreement of our theoreti- 
cal predictions with experiment. 

3.8 Parameter differences of the charged and neutral p-meson family 

Isotopic spin is to a very good approximation the conserved quantum number in strong interac- 
tions. A breaking of the corresponding symmetry occurs as a consequence of EM interactions 
and the mass difference of the up and down quarks. Practically, it is demonstrated in nature by a 
splitting of hadrons into isomultiplets. 
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In this paragraph we are concerned with p-meson resonances. The cleanest determination 
of their parameters comes from the e + e~ annihilation and r-lepton decay. However, as it is 
declared by the Review of Particle Physics [12], experimental accuracy is not yet sufficient for 
unambiguous conclusions. The difference for p° and are presented in averaged to be m p o — 
m p ± = -0.7 ± O.SMeV and r> - T p ± = 0.3 ± 1.3MeV, respectively. 

Nowadays the situation is changed in two aspects. 

On one hand, new very accurate KLOE data [43] on the pion EM form factor (FF) at the 
energy range 0.35GeV 2 < t < 0.95GeV 2 were obtained by the radiative return method in 
Frascati. Also corrected CMD-2 [45] and SND [46] Novosibirsk e + e~ -> tt + tt~ data have 
appeared recently. 

On the other hand the weak pion FF accurate data [63] from the measurement of the t~ — > 
ir~-K°v T decay by Belle (KEK) experiment were published recently. 

The Uk,A model (3.66) of the pion EM FF, to be represented by one analytic function for 
— oo < t < +oo was elaborated, which is always successfully applied for a description of 
existing data on the pion EM FF from e+e~ -> 7r + 7r~ and due to the CVC hypothesis [60] 
equally well also for a description of existing data on the weak pion FF from t~ — > 7r~7r tv 
decay. 

As a result more sophisticated evaluation of a difference of the p-meson families parameters 
can be achieved [64]. 

For the pion EM FF there is almost continuous interval of 38 1 experimental points for 
-9.77GW 2 < t < 13.48GeV 2 , which all are described by the pion EM FF model (3.66) 
in the space-like and the time-like regions simultaneously. The most important from them are 
accurate KLOE data [43] at the energy range 0.35GeV 2 < t < 0.95GeV 2 obtained in Fras- 
cati by the radiative return method and also the corrected Novosibirsk CMD-2 data [45] at the 
range 0.36GeV 2 < t < 0.9409Ge^ 2 and SND data [46] at the range 0.1521Gey 2 < t < 
0.9409Gey 2 , which can influence the finite results substantially. They are supplemented at the 
interval -9. 77GeV 2 < t < 0.3364Get/ 2 and 0.9557Gey 2 < t < 13.48GeF 2 by other existing 
data. 

An application of the pion EM FF USzA model (3.66) with three lowest resonances, in order 
to take into account the fact that the mass of the r lepton in decay t~ — > 7r _ 7r°^ T allows to 
reach not more than the energy corresponding just to the second radially excited state p"(1700) 
meson, leads to the best description of all existing 38 1 experimental points with the values of the 
parameters presented in the second column of the Tab. 3.1. 

The weak decay t~ -4- ir~n°is T is, like the e + e~ — > tt + tt~ process, dominated by the p~ - 
meson family resonances and thus it can be used to extract information on the charged p-meson 
family properties. 

From the conservation of vector current (CVC) theorem it follows (3.106), i.e. the 7r~7r° 
mass spectrum in the t~ — > 7r _ 7r°^ T decay can be related to the total cross section of the 
e + e~ — > 7r + 7r~ process. As a result the same USzA pion EM FF model (3.66) can be applied 
to a description of the weak pion FF data, which can be drawn from the measured normalized 
invariant mass-spectrum. 

Though there are measurements of t~ — > 7r _ 7r°^ T decay to be carried out previously by 
ALEPH [65] and CLEO [66], here we are concentrated only on the high-statistics measurement 
[63] of the weak pion FF from t~ — > Tr~n°v T decay with the Belle detector at the KEK-B 
asymmetric-energy e+e~ collider as they are charged by the lowest total errors. 
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Table 3.1. The values of the fitting parameters for the fit of pion FF. The values are shown for two cases, 
the result of fitting e + e _ data to the U&A model of pion EM FF (second column), the result of fitting t~ 
data to the U&iA model of weak pion FF (third column). The differences of the values for both cases are 
presented in fourth column. 



Parameter 


o 

P 


+ 

P 




kn [Gev 2 ] 


1.3646 ± 0.0198 


1.2432 ± 0.0157 




W z 


0.1857 ± 0.0004 


0.4078 ± 0.0013 




W P 


0.2335 ± 0.0005 


0.6197 ±0.0007 




m p [ MeV] 


758.2260 ± 0.4620 


761.6000 ± 0.9520 


3.3740 ± 1.0582 


m p > [ MeV] 


1342.3060 ± 46.6200 


1373.8340 ± 11.3680 


31.5280 ±47.9860 


nip" [ MeV] 


1718.5000 ± 65.4360 


1766.8000 ± 52.3600 


48.3000 ± 83.8060 


T p [ MeV] 


144.5640 ± 0.7980 


139.9020 ± 0.4620 


-4.6620 ± 0.8502 


I> [ MeV] 


492.1700 ± 138.3760 


340.8720 ± 23.8420 


-151.2980 ± 140.4150 


T p „ [ MeV] 


489.5800 ± 16.9540 


414.7080 ± 119.4760 


-74.8720 ± 120.6729 


fpTvn/fp 


1.0009 ± 0.0001 


0.9998 ± 0.0002 




X 2 /NDF 


1.78 


1.96 





An application of the pion EM FF U&.A model (3.66) to the best description of 62 experi- 
mental points [63] on the weak pion FF leads to the values of the parameters also presented in 
Tab. 3.1, however, in the third column. 

There may be a question about the fact that the obtained p-meson family parameters are dif- 
fering from those presented at the Review of the Particle Physics [12], especially of the p(770)- 
and //(1450)-mesons. As it is well known, the resonance parameters depend on the parametriza- 
tion used in a fit of data. However, in this contribution it plays no crucial role as finally we are 
interested in a difference of the corresponding parameters. 

Moreover, in a determination of the parameters we exploit the UhA model of the pion FF, in 
which any resonance is defined as a pole on the unphysical sheets of the Riemann surface to be 
considered as the most sophisticated approach in a description of resonant states. 

The difference of the p-meson family resonance parameters is presented in Tab. 3.1 in the 
fourth column. 

On the base of these results one can declare that the masses of the neutral p-meson family are 
lower than the masses of the charged p-meson family. 
Their widths are just reversed. 

Considering the evaluated errors, one can confidently affirm only in the case of the p(770)- 
meson parameters and also in the case of the p'(1450)-meson widths, that for the charged and 
the neutral states they are different. 

For other p-meson family parameters one can say nothing definitely and one has to wait for 
even more precise experimental data. 
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4 Electromagnetic structure of 1/2+ octet baryons 

Experimentally we know eight baryons [p, n];[/l];[Z ,+ , S°, E~] and [S^S 1- ] of the spin 1/2 
and the positive parity, which according to the SU(3) classification belong to the same multiplet 
and are compound of 3 light (up, down and strange) quarks. 

After more than half of the century of the discovery of EM structure of hadrons, a knowledge 
about the EM structure of octet l/2 + baryons is unsatisfactory up to now. Almost all experimen- 
tal investigations are concentrated to the proton and mainly in the space-like region. Less precise 
information exists on the neutron. Concerning other members of the baryon 1 /2 + octet, there is 
only one experimental point on the total cross section of the A A production in e + e~ annihilation 
and an upper limit on the cross-section of the production at t = 5.693GeV 2 [67]. There is 
no experimental information on the electromagnetic structure of T,^ 1 and S-hyperons up to now. 

The electromagnetic structure of each 1/2+ baryon is completely described (see 2.1) by 
two scalar functions of one variable, commonly to be chosen in the form of the Sachs electric 
G_E_e(t) and magnetic Gmb(£) FFs. They are defined by the relation (2.1 1), which have no more 
the asymptotic behavior of VMD model (2.29) and according to the quark structure of baryons it 
takes the form 



GEB{t)\t\^oo — G M B{t)\t\^oo ~ ^2- (4-1) 

Inspired by a good experimental situation, many more or less successful phenomenological 
models have been constructed for a global description of the nucleon EM structure. However, 
missing are analogous attempts to predict a behavior of EM FFs of other baryons of the same 
octet. The main reason is that almost all constructed phenomenological models depend on some 
number of free parameters which, however, have to be fixed in a comparison with experimental 
data. 

Further we show how it is possible to go round this problem to some extent. 



4.1 Nucleon electromagnetic structure 

The EM structure of the nucleons (isodublet compound of the proton and neutron), as revealed 
first time in elastic electron-nucleon scattering more than half century ago, is completely de- 
scribed by four independent scalar functions of one variable called nucleon EM FFs. They de- 
pend on the squared momentum transfer t = — Q 2 of the virtual photon. 

Nucleon EM FFs can be chosen in a divers way, e.g. as the Dirac and Pauli, Ff(t), F™(t) and 
F$(t), Fg(t), FFs, or the Sachs electric and magnetic, GE P (t), G_E„(i) and GM P (t), GMn(t), 
FFs, or isoscalar and isovector Dirac and Pauli, F^ N (t), F" N (t) and F^ 2 (t), i^JvM' ^ s an ^ 
isoscalar and isovector electric and magnetic, G s EN (t), G v EN {t) and G s MN (t), G v MN {t), FFs, 
respectively. 

All these always four independent sets of four nucleon EM FFs are related by 



76 



Analyticity in a phenomenology of electro-weak structure of hadrons 



GEp{t) — 

GMpif) = 
GEn(t) = 



G s EN (t) + G EN (t) = F lp (t) + ^F 2p (t) = 

= [F? N (t)} + F? N (t)} + + FS N (t)]; 

= G s MN (t) + G v MN (t) = F lp (t) + F 2p (t) = 

= [F? N (t) + F" N (t)} + [Fi N (t) + F% N (t)] ; 

= G% N (t) - G v EN (t) = F ln (t) + ^—jF 2n (t) = 

= [F? N (t) - F? N (t)] + foV*) - *2n(% 

= G s MN (t) - G v MN (t) = F ln (t) + F 2n (t) = 

= [F? N (t)-F^(t)] + [Fi N (t)-F^ N (t)}, 



GMnif) 

and at the value t = normalized as follows 



(4.2) 



(0 G Ep (0) = 1; G Mp (0) = 1 + n P ; G En (0) = 0; G Mn (0) = fi n ; 

(ii) G S EN (0) = G EN (0) = G S MN (0) = \(l + ^ p + //„); G V MN (0) = 1(1 + Mp - M „) ; 

(Hi) F lp (0) = 1; F 2p (0) = Mp ; *i n (0) - 0; F 2 „(0) - M «; 

(it;) ^(0) = f^(0) = 1; F 2 V(0) = l( Mp + Mn); *3W(0) = - Mn), (4.3) 

where ^ p and \i n are the proton and neutron anomalous magnetic moments, respectively. 

The Dirac and Pauli FFs are naturally obtained (2.10) in a decomposition of the nucleon 
matrix element of the EM current into a maximum number of linearly independent covariants 
constructed from the four-momenta, 7-matrices and Dirac bispinors of nucleons 



(N\ J* M \N) = eu(p'){ ltl F 1N (t) 



2wln 



OjXv(v' - P)vF 2N (t)}u(p) 



(4.4) 



with m N to be the nucleon mass. 

On the other hand, the electric and magnetic FFs are very suitable in extracting experimental 
information on the nucleon EM structure from the measured cross sections 



dcr lab {e-N e~N) 



a 2 cos 2 (9/2) 



1 



4E 2 sin 4 (0/2)l + (^)sin 2 (#/2) 



'EN 



t r<i 

Ami, MN 



t 



1 t_ 

4m 2 . 



'Am 2 N 



G 2 Miv tan 2 (0/2), 



(4.5) 
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a = 1/137, E-the incident electron energy, and 



ot T{e + e- -+ NN) = ^IL[\G MN {t^ + 2 -^\G EN {t)\% p N = yjl-^- 

(4.6) 

or 



af TiPP e+e") = 7=[|G M p(t)| 2 + — ^|G £p (t)| 2 ], (4.7) 

3Pc.m.vi 1 

(Pc.m.-antiproton momentum in the cm. system) as there are no interference terms between 
them. 

In the Breit frame, the Sachs FFs give the distribution of charge and magnetization within 
the proton and neutron, respectively. From all four Sachs FFs the neutron electric FF plays a 
particular role. Though the total neutron charge is zero, there is a nonvanishing distribution of 
charge inside of the neutron, which leads to the nonvanishing neutron electric FF. 

The isoscalar and isovector Dirac and Pauli FFs are suitable for a construction of various 
phenomenological models of the nucleon EM structure, however with the correct asymptotic 
behavior 

^l S (*)|t|->oo = ^1 (*)|t|-KX> ~ ^2 (4-8) 



*2 (%|-oo = n (*)|t|-oo ~ ^ (4-9) 

in conformity with (2.1 1) and (4.1). 

The most attractive of them is the Vector Meson Dominance (VMD) picture in the framework 
of which FFs are simply saturated (see the paragraph 2.3) by a set of isoscalar and isovector vec- 
tor meson poles on the positive real axis. However, this turns out to be practically an insufficient 
approximation and in a more realistic description of the data (especially in the time-like region) 
instability of vector-mesons has to be taken into account and the contributions of continua, to 
be created by n-particle thresholds, like, e.g., 27r, 37r, KK, NN etc., together with the correct 
asymptotic behaviours (4. 8), (4. 9) and normalizations (4.3) have to be included. 

In recent years, abundant and very accurate data on the nucleon EM FFs appeared. Most of 
the references concerning the nucleon space-like data can be found in [68]. More recent precise 
measurements are presented in [69]- [77]. Besides the latter, there are also new data on the 
neutron electric FF from BATES [78], MAMI [79]- [82] and NIKHEF [83]. 

For the time-like region data see [84]- [92]. There is, in particular, the FENICE experiment 
in Frascati (Roma, Italy) measured, besides the proton EM FFs [91], the magnetic neutron FF in 
the time-like region [92] for the first time. There are also valuable results on the magnetic proton 
FF [87, 88] at higher energies measured in FERMILAB (Batavia, USA). 
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Figure 4. 1 . Compiled proton electric and magnetic FF data. 

All this stimulated recent dispersion theoretical analysis [93, 94] of the nucleon EM FF data 
in the space-like region and in the time-like region [95] as well. 

The latter works are an update and extension of historically the most competent nucleon FF 
analysis carried out by Hohler with collaborators [96]. However, the model does not allow one 
to describe all the time-like data consistently, while still giving good description of the data in 
the space-like region. 

In the next paragraph, we construct a ten-resonance U&cA model of the nucleon EM struc- 
ture [97], defined on the four-sheeted Riemann surface with canonical normalizations and QCD 
asymptotics, which provides a very effective framework for a superposition of complex conju- 
gate vector-meson pole pairs on unphysical sheets and continua contributions in nucleon EM FFs 
. The model contains, e.g., an explicit two-pion continuum contribution given by the unitary cut 
starting from t = 4mJ and automatically predicts the strong enhancement of the left wing of the 
p(770) resonance in the isovector spectral functions to be consistent with the results of [93,98]. 

Another result of the presented model is the prediction of parameters of the fourth excited 
state of the p(770) meson and the automatic prediction of isoscalar nucleon spectral function 
behaviors. At the same time, a description of all existing space-like and time-like nucleon EM 
FF data, including also FENICE (Frascati) results on the neutron from the e+e~ — > nh process, 
is achieved. 



4.2 Unitary and analytic model of nucleon electromagnetic structure 

There are more than 500 experimental points with errors (see Figs. 4.1, and 4.2) collected for 
qualified analyses to be carried out by the 10 resonance U&cA model of the nucleon EM FFs. 

The model will represent, as we have mentioned previously, a consistent unification of the 
following three fundamental features (besides other properties) of the nucleon EM FFs 
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Figure 4.2. Compiled neutron electric and magnetic FF data. 



1 . The experimental fact of creation of unstable vector-meson resonances in the e + e -annihi- 
lation processes into hadrons. 

2. The hypothetical analytic properties of the nucleon EM FFs on the first (physical) sheet of 
the Riemann surface, by means of which just the contributions of continua are taken into 
account. 

3. The asymptotic behavior of nucleon EM FFs (4.1) to be proved also in the framework of 
the QCD [99]. 

Here we would like to note that a further procedure will not mean any mathematically correct 
derivation of the USzA model, but only an (noncommutative) algorithm of its construction which 
is, however, generally valid also for any other strongly interacting particles. 

It has been practically manifested to be optimal to saturate the isoscalar and isovector Dirac 
and Pauli nucleon FFs by 5 isoscalars (u>, <f>, u>', ui", </>') and 5 isovectors (p, p' , p" , p'" , p""), re- 
spectively, in the form of the VMD parametrization (2.96) to be automatically normalized with 
required asymptotic behavior. For F( N (t) and F" N (t) one takes the formula (2.96) with n = 
5, m = 2, and for ^^(t) and f^iv W one takes the general solution (2.96) with n = 5, m = 3. 
As a result the following zero-width VMD expressions are obtained 

9 9 

; m /, m , 
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however, they are already automatically normalized (4.3) and they govern the asymptotics (4.8) 
and (4.9), respectively, as predicted by QCD up to the logarithmic corrections. 

Despite of the latter properties the model is unable to reproduce the existing experimental 
information properly and only its unitarization (see paragraph 2.7), i.e., inclusion of the contribu- 
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tions of continua and instability of vector-meson resonances, leads to a simultaneous description 
of the space-like and time-like data. 

It is well known that the unitarity condition requires the imaginary part of the nucleon EM 
FFs to be different from zero only above the lowest branch point to and, moreover, it just predicts 
its smoothly varying behavior (see e.g. [93,98]). 

The unitarization of the model (4.10)-(4.13) can be achieved by application of the following 
special non-linear transformations 
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(4.14) 



respectively, and a subsequent incorporation of the nonzero values of vector meson widths. 

Here tg = 9m 2 ,, tg = 4m 2 ., 4 ; 4' 4' 4 &K square-root branch points, as it is transparent 
from the inverse transformations to (4.14), e.g 



V(t) = i 



1/2 



1/2 



1/2 



1/2 



1/2 



(4.15) 



and similarly for W(t), U(t) and X(t). 

The interpretation of tg = 9m 2 and ig = 4m 2 is clear. They are the lowest branch points 
of isoscalar and isovector Dirac and Pauli nucleon FFs on the positive real axis, respectively, as 
in the isoscalar case the 3-pion states and in the isovector case the 2-pion states are the lowest 
intermediate mass states in the unitarity conditions of the corresponding FFs. 

However, as it follows just from the unitarity conditions of FFs, there is an infinite number 
of allowed higher mass intermediate states and as a result there is an infinite number of the 
corresponding branch points (and thus, an infinite number of branch cut contributions) in every 
of the considered nucleon FFs. 

Since, in principle, an infinite number of cuts cannot be taken into account in any theoretical 
scheme, we restrict ourselves in every isoscalar and isovector Dirac and Pauli FF to the two-cut 
approximation. The second one, an effective inelastic cut, in every isoscalar and isovector Dirac 
and Pauli FF is generated just by the square-root branch points 4'4>4>4' respectively. 
They are free parameters of the model and the data themselves, by a fitting procedure, will 
choose for them such numerical values that the contributions of the corresponding square-root 
cuts will be practically equivalent to the contributions of an infinite number of unitary branch 
cuts in every considered FF. 
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Some experts privately have been suggesting us to fix these square-root branch points at the 
NN threshold. However, it has been demonstrated practically that this can be done only in the 
case of isovector parts of Dirac and Pauli FFs, but in none of the cases of the isoscalar parts of 
the Dirac and Pauli FFs. The latter are found at much lower values than the NN threshold. This 
result indicates that between the lowest = 9mJ branch point and the NN threshold there 
is some allowed intermediate mass state in the unitarity condition generating an important cut 
contribution which cannot be neglected in a description of the nucleon e.m. structure. We know 
from other considerations that it is just the KK threshold. 

So, by the unitarization of (4.10)-(4.13), one gets one analytic function at the interval — oo < 
t + oo for every FF 
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where 
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with tg = 9m 2 ., tg = 4m 2 ., = if^ 1 = 4m 2 N and t^,tf^ effective square root branch points, 
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Figure 4.3. The predicted different behaviors of Ge p and Gjvfpin t < region dependent on the fact if 
Rosenbluth technique data (dashed line) or JLab proton polarization data (full line) are used in the analysis 



which however can not be fixed at two-nucleon treshold as in the isoscalar case there is a remark- 
able contribution of KK intermediate state in the unitarity condition. 

There are no data in the unphysical region < t < Am 2 N in order to determine the parameters 
of p, u, <j), p',u)', <p' , p", uj" by the fitting procedure. Therefore, they are taken from Review of 
Particle Physics [12]. The parameters of p'" are taken from [53] and the parameters of p"" are 
left to be free. Then we are finally in [V(i)] with 4 free parameters, in [W(t)] with 3 free 
parameters, in -F| [£/(£)] with 3 free parameters and in F% [X(t)] with 4 free parameters. All 
of them are determined in comparison of the model with the collected 512 experimental points 
on the nucleon EM FFs in space-like and time-like region simultaneously [97]. As a result, the 
fourth excited state of p-meson with the slightly lower mass m p »« = 1455 ± 53 MeV than 
in [100] and Y p im = 728 ± 42 MeV is determined. Its existence, however, has to be confirmed 
also in other processes and not only in the e + e~ — > NN process. 

The quality of the achieved descriptions is graphically presented in Figs. 4.3 and 4.4 by 
dashed lines. 

Such situation has existed up to the year 2000, when measurements of the ratio of the electric 
to magnetic proton form factors were carried out by a completely different method at Jefferson 
Lab. in Newport News (Virginia). 



4.3 JLab proton polarization data puzzle 



Starting from the year 2000 a large progress has been done [101-103] in obtaining of the ratio 
(see Fig. 4.5) 
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in the space-like (t = —Q 2 < 0) region by measuring simultaneously transverse 

P t = (-2) v/r(l + T)G Ep G Mp tan(0/2) (4.21) 
Jo 

and longitudinal 

P = ] v^(TT7)GL p tan 2 (0/2), (4.22) 

components of the recoil proton's polarization in the electron scattering plane of the polarization 
transfer process er-p — > e~p, where h is the electron beam helicity, I is the unpolarized cross- 
section excluding ouott and r = Q 2 /Am 2 . As one can see from Fig. 4.5, these ratio data are in 
strong disagreement with the data obtained by Rosenbluth technique. 
It could be understandable as follows. 

Due to the fact that G 2 Mp {t) in (4.5) is multiplied by — t/Am 2 factor, the measured cross- 
section with increased —t becomes dominant by G 2 Mp (t) part contribution and the extraction of 
G 2 Ep (t) is more and more difficult. As a result the extraction of Ge p at higher values of —t by 
Rosenbluth technique is not promising. 

We have carried out a test of this hypothesis [104] in the framework of the ten-resonance 
USzA model of nucleon EM structure [97], which is formulated in the language of isoscalar 
F' 2 (t) and isovector 2 (*) parts of the Dirac and Pauli FF's and comprises all known nucleon 
FF properties. 

First, we have carried out the analysis of all proton and neutron data obtained by Rosenbluth 
technique together with all proton and neutron data in time-like region. 

Then all \GE P {t)\ space-like data obtained by Rosenbluth technique were excluded and the 
new JLab proton polarization data on [i p G e p (Q 2 ) / 'G m p (Q 2 ) for 0.49 GeV 2 < Q 2 < 5.54 
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Figure 4.5. The JLab proton polarization data on ratio fi v G e p / 'G m p (circles) and the same ratio calculated 
from data on Ge p and Gm p obtained by Rosenbluth technique (quadrangles). 



GeV 2 were analyzed together with all electric proton time-like data and all space-like and time- 
like magnetic proton, as well as electric and magnetic neutron data. 

The results of the analysis are presented in Figs. 4.3 and 4.4 from where it is seen that almost 
nothing is changed in a description of Gm p (<), GetiW an d GMn(t) in both space-like and time- 
like regions, and also \GE P (t)\ in the time-like region. There is only a difference in behaviors 
of GEpif) in t < region dependent on the fact if old data obtained by Rosenbluth technique 
are used (dashed line) or the new JLab proton polarization data are analysed (full line). The new 
JLab proton polarization data require in GE P (t) an existence of the zero (diffraction minimum) 
around the value t = — 13GeV 2 of the momentum transfer squared, which could change the 
charge distribution behavior within proton. 

The proton charge distribution (assuming to be spherically symmetric) is an inverse Fourier 
transform of the proton electric FF 



from where 



1 



(27T) 



, e'^ r G Ep (Q 2 )d 3 Q 



(4.23) 



sin(Qr) 
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Substituting for the Ge p (Q 2 ) under the integral either the older behavior from Fig. 4.3 given 
by dashed line, or the new behavior with the zero (see Fig. 4.3 given by full line) following from 
the new JLab polarization data, one gets different charge distributions within the proton given 
in Fig. 4.6 by dashed and full lines, respectively. That all leads also to different mean square 
proton charge radii. The old proton charge radius takes the value (r 2 ) = 0.68fm 2 and the new 
one (r 2 p ) = 0.72fm 2 . 

In order to distinguish which of these behaviors of GE P (t) in the space-like (t < 0) region is 
correct, we suggest to employ also the new sum rule [105] (to be derived later in this review) 

FU-Q 2 ) + ^2 F U-Q 2 ) ~ (4-24) 
- F 2 n (-Q 2 )-^F 2 n {-Q 2 ) = 



ira 2 y dQ 2 dQ 2 J 1 

giving into a relation proton and neutron Dirac and Pauli FFs in the space-like region with a 
difference of the differential proton and neutron cross-sections describing Q 2 distribution in DIS. 

Evaluating Dirac and Pauli FFs on the left-hand side corresponding to the old (dashed line) 
and new (full line) space-like behavior of GE P {t) in Fig. 4.3 one predicts the corresponding 
behaviors of the difference of deep inelastic cross sections in Fig. 4.7. By a measurement of the 
latter the true t < behavior of G_E P (i) can be chosen. 
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Figure 4.7. A prediction of two different behaviors of the right-hand side in (4.24) following from two 
different behaviors of Ge p in Fig. 4.3. 



4.4 Prediction of a tot (e + e -> YY) behaviors 



It is well known, that according to SU(3) classification of hadrons there is l/2 + octet of baryons 
including nucleon doublet [n,p] together with 6 other hyperons [A ], and 
Though there is almost zero experimental information on the hyperon EM structure, by using 
special nine-resonance USzA models of EM FFs of all members of the 1 /2 + octet of baryons, 
the experimental information on the nucleon EM FFs and SU(3) symmetry, one can predict 
behaviors of all hyperon EM FFs and as a consequence also behaviors of 



<7tot(e+e- -+ YY) 



4ira 2 /3 Y 



\G MY (t)\ 2 + 



-'"V r ,,,|2 



(*)f 



(4.25) 



Practically, one has to start with a specific nine-resonance U&iA model of the EM FFs of 
l/2 + octet of baryons, unifying compatibly all known properties of FFs, which renders just p-, 
co- and (/(-meson coupling constant ratios as free parameters. 

For nucleons, these free parameters are evaluated numerically by a comparison of the nucleon 
U<kA model with existing nucleon FF data and then it is straightforward to find numerical values 

Of fpNN, fuiNN and f<pNN- 

On the other hand, the trace of SU(3) invariant Lagrangian for vector-meson-baryon-antibaryon 
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vertex 

Tr(Lv B B) = ^f F [B^BP-B^B%](V„)Z + 

+ ^fB^B^l (4.26) 

with u> — <f> mixing 

<fp = cos i3 — lji sin i? 

uj° = (f> s sin d + wi cos ■& (4.27) 

B,B and V baryon, antibaryon and vector-meson octuplet matrices, to omega-meson singlet, 
f F , f D and f s SU(3) coupling constants and i9 mixing angle, provides the following expressions 
for vector-meson-baryon coupling constants 

fpNN — -^{fia^ -fia) 

fbjNN = ^costf/^-^sin^/^-/^) (4.28) 



f(1.2) _ 
JwAA 


-^COS0/£ 2 + 


1 


sin 0/£ 2 


f(l>2) _ 
— 


-^sintf/^- 


1 

71 




f(1.2) 
■/pEE — 


— /l^2 






f(1.2) _ 




i 

71 


sin ^ 


f(1.2) _ 


^sin^ 2 + 


i 

71 




f (l,2) _ 


2 (/i?2 — /l^) 







(4.29) 



(4.30) 



-j= cos <?A,2 + ^-j= sin 0(3/^ + /ft) (4.3 1) 
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Figure 4.8. e + e -annihilation cross-sections to A and S + hyperons. 




2.0 6.0 10.0 14.0 t [deV a ] 

t [GeV s ] 



Figure 4.9. e + e -annihilation cross-sections to S and E hyperons. 



Then the solutions of the system of algebraic eqs. (4.28) according to /f 2 , fi2< /l 2 w i tn nu " 
merical values of J p nn, JuNN and f^NN, by means of the expressions (4.29)-(4.31), enable to 
predict all free vector-meson-hyperon coupling constant ratios in the EM FFs of hyperons [A ], 
[£+, E°, £~] and S~] and, as a result, also behaviours of the total cross-sections (4.25) (see 
Figs. 4.8-4.10). 
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5 Deuteron electromagnetic structure 

The deuteron is the simplest object of all the existing nuclei and its EM structure provides an 
illustration of continuity between nuclear and particle physics at the microscopic level. The 
latter is completely described by three scalar functions 2.1 to be normalized as follows 

G cd (0) = l; G Md (0) = — Md! G Qd (0) = m 2 d Q (5.1) 

TTlp 

where [i d = 0.857406±0.000001[^tv] is the magnetic moment and Q = 0.2860±0.0015[/m 2 ] 
the quadrupole moment of the deuteron d. 

Experimentally, the EM structure of the deuteron is practically measured in the elastic scat- 
tering of electrons on deuterons described by the cross-section (2.21), from which exploiting the 
Rosenbluth methods one obtains the data on the elastic deuteron structure functions (2.22) to be 
compiled in [106]. 

In principle, by analogy with mesons and baryons, there is another source of experimental 
information on the EM structure of the deuteron in the time-like region for t > Am 2 d by means 
of the e + e~ annihilation into deuteron-antideuteron pairs, which is described by the following 
cross-section 



a tot (e+e- ^dd) = ^*\3\G cd (t)\' , , 



\G Md (t)\ 2 + 6^2l G <w(*)l S 



(5.2) 



where /3d = [1 — -^-J 1 / 2 is a velocity of produced deuterons or antideuterons. However, there 
are no suitable e+e~ colliding beam accelerators giving such an experimental information up to 
now. 

A lot of work has been done in theoretical attempts, like the traditional nonrelativistic impulse 
approximation [106, 107], sometimes also augmented by meson-exchange currents and isobar 
contributions [108], as well as the relativistic impulse approximation [106, 109], the hybrid quark 
model [110] and the Skyrme model [111], to describe the existing data on A(t) and B(t) in the 
space-like region. But in the framework of these models there is no concept of the time-like 
region behavior of the nuclear EM FFs at all. 

In this section we demonstrate [112] that the universal approach of the USzA model of the EM 
structure of strongly interacting particles, elaborated in the Chapter II., is suitable also for a de- 
scription of existing data on the elastic deuteron EM structure functions A(t), B(t) in the space- 
like region and simultaneously also for a prediction of the deuteron EM FFs Gcd(t), GMdif) 
and GQd{t) in the time-like region for the first time. 



5.1 Unitary and analytic model of deuteron electromagnetic form factors 

Any model, describing the EM structure of deuteron, should satisfy general properties of the 
deuteron FFs. As deuteron is a spin 1 particle, the ratios of the deuteron EM FFs at large space- 
like and time-like momentum transferred squared hold [113] the relations 



G cd (t) : G Md it) : G Qd (t) = (1 - §77) : 2 : -1, 



(5.3) 
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Figure 5.1. The Feynman diagram of the deuteron EM vertex generating the lowest anomalous threshold in 
deuteron electromagnetic form factors. 



which together with the QCD [1 14], [1 15] predictions for the asymptotic behavior of the deuteron 
EM structure function A(t) 

m)] 1 ' 2 -^-^ (5-4) 

imply the asymptotic behaviors of all deuteron EM FFs in both space-like and time-like regions 
to be 

G C d(t) ~ £ m Vco 

G M d{t) ~ *jt|->-oo ( 5 ' 5 ) 

G Qd (t) ~ ^-co- 
in the time-like region above the lowest branch point the deuteron EM FFs are complex 
functions of t- variable, i.e. they have imaginary parts different from zero to be given explicitly 
by the unitarity conditions. In other words the deuteron EM FFs are analytic functions in the 
whole complex t-plane besides cuts on the positive real axis starting from the lowest threshold 
(see paragraph 2.6) 

t = 4m^ - (^~ m p -0 2 = L7298m 2 (5 6) 

P ml 

to be anomalous [116] and generated by the Feynman diagram in Fig. 5.1. Its position is calcu- 
lated from the dual diagram presented in Fig. 5.2 by methods of elementary geometry. 

In order to obtain U &A parametrization of deuteron EM FFs with the asymptotic behavior 
(5.5) one needs to utilize the modified VMD parametrization with correct asymptotic behavior 
and normalization 



n- 



in 



2 



m) ~ f ° ii£hV (5 ' 7) 

/ , \ / , i 9. ,\ T-im / o -A TT m / 9 9\ T-rm r o ,\ < u k 
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Figure 5.2. The dual diagram from which by means of the methods of elementary geometry the position of 
the lowest anomalous threshold is calculated. 



discussed in the paragraph 2.5, where F(0) is normalization of the FF, m^'s are masses of the 
vector mesons, a x s are ratios of coupling constants (f x dd/fx), n is the number of saturated 



vector mesons and F h (t) 



t m is demanded asymptotic behavior of the examined FF. 



In what follows we apply the procedure of the unitarization (see paragraph 2.7) of every pole 
term in all three deuteron EM FFs in the form of the modified VMD (5.7) parametrization, where 
to is the lowest threshold (5.6) and tj„ is a free parameter independent for every deuteron EM FF 
to be simulating contributions of all higher branch points effectively. 

As a result one gets the following U SzA parametrization of the FF 



Fh(t) = 



i-*(*) s 



+ LH ( x k) 



k—m-\-l 



F - ]T a k I f[ LH(X V ) + 

\ k=m+l I v = \ 

V TT irn(x) C{Xv) - c{Xk) 
h tt\ LH{Xv) c(x v )-c(x t) 



, (5.8) 



where LH(X V ) equals 
in case of: 

LH(X X ) - 



in case of: m 2 x — T 2 X /A < t in 

(X N - X X )(X N - X*){X N - l/X x )(X N - l/X*) 
(X(t) - X x ){X(t) - X*)(X(t) - l/X x )(X(t) - l/X*) 
in case of: m 2 x — T x /A > U n 

(X N — X X )(X N - X*)(X N + X X )(X N + X*) 



LH[Xx) (X(t) X x ){X(t) X*)(X(t) + X x ){X(t) + X*) 
and C{X X ) equals 

in case of: m x — T x /A < t in 



(X N - X X )(X N - X*)(X N - l/X x )(X N - l/X*) 
-(X X -1/X X )(X*-1/X*) 
in case of: ■m 2 x — T 2 x /A> t ln 

(X N — X X )(X N — X*)(X N + X X )(X N + X*) 



C(X X ) = 
in case 

C(X X ) 



-(X x l/X x )(X* l/X*) 
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The minimal number n of isoscalar vector mesons needed to describe deuteron EM FFs 
depends on the asymptotic behavior, as n > m should be valid for every deuteron EM FF. 
However, it also depends on the existence of a node in the FF behavior. It can be shown, that 
special solution of Eq. (5.8) for n — m 

/i _ x(t) 2 \ 2m m 

F h (t) = F ( — l[LH(X v ), (5.9) 

^ N ' v=l 

is nonzero for any real value of t. As there are nodes in the behaviors of deuteron EM FFs 
Gc(t) and Gm(£), according to (5.5) the minimal number of vector mesons has to be 7(=6+l). 
Moreover, the positions [117] of the nodes t oc ^ -0.7 GeV 2 and i M — —2.0 GeV 2 give us 
additional constraints for the fitting procedure. 

The masses m x and corresponding widths T x of the first 5 vector mesons can be fixed to 
the world average values of the Tight' vector mesons oj, <p, uj' , lo" , <f>' '. Another ones can be fixed 
to the parameters of the charmed vector meson J/ty and m x ,Tx of the last vector meson will 
remain to be free, as we expect it to correspond to yet unknown Tight' vector meson resonance 
<f>" following from the SU(3) classification to be valid also for excited states of vector mesons. 

Finally, this procedure leads us to the U&lA parametrizations of the deuteron EM FFs 



G C d(t) 



GMd{t) 



G Qd (t) = 



1 - W(t) 
l-W% 



2\ 10 



(1 - «C:J/* - a C:x 



+ CH{W J/9 )ac:j/* + CH(W X ) 



(5.10) 



i-n*) s 



m d 
m 



-Hd-a M:x ) I] LH(V v )+CH'(V x )a M 

v=u> ,4>,<*j, 

w" ,J/<a 

(m 2 d Q-a Q .. x ) 11 LH(U v )+CH'(U x )a Q: 



1 A 1 ,J/* 



where 



CM{X W ) 
CH'(X W ) 



LH(X W ) ]T 
LH(X W ) ]T 



n [lh(x v ) 



n 



i=u>,0,u', v=u,4>,u',w", 
oj" ,4,' ,.J/<S 4>>,. I/<S, v^ti 



LH(X V ) 



C(X V ) - C(X W ) 
C(X V ) - C(Xi) 

C(X V ) - C(X W ) 



C(X V ) - C{Xi) 



(5T1) 
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with two conditions for ratios of coupling constants arising from the existence of nodes 

(1 - ac-. x ) IU^,^ LH(W v (t QC )) + CH(W x (t QC ))ac: X 
aC:J/ * " Y[ v = M ^ t¥ LH{W v {t Q c))-mWj,^c)) 

^Mdrit,= w ,0, W '^"^',j/*- L - ff ( i ^(*oM)) 

GM:X = T\7=^>,u"^,j/* LH(V v (t QM )) CH(V x (t 0M )) • 

The function X(t) takes different form for each deuteron FF, W(t) for Gc(t), V(i) for 
G M (t) and U(t) for G Q (t) 




5.2 Analysis of data on deuteron electromagnetic structure 



In order to have experimental information on all three deuteron FFs, one needs another observ- 
able, usually the component t 2 o of the tensor polarization [106] of the recoil deuteron, which 
contains the following combination of all three FFs 



£20 — — " 



1 



f2R{t) 



-■qGcdGQd + 



^ 2 G^ + ^(l + 2(l + J? )tan 2 §)G 



2 

Md 



(5.13) 



where R(t) = A{t) + B(t) tan 2 §. 

There are experimental data on the structure function A(t) [118]- [125] in the range —5.9536 
GeV 2 < t < -0.0016 GeV 2 , on the structure function B(t) [126]- [133] in the range -2.7556 
GeV 2 < t < -0.01 GeV 2 and on the t 20 [134]- [143] in the range —1.7161 GeV 2 < t < 
-0.0289 GeV 2 . 

The constructed USzA model of deuteron EM FFs depends on physical parameters like 

m v , T v for v = u>, <j>, w', to" , <j>', J/^, x, 

on the effective thresholds t in c,t in M, U n Q and on unknown ratios ac :x and aQ :x . 
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t [GeV 2 ] t [GeV 2 ] 

Figure 5.3. Behavior of the deuteron elastic structure functions A(t) and B(t) obtained by a simultaneous 
comparison of our U&iA model of the deuteron EM structure to existing experimental data on A(t), B(t) 
and polarization observable £20 



Naturally, we can fix masses m v and widths T v of all known vector mesons (w, <fi, w', u", <p' , 
J what means, that there will remain seven free parameters, which are numerically evaluated 
in the optimal description of all available experimental data on the deuteron elastic structure 
functions A(t), B(t) and the polarization observable t2o(t). The results were obtained by using 
CERN program ROOT [144] and corresponding behaviors are shown on Figs. 5.3, 5.4 and 5.5, 
together with data. The values of free parameters are given in Table 5.1. 

As one can see, we have obtained a quite reasonable description of the deuteron structure 
functions A(t), B{t) and polarization observable t 2 o(t), however with x 2 /n.d.f. = 4.61 indi- 
cating an inconsistency among independent sets of data. 

Moreover, the constructed model with the same values of free parameters can be used for the 
description of experimental data on the additional deuteron EM FFs Gcd{t) and GQd(t), as well 
as for estimation of all three deuteron EM FFs time-like region behavior ( see Fig. 5.6). They 
allow us to estimate the total cross section of the electron-positron annihilation into deuteron- 
antideuteron (5.2), which is planned to be measured on BES3 in Peking for the first time. 

Its predicted behavior is given in Fig. 5.7. 



Table 5.1. Fitted parameters of the U&A model of the deuteron EM structure. 





T^MeV 


t inC GeV 2 


tinMGeV 2 


t inQ GeV 2 


ac-.x 




504.9 ±0.1 


677.6 ±0.2 


18.2 ±0.1 


20.2 ±0.2 


7.8 ±0.1 


3.43 ±0.01 


28.14 ±0.03 
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t [GeV 2 ] 

Figure 5.4. Behavior of the deuteron tensor polarization observable t2o(t) obtained by a simultaneous 
comparison of our U&iA model of the deuteron EM structure to existing experimental data on A(t), B(t) 
and £20 




t [GeV 2 ] t [GeV 2 ] 

Figure 5.5. The predicted behavior of the space-like deuteron charge form factor Gc(t) and quadrupole 
form factor Gg(i) by our U&A model of the deuteron electromagnetic structure and their comparison to 
existing experimental data. 
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t [GeV 2 ] 

Figure 5.6. The predicted time-like behavior of the absolute values of the deuteron complex EM form factors 
Gc(t), GM(t), Gq{€) by our U&A model of the deuteron electromagnetic structure. Full line corresponds 
to the charge FF | Gc (t) | , dashed line corresponds to the magnetic FF | Gm (t) | and dotted line corresponds 
to the quadrupole FF j Gq (t) \ . 




t [GeV 2 ] 

Figure 5.7. The total cross section behavior of the e~e + — > dd annihilation process to be estimated by our 
U&A model of the deuteron EM structure. 

5.3 Study of two-photon contribution relevance into deuteron electromagnetic structure 

by impulse approximation 

As a result of a consideration of one-photon exchange approximation in elastic scattering of 
electrons on any hadron with nonzero spin, one obtains the same expression in the form for the 
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differential cross-section in the laboratory system to be expressed through the elastic structure 
functions A(t) and B(t). For comparison see (2.17) and (2.21). 

If also two-photon contributions are considered, no more simple form of the differential 
cross-section of elastic scattering of electrons on hadrons with nonzero spin is obtained. 

In the previous chapter we have analysed the new JLab proton polarization data [101-103] in 
the framework of the ten-resonance U&A model [97] of nucleon EM structure. The parameters, 
in comparison with those in [97] have been found to be changed very little. But a reasonable 
description of the new JLab data was achieved, whereby almost nothing has been changed in 
the description of GM P (t), Ge„(£), GMn(t) in both the space-like and time-like regions and 
\GEp(t) | in the time-like region. However, the existence of the zero (see the full line in Fig. 4.3), 
i.e. a diffraction minimum in the space-like region of \GE P (t)\ at t = —Q 2 rts —13 GeV 2 is 
predicted by such ten-resonance unitary and analytic model and so, \GE P (t) \ has no more dipole 
behavior in the space-like region. 

Recently it was suggested [145]- [147] that the two-photon corrections could be responsible 
for the found non-dipole behavior of \GE P (t) \ in the space-like region. Then there is a natural 
question arisen, if two-photon corrections play so important role in the elastic electron-proton 
scattering, what about a size of two-photon corrections in elastic electron-deuteron scattering. 

In this paragraph [148] the non-relativistic impulse approximation (NIA), which requires a 
knowledge only of the deuteron wave functions and the nucleon EM FFs, is used to study the 
previous question. 

As deuteron can be found in S- (« 96%) and D-state (« 4%), then NN non-relativistic full 
wave function of the deuteron can be written in terms of two scalar wave functions 

^ ' abm = ~ ^,m-ro s ( r )X a b " 

(I, l,m — m s ,m s \l, to) 
= ^y ,o(?)Xar + (5-14) 

+ ~~ H Y 2 ,m-m s (rkab' 3 (2,1, TO TO S , TO S 1 1 , to) , 

where {I, l,m— m s ,m s \l,m) are Clebsh-Gordan coefficients, Y^ mi are spherical harmonics 
normalized to unity on the unit sphere, XaT* * s tne s P m P 31 * °f tne wave function and zq = u, 
z 2 = w are reduced S— and D-state wave functions, respectively. 
The normalization condition 

Jd 3 r^ abm ,^ abm =<Wm 

implies normalization 

dr [u 2 (r) + w 2 (r)] = 1, (5.15) 

which could be understood as the probability of finding deuteron in S— or D-state. The D-state 
probability 

poo 

P D = drw 2 (r) 
Jo 
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Figure 5.8. The S and D state deuteron wave functions (u,w) behaviors for Paris potential [149]. 



is an interesting measurement of the strength of the tensor component of the NN force. 

The non-relativistic wave functions are calculated from the Schrodinger equation using po- 
tentials adjusted to fit NN scattering data for laboratory energies from to 350 MeV. In this paper 
we will use one of the most common potentials called Paris potential [149], which depends on 
the minimal number of free parameters and it was among the first potentials to be determined 
from such realistic fit. The S— and state wave functions determined from this model are 
presented in Fig. 5.8. 

The deuteron is a pure iso-scalar target, therefore within NIA its FFs depend only on the 
iso-scalar nucleon form factors G S EN and G S MN 



Gen — Gep + Geu 

Gmn = Gmp + Gmu (5.16) 
in the following way 

Gcd = G S EN D C 

GmcI = -z— [G s mn Dm + G s en De] (5.17) 
2m p 

GQd = G s en Dq, 



where the body form factors Dq, Dm, De and Dq are functions of the momentum transfer 
squared t. The non-relativistic formulas for the body form factors D involve overlaps of the 
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wave functions u(r), w(r), weighted by spherical Bessel functions 



D c (q 2 ) = / dr[u 2 (r)+w 2 (r)]jo{K) 
Jo 

/•OO 

D M (q 2 ) = / dr [2u 2 (r)-w 2 (r)] Jo ( K ) + 
Jo 

V2u(r)w(r) + w 2 {r) J2{k) 

drw 2 {r) [j'o(k) + ^(k)] 
w(r) 



(5.18) 



+ 



D E (q 2 ) 
D Q (q 2 ) 



3 



drw{r) 



u(r) 



V8 



where k = qr/2. At g 2 = 0, the body form factors become 



D c (0) 
D M (0) 
D e (0) 
Dq(0) 



= dr [u 2 (r) + w 2 (r)] = 1 

Jo 



dr [2u 2 (r) - w 2 (r)] 



'iP D 



f 

Jo 



drw 2 {r) 



drw(r) 



D 



(5.19) 



u(r) ~ 



w(r) 



giving the non-relativistic predictions 



Dq(0) 

» S N D M (0) + D E (0) = » S N (2 - 3P D ) + l.bP D , 



(5.20) 



where is the quadrupole moment of the deuteron, fi^ is the magnetic moment of the deuteron 
and /j, s n = |(/i p + Hn — 1) is the isoscalar nucleon magnetic moment. The experimental value 
of the deuteron magnetic moment fid = 1.7139, leads to the probability of D-state Pjj = 4.0%. 
But this is only approximate value, because the magnetic moment is very sensitive to relativistic 
corrections. 

Experimentally the EM structure of the deuteron is measured in the elastic scattering of 
electrons on deuterons, described by the differential cross-section (2.21) to be calculated in the 
one-photon-exchange approximation with the deuteron elastic structure functions. In order to 
see predicted behaviors of the deuteron elastic structure functions A(t) and B(t) to be caused by 
the non-dipole behavior of the proton electric FF with the zero around t = —Q 2 = — YiGeV 2 , 
we use the Gspif) (see full line in Fig. 4.3) together with all other nucleon FF behaviors in the 
calculation of the nucleon EM FF isoscalar parts by means of the relation (5.16). Then through 
the relations (2.22) and (5.17) one comes to the behavior of A(t) and B(t) as presented in Fig. 5.9 
and Fig. 5.10, respectively, by full lines. 

For comparison the deuteron elastic structure functions A(t),B(t) are predicted also by 
means of the nucleon EM FFs [97] obtained in the elastic scattering of unpolarized electrons 
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t [GeV 2 ] 



Figure 5.9. Deuteron structure function A(t) data and their comparison with NIA predictions using non- 
dipole (full line) and dipole (dashed line) Ce p (£) behavior. 




t [GeV 2 ] 

Figure 5.10. Deuteron structure function B(t) data and their comparison with NIA predictions using non- 
dipole (full line) and dipole (dashed line) GE P {t) behavior. 

on unpolarized protons by Rosenbluth technique, i.e. by using the proton electric FF in t < 
region having more or less dipole behavior (see dashed line in Fig. 4.3). The results are presented 
in Fig. 5.9 and Fig. 5.10 by dashed lines. As one can see from figures the examined difference 
is negligible, indicating that the two-photon-exchange contribution is unimportant for the unpo- 
larized electron-deuteron elastic scattering in comparison with the unpolarized electron-proton 
elastic scattering. 
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6 Sea-strange quark contributions to electromagnetic structure of hadrons 

Hadrons are classified according to the SU(3) symmetries into multiplets, the most of them are 
known nonet of pseudoscalar mesons and 1/2+ octet of baryons to be already treated in this 
review. The mesons in the quark model are qq bound states and members of octuplet are 3-quark 
qqq configurations plus any number of quark-anti quark pairs of light quarks u, d and s. 

At the nonet of pseudoscalar mesons the ir and -q are compound from only up and down 
quarks, whereas X-mesons and rf to some extent are compound also from strange quark. 

At the octet of baryons only proton and nucleon are compound from up and down quarks. 
The hyperons contain one or two strange quarks. 

This section we start with experimental indications on sea ss pairs contributions to the nu- 
cleon structure. 

They are as follows 

i) Pion-nucleon sigma term a^N 

The er-term of irN scattering is defined by the relation 

a nN {t = 0) = {p\(uu + dd)\p) (6.1) 

and it can be determined from extrapolation of scattering data on the isospin even nN on 
shell amplitude S = f 2 D{v, t) to the unphysical Cheng-Dashen point v = 0, t = 2fi 2 in 
the chiral limit m n = m u = nrid = 0. 

The results obtained from the analysis of irN data are now 

E(2 M 2 ) = a vN (2fi 2 ) = 60MeV (6.2) 
and consequently a nN (0) = 45MeV. 

On the other hand, the SU(3) mass splittings are due to the quark mass term H m (the 
isospin symmetry is assumed) in the QCD Hamiltonian 

H m = m(uu + dd) + m s ss (6.3) 

in which m s is the strange quark mass. The latter can be written in terms of singlet and 
octet contributions: 

H m = l/3(m s + 2rh)(uu + dd + ss) - l/3(m s - rh)({uu + dd - 2s s). (6.4) 

The octet matrix element 

6 = (N\(uu + dd - 2ss)\N) = (N\(uu + dd)\N) - 2(N\ss\N) (6.5) 

is determined from the octet masses by current algebra and SU(3) symmetry. As a result 
one finds 

S ~ 26MeV (6.6) 
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Since ovjv 7^ S, the scalar strange quark density (7V|ss|iV) is non-vanishing, leading to 
the following fraction of strange quarks in the proton 

Y=- \ N ^ N \ =0.21 (6.7) 
(N\uu + dd\N) 

ii) Proton spin crisis 

Deep inelastic scattering of polarized charged leptons (e~ or p,~) by polarized protons has 
been investigated at CERN and SLAC and it was found so-called EMC effect, i.e. quarks 
appeared to contribute very little to the proton spin. 

If A £ is the fraction of proton spin contributed by light quarks, then 

A£ = Aw + Ac?+ As, (6.8) 
where Aq's are integrals 

Aq = f dx{[q\x) - q\x)] - [q\x) - q+(x)]} (6.9) 
Jo 

of the difference of differences between parallel and antiparallel probability quark dis- 
tributions, q^-^(x) and q^(x) (q=u,d,s), for quarks and antiquarks with spin parallel and 
antiparallel to the proton spin and x is the Bjorken scaling variable. 

In the non-relativistic constituent quark model 

Ait = 4/3, Ad =-1/3, As = (6.10) 

and 

AS = 1. (6.11) 

However, the Am, Ac? and As can be determined from DIS spin-dependent proton struc- 
ture function g\{x) and baryon /3-decays. 

Really 

1. The integral of g\ (x) is a sum over Aq's, each weighted by the square of that flavor 
quark's electric charge 

Tf = J g[(x)dx = 1/2{4/9Au + 1/9Ag? + 1/9As}; (6.12) 

2. The matrix elements F and D measured in hyperon decays are related to A« - Ag! 
and Au + Ad - 2 As as follows 

g A = F + D =Au-Ad (6.13) 
3F-D = Au + Ad - 2As. (6.14) 



108 



Analyticity in a phenomenology of electro-weak structure of hadrons 



A recent analysis of T^, F and D finds 



Au = 0.83±0.03; 



Ad= -0.43±0.03; 



As = -0.10±0.03 (6.15) 



and 



AS = 0.30 ±0.05 



(6.16) 



indicating that very little of the proton spin is carried by the light quarks and that the 
strange sea is polarized antiparallel to the proton spin. 

iii) OZI rule violation. The OZI rule states that quark diagrams with disconnected quark lines 
are strongly suppressed. As a result the process pp — »■ 4>X has to be remarkably suppressed 
in comparison with the process pp — > coX. 

Recent experiments at LEAR at CERN have found violations of the OZI rule of up to two 
orders of the magnitude in pN — > <pX annihilations. 

Such data could be explained if there are ss pairs present in the nucleon (and antinucleon). 

iv) Neutrino experiments 

Elastic v and v scattering on protons is sensitive to strange axial-vector form factor, indi- 
cating on the presence of ss pairs in the proton. 

On the other hand, in deep inelastic scattering of v and v on protons the charmed particles 
can be produced in d — > c and s — > c transitions. The probability 



where # c =0.23 is the Cabibbo angle. 

Due to the smallness of 6 C the s — > c transition is dominant and gives the probability of 
studying the strange sea in N. 

All these evidences for strangeness in N mean that various elements of 

• scalar ss 

• pseudoscalar sj 5 s 

• vector s7 M s 

• axial-vector s-f^j^s 

• tensor sa^s 

currents, appearing in various processes with nucleons N, are non vanishing and it has sense 
to investigate them theoretically. 

Further we concentrate [150] only to the nucleon matrix element of the strange-quark vector 
current = s-f^s, which is experimentally accessible in parity- violating elastic and quasi- 
elastic electron scattering from the proton and light atomic nuclei, where the strange electric and 
the strange magnetic nucleon form factors (FFs) (or their combinations) are measured. 



d ->• c 



s — > c 
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6.1 Prediction of strange nucleon form factors behaviors 

The momentum dependence of the nucleon matrix element of the strange-quark vector current 
J*=S7 M s is contained in the strange Dirac F-f(t) and Pauli F^t) nucleon FFs 



(Nls^slN) = u(p ) 



«(p), (6.17) 



G' E (t) = J?(t) + G s M (t) = F?(t) + Fi(t). (6.18) 



by means of which the strange electric and strange magnetic nucleon FFs are defined 

t 

The latter can be measured in parity-violating elastic and quasi-elastic scattering of electrons on 
protons and light atomic nuclei. 

Since the strange-quark vector current = sj^s carries the quantum numbers of the 
isoscalar part of the EM current J^ =a and thus both currents couple to the nucleon through 
the same intermediate states, it is natural to expect [151] that one can extract the behavior of 
strange Dirac and Pauli nucleon FFs just from the isoscalar parts of Dirac and Pauli nucleon EM 
FFs. 

The main idea of a prediction of strange nucleon FFs behaviors from the known isoscalar 
parts of the Dirac and Pauli nucleon EM FFs is based on two assumptions 

• the u — (f) mixing is valid also for coupling constants between EM current (the strong 
strange quark current as well) and vector-meson 

11 1 11 1 

— = — cose -— sine; — = — sine +— cose, (6.19) 

Jul Jui a J<t>o J<t> J^a J<t>o 

where e = 3.7° is a deviation from the ideally mixing angle 8 = 35.3°; 

• the quark current of some flavor couples with universal strength n exclusively to the vector- 
meson wave function component of the same flavor 

{0\q r -fq r \(q t q t )v) = tm^ri^, (6.20) 

where my and are the mass and the polarization vector of the considered vector-meson. 

Starting from a definition of the virtual-photon vector-meson transition coupling constants 
1/fy by the relation 



<0|J^>=^e, 
Jv 



(6.21) 



and the second assumption for the isoscalar EM current J^ _0 to be expressed by quark fields, 
one comes to the equations 

<0|J^ =f Vo> = <Q\^{u 1 ^u + d 1 ^d)-^s 1 ^s\^{\uu> +\dd>) = 

11 1 m 2 

- s<7S + <6 - 22) 
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<0\J*-°\4>o> - <0\^(u ltl u + dj l td)-^s^s\ss>) = 

\ 777, 

= -3 Km l £ M = "77% ( 6 - 23 ) 
from where expressions for EM coupling constants follow 

1 1/ 1 1 n 11 1 1 1 /2 

- — = -(— — H =k=— -—=n: — — = — k = =W-k. (6.24) 

Substituting the latter into the first assumption, together with identities ^ = sin O and 
cos 0o, one obtains coupling constants of real co and (ft 

"7J = -^sin(£ + 6»o; = — ^= cos(e + O ). (6.25) 
These relations, together with — -^k following from 

<0|J^ =1 |p> = < 0|i(S7 M u~ d^d)\-^=(\uu > -\dd>) = 

11 1 TTV 2 

= 2 ( 7! + 7! )Km ' eAI = (6 ' 26) 

give for the ratios of the universal vector-meson coupling constants the values 
js : ■ 7V = 0.71 : 0.25 : (-0.32) 
in a very good agreement with experimental values 
h ■ 7^ : je = 0.79 : 0.23 : (-0.31) 

J p J <_u J 

obtained from leptonic widths r(V — > e + e~) of considered vector- mesons. Just this agree- 
ment demonstrates the previous two assumptions to be compatible with physical reality and one 
can extend their validity also for strong strange-quark current vector-meson transition coupling 
constants 1/fy. 

Then analogically one can write for the strong strange-quark current the equations 

< 0|J>„ >=< 0|(«7 M s|-^(|«u > +\dd >) = = -^e M (6.27) 

V 2 Juin 



< O|j;|0 o >=< 0|(s 7mS ||s S >= l.™; £ , = (6.28) 

ha 

from where one gets -X- — and -A- = I.k. Substituting them into u> — (ft mixing relations one 
comes to the strange coupling constants of the real to and (ft 



-jj = — Ksine; = +/«cose (6.29) 

J u J ch 
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Bringing these expressions for ui and <j> vector mesons into ratios with EM coupling constants 
(6.25), respectively, one gets rid of the unknown parameter n and comes to the relations 



(/wiV"iv//w) — ^o4«^ _i_ a ^ (fuNN/fu) 

sm(e + Wo J 

(fiL/ft) = -^ cos ( e ° S + £ g o) (/^iv//|) (i = l,2) (6.30) 



giving a possibility to calculate the unknown strange coupling constant ratios (parameters of 
strange nucleon FFs) from the known EM coupling constant ratios (parameters of EM nucleon 
FFs) to be determined in a description of all existing nucleon EM FF data by a suitable model of 
the EM structure of nucleons. 

The derived relations (6.30) are valid also for any pairs of excited states co 1 ',<f>' ', co", <p" ', etc. 
of the ground state of lo and </> isoscalar vector- mesons. 

Then for a prediction of strange nucleon form factors the U &A model of the nucleon EM 
structure [97] will be used, which comprises all known nucleon FF properties to be contained in 
the following models of isoscalar parts of Dirac and Pauli nucleon EM FFs 



- L(VMV„) ff"' - C »\ - L{V u .,)L{V u ,)]{fV N /r u ) + 
(G w " — L/ U i) 

+ [L(V^)L(V^) ^"~^\ - L(V^)L(V^) _ °^ 



(C<»)" — Cuj') (C«" — Cw') 

- L(K,»)L(K,0](4^//|)} (6.31) 



Fj=°[V(t)] = 



l-V 2 
i _ v 2 



L{V U „)L{V U ,)L{V U ) 



1 - 



{Cuj" — C w ) (Cui' C u ) 



Cu," 



+ L(V^)L(V^)L{V^) 



(Cuj" — C w <) 
(fuNN I ' fu:) + 



(Cw" — Cw') \ Co; 



(C u > — C$) 

Cu>" 



(^nn/ f<j>) 



(6.32) 
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and of the Dirac and Pauli strange nucleon FFs 



Fi[V(t)} 



where 



1 - V 2 
1 _ v 2 

1 V N 



L(V U ,)L(V U ) C "\ L{V„„)L{V U ,) 

L{V U ,,)L{V ^"~^\ - L{V U ,)L{V<,) 



{C^-Cu.) 



(/(iiVJv/ /«) + 

(Cy - C<f,) 
(C u » - C w /) 



- L{V u „)L{Vu,) 



1-V< 



1 - T/ 2 

1 N 



L{V„„)L{V U ,)L{V U ) 



(6.33) 



(6.34) 



(/^/.^) + L{V u „)L{V u ,)L{Vt) x 



Coj" 

( (C w n - C^) _ (Cq,/ - C^) 



Ojv - Vr)(v N - v;)(v N - i/Vr)(v N - i/v;) 

(V-V r )(V-V*)(V-l/V r )(V-l/V*) ' 

(Vn - Vr)(V N - V;)(Vn - IJVtKVn - i/v*) 
-(V r - l/V r )(V* - l/V*) 

V N = V(t)\ t=0 ;V r = V(t)\ t= ( mr -ir r /2)^,{r = u, </>,(*>', u"), 



(/<j>JVJv/ f<t>) 



L(V r ) = 

C T = 



V(t) = i 




(6.35) 



and t NN = Am 2 N is a square-root branch point corresponding to NN threshold. 

The expressions (6.31) and (6.32) for F{ =0 ,Fj =0 , respectively, together with similar ex- 
pressions for F/ =1 , F^ =1 [97] were used for a description of all solid nucleon EM FF data and 
a reasonable description of them has been achieved. The found numerical values of (f^ NN / f^), 
(f<j>NN I 'ftp i = 1, 2 in the relations (6.31), (6.32) are used to calculate the unknown strange 

coupling constant ratios {f^ NN / f^), {f$ NN /f$) i = 1, 2 in (6.33), (6.34) by means of the 
relations (6.30). As a result the behaviors of the strange nucleon FFs are predicted (see Fig. 6.1) 
without the use of any experimental point obtained in parity-violating elastic and quasi-elastic 
scattering of electrons on protons and light atomic nuclei. 

Due to the fact, that the expressions (6.33) and (6.34) for the strange nucleon FFs are parame- 
trized also by the UhA model and so, both FFs are analytic functions for — oo < t < +oo, 
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Figure 6.2. Predicted behaviors of the strange nucleon FFs in the time-like region 



naturally we are predicting the time-like behavior of strange nucleon FFs (see Fig. 6.2), though 
there is no method to be known for their experimental determination until now. 

Moreover, as the USzA model represents compatible unification of the pole and continuum 
(given by cuts on the positive real axis of the analytic strange nucleon FFs) contributions, one 
can predict even the imaginary parts behaviors of the strange nucleon FFs to be given by unitarity 
conditions of FFs under consideration. 

Since the first results reported by the SAMPLE Collaboration in 1997 [152], approximately 
ten independent measurements of the parity-violating contribution to the elastic EM FFs of the 
nucleons have now been completed. But only four of them [153]- [156] declare clearly nonzero 
experimental values of the strangeness within the proton and our theoretical predictions are com- 
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Figure 6.3. Compatibility of our theoretically predicted curve with GO Collab. data. 



patible with them. 

The A4 Collaboration [153] result at Q 2 = 0.2'iQGeV 2 is G% + 0.225G| f = 0.039 ± 0.034 
to be in accordance with our prediction 0.055; the SAMPLE Collaboration [154] result at Q 2 = 
O.lGeV 2 is G S M = 0.37 ± 0.34 to be also in accordance with our result 0.18; another A4 
Collaboration [156] result at Q 2 = Q.l08GeV 2 is G% + 0.225G S M = 0.039 ± 0.034 to be again 
in accordance with our theoretical prediction 0.030. 

May be the most impressive is a compatibility of our theoretically predicted curve (see 
Fig. 6.3) with the recent data [155] obtained by GO Collab. on the combination G S E (Q 2 ) + 
rj(Q 2 )G a M {Q 2 ) for the interval Q.\2GeV 2 < Q 2 < l.OGe^ 2 of momentum transfer squared 
values, which strengthen our belief in the strangeness in the nucleon. 



6.2 Strange vector form factor of if -mesons 

We have extended the method elaborated for prediction of the strange nucleon form factors be- 
haviors to Zf-mesons. 

Taking into account the splitting of FFs for charge kaon (3.78) and for neutral kaon (3.79) 
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Figure 6.4. Predicted behaviors of the strange kaon FFs in the space-like and the time-like region 

into isoscalar and isovector FFs and having the U&.A models of the kaon EM structure 

FfXt)} = (l^J) {l H (Vu>) + [HV u )-H(V u .)](f uKK /fZ)+ 

+ [L(^)-J?(^)](W//;)} (636) 

Fk =1 [W(t)} = ([^) {lH(W p ,,,) + [L(W p )-H(W p ,,,)}(f pKK /f;)+ 

+ [L(W p ,)-H(W p ,„)](f p , KK /r p ,)} (6.37) 
and also for the strange FF of kaon with the inner structure of F^"°[V(t)] 

F Kiv(t)} = f^T) {-H{v u ,) + [Hv u )-H<y u ,)]{f uKK /K)+ 



i — v 2 

1 V N 



+ [L(Vt)-H(V u ,)](U K K/f%)j (6.38) 

set up, one predicts [157] the behavior (see Fig. 6.4) of the strange FF of if-mesons (6.38) by 
means of an evaluation of {f uK K I '/£)> {f^Kx/f^) from {f uK K/fZ,)> (Ukk//^) determined 
in a comparison of (3.78) and (3.79) with data on charge and neutral kaon EM FFs. 

The obtained result is interesting from the point of view that in the case of the i-C -meson, 
which is naturally compound from one non-strange (up ore down) and one strange valence quark, 
the nonzero behavior of the strange vector FF of kaons presented in Fig. 6.4 means that besides 
explicit contribution of the strange valence quark into the kaon structure there are nonzero con- 
tributions to the if -meson structure also from sea strange quark ss pairs. However, we have no 
idea how to measure it experimentally. 
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7 Polarization phenomena in electromagnetic interactions of hadrons 

The polarization phenomena appear if beam particles or target particles, or both are considered 
to be polarized. If colliding particles are unpolarized, nevertheless final particles can be found in 
polarized states. 

An investigation of the polarization phenomena is very important in particle physics phe- 
nomenology as they provide more wealthy information on the EM interactions of hadrons. Some- 
times they reveal new unexpected results. Just typical example is JLab proton polarization data 
puzzle, discussed at 4.3 

Prior to the year 2000 all data on proton EM FFs GE P {t) and GM P {t) in the space-like 
(t < 0) region were obtained by measuring the differential cross-section of elastic scattering of 
unpolarized electrons on unpolarized protons in the laboratory system, utilizing the Rosenbluth 
technique. Both FFs manifested more or less dipole behaviors and their ratio in error bars is 
approximately equal one. 

Though polarization techniques have been suggested long time ago [158], only at the begin- 
ning of 21st century they were used [101-103] in obtaining of the ratio (see Fig. 4.5) 

which revealed a non-dipole behavior of the proton electric FF GE P (t) leading to the JLab proton 
polarization data puzzle. 

Further we present some other useful exploiting of polarization effects in particle physics 
phenomenology. 



7.1 Prediction of polarization observables in e+e -> pp process 

This paragraph is devoted to the analysis of polarization effects [159] in the process e+e~ — ► pp 
calculated in the framework of the one-photon exchange approximation. 

The above-mentioned process is interesting as it has noticeable polarization effects even if 
there are no polarized particles in the initial state. The appearance of such polarization effects is 
due to GE P (t) and GM P (t) being complex with non-zero relative phase. On that account there 
are also nontrivial polarization effects in the scattering of longitudinally polarized electrons on 
unpolarized target. 

The matrix element of the process e + e~ — > pp in the framework of the one-photon exchange 
approximation to be presented in Fig. 7.1 is defined by the formulae 

e 2 

M = p^ J M> 

U = «(-p)7a» u (p')> ( 7 - 2 ) 
J, = u{k)[F lp {t) 1 ^-F 2p {t)' J ^]u{-k% 

where t = q 2 > Am 2 . We note, that the cm. system of the reaction e+e~ — » pp is the most 
suitable for the analysis of polarization effects. 
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Figure 7.2. Definition of the scattering plane of the e e — > pp process 



The EM currents and J M are conserved q ■ j — q ■ J = and the matrix element M is 
completely determined by the product of spatial components of the currents j and J. 

The electromagnetic current J can be expressed through two-component spinors ipi and tp 2 



J = Vipf 
where we denote 



2m r 



GMp(t)(& — na ■ ft) H jH-G Ep{t)n& • n 



2m r 



F = Vt[GM P {t)(cr — na.n) H G E P {t)n a \n], 

yt 



(7.3) 



(7.4) 



a are Pauli matrices, n = (0, 0, 1) is the unit vector along the three momentum q of the proton, 
fh = (— sin^, 0, cos $) is the unit vector of incoming electron (see Fig. 7.2) and Gsp(i) and 
GMp(t) are defined by the relations (4.2). 

In order to find the corresponding cross-section, one has to calculate |Af| 2 . The differen- 
tial cross section of the reaction e + e~ — > pp in terms of EM FFs, for the case of unpolarized 
particles, has the form 



da 

dn 



a 
4s 



-\G Ep \ 2 sin 2 d + \G Mp \ 2 [l + cos 2 0] , a 



e 

4tt' 



t 



4m 2 ' 



(7.5) 



Further single and double spin polarization observables are calculated explicitly. 

First the single spin polarization observables (either for proton or for antiproton) in the case 



of 
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• unpolarized incoming leptons 

• incoming electron to be longitudinally polarized, 
are calculated. 

In the case of unpolarized initial leptons the corresponding lepton tensor takes the form 

jij = 2t(Sij - mirrij) (7.6) 
and the vector polarization is 

P = 1 \ . (7.7) 

jijTrlFiFj] 

Calculating the corresponding trace in the numerator, and similarly in the denominator, one 
obtains only P y component of the vector polarization to be nonzero 

P x = 

sm20 ■ Im[G* Mp (t)G Ep (t)} 
^{1/t\G Ep \ 2 sin 2 + |G Mp | 2 (l + cos 2 0)] 
P z = 0. 



F V - r=u 19.- 2„a , ,^ 19/-, , _.9.aM (l -*> 



The y-axis is orthogonal to the scattering plane (see Fig. 7.2) defined by the unit vectors m 
and ft, along the three -momentum of the electron and along the three-momentum of the created 
proton, respectively. 

The contributions of P x and P z in proton polarization are different from zero only if the 
electron (or positron) is longitudinally polarized, i.e. the lepton tensor takes the following form 

jij = 2t( Sij - mirrij + Xieijimi ). (7.9) 

Then the components of the vector polarization P of the created proton (or antiproton) in the 
• pp are 



2sm#-Re[G Ep (t)G* Mp (t)} 
p(*)I : 

sin2tf ■ Im[G* Mp (t)G Ep (t)} 



p = ' y ' " (7 10) 

V^[l/T\G Ep (t)\ 2 sin 2 e + \G Mp {t)\ 2 {l + cos 2 9)] 

P* = 



^[l/T\G Ep (t)\ 2 sin 2 i? + \G Mp (t)\ 2 (l + cos 2 0)] 

2costi\G Mp (t)\ 2 

[l/T\G Ep (t)\i sin 2 9 + |G Mp (t)| 2 (l + cos 2 6)] ' 

assuming 100% longitudinal polarization of one of the leptons. 

In a similar procedure one can find explicit forms of double spin polarization observables 
in the e+e~ — »■ pp process, where we are interested for polarizations of created proton and 
antiproton simultaneously. The corresponding polarization tensor is 

jyTrlF.Fh 
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where k,l = x, y, z. Calculating the trace in the numerator and similarly the trace in the denom- 
inator, considering unpolarized incoming leptons, one finds 



sin 2 t? • 



\G Mp (t)\ 2 + \G Ep {t)f 



Pyy = sin 2 M , , r , t ~' > — (? ' 12) 



P = P = 

1 xy yx u 



R 



r [1 It I G Ep (t ) 1 2 sin 2 + |G Mp (*)| 2 (l + cos 2 0)] ' 

gg P wi 2 -TiG MP ft)r 

r[l/r|G £p (t)| 2 sin 2 + |G Mp (*)| 2 (l + cos 2 <?)] : 
t(1 + cos 2 0)|G Mp (*)| 2 - sm 2 0|G gp (t)| 2 
r [1 /r | G Bp (i ) | 2 sin 2 + | G Mp (t ) | 2 ( 1 + cos 2 0)] 5 

sin20i?e[G^ p (i)G £p (i)] 



V^[l/T|G Bp (t)| 2 sin 2 + |G MP W| 2 (1 + cos 2 0)] 



p — p — n 



Now, if the expression (7.9) for the lepton tensor is used with the longitudinally polarized 
electron or positron, for components P xx , P yy , P zz , P xy , P yx , P xz , P zx one obtains identical 
expressions with (7.12), but the last two components are now nonzero as well 



smdIm[G* Mp {t)G Ep (t)} 
V^[l/r\G Ep (t)\ 2 sin 2 + |G Mp (t)| 2 (l + cos 2 0)] ' 



P - p - l 1V1 p \ y ^ 1 J \ y i , m 



Every of components Pki characterize a polarization of the proton p in the direction k, if 
antiproton p is polarized at the direction I and they all are calculated for 100% polarization of 
one of the initial leptons. 

As one can see from explicit formulae the vector polarization component P y and the tensor 
polarization components P X z=Pzx depend on the imaginary parts of electric and magnetic proton 
FFs. So, only our UhA models of the nucleon EM structure can give a sophisticate prediction of 
behaviors of P y and P xz =P zx - The eight and ten resonance U&A models represent compatible 
unification of pole and continua (in the language of the analyticity given by cuts on positive real 
axis) contributions, which give just imaginary parts of FFs different from zero starting always 
from the lowest branch point representing the lowest opened physical threshold. 

Now, exploiting behaviors of proton electric and magnetic FFs in the time-like region as 
predicted by our eight and ten-resonance UhA models one finds [160] behaviors of single and 
double spin polarization observables of the e + e~ — > pp process as they are presented in figures 
Fig. 7.3, 7.4. 



7.2 Polarization effects in e+e -> dd and experimental determination of time like 

deuteron form factors 



In the present paragraph the polarization observables in the reaction 

e-(fci) + e + (k)->d(pi) + d(p2), (7.14) 
where the momenta of the particles are indicated in brackets, are calculated [161] explicitly. 
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Figure 7.3. Prediction of the single polarizations observables by ten-resonance (left-hand) and eight- 
resonance (right-hand) U&A model 
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Figure 7.4. Prediction of the double polarizations observables by ten-resonance (left-hand) and eight- 
resonance (right-hand) U&A model 



With this aim one considers the case of unpolarized and longitudinally-polarized electron 
beam with production of vector- and tensor-polarized deuterons. The expressions of polariza- 
tion observables are given in terms of the deuteron EM FFs in time-like region, where they are 
complex functions of the momentum transfer squared. 

In the one-photon approximation, the differential cross section of the reaction (7. 14) in terms 
of the leptonic L^ v and hadronic W^ v tensors contraction (in the Born approximation one can 
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neglect the electron mass) is written as 

da _ a 2 f3 L^W^ 



(7.15) 



where a = 1/137 is the electromagnetic constant, (3 = y/l — 4m 2 ,/ q 2 is the deuteron velocity in 
the reaction center of mass system (CMS), nid is the deuteron mass and q is the four momentum 
of the virtual photon, q = k\ + k 2 = P1+P2 (note that the cross section is not averaged over the 
spins of the initial beams). 

The leptonic tensor (for the case of longitudinally polarized electron beam) is 

= -q 2 g pv + 2(k lt ,k 2 „ + k 2li k lu ) + 2i\e l _ ivap k la k 2p , (7.16) 

where A is the degree of the beam polarization (further we assume that the electron beam is 
completely polarized and consequently A = 1). 

The hadronic tensor can be expressed via the deuteron electromagnetic current J p , describing 
the transition 7* — > dd, as 

= J P J; . (7.17) 

As the deuteron is a spin-one nucleus, its electromagnetic current is completely described by 
three FFs. Assuming the P and C invariance of the hadron EM interaction this current can be 
decomposed (2.19) into three deuteron EM FFs Fi, F 2 and F 3 , by means of which analogically 
to the nucleon Sachs EM FFs Gcd, GmcI and GQd (2.20) are defined . 

When calculating the expression for the hadron tensor W^„ in terms of the deuteron electro- 
magnetic FFs, the spin-density matrices of the deuteron and antideuteron to be defined by means 
of deuteron polarization vectors £ and £* are 

£l/i£h, = - - Pl ^ V J + -^—E^paSpPla + ZQ^v (7.18) 

V m d J ^ m d 

c , * ( P2 l iP2u 



rn 



d 



if the deuteron polarization is measured and the antideuteron polarization is not measured. Here 
s M and Q^ v are the deuteron polarization four vector and quadrupole tensor, respectively. The 
four vector of the deuteron vector polarization s p and the deuteron quadrupole-polarization ten- 
sor Q pv satisfy the following conditions 

s 2 = -1, spi = 0, Q^ v = Q VfJt , Q w = 0, pi^Q^ = 0. 

Taking into account Eqs. (7.17) and (7.18), the hadronic tensor in the general case can be 
written as the sum of three terms 

= + w i>»(Y) + W^(T), (7.19) 

where Wu„(0) corresponds to the case of unpolarized deuteron and W^(V) (Wu V (T)) corre- 
sponds to the case of the vector (tensor) polarized deuteron. The explicit form of these terms 
is 
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- the unpolarized term W^ v {0): 

W (a 2 ) 
m d 

_ qpq v ~ _ Pig 

W^q 2 ) = 8m 2 T(l-T)\G Md \ 2 , 

W 2 (q 2 ) = 12m 2 (\G cd \ 2 -\r\G Md \ 2 + \r 2 \G Qd \ 2 ) (7.20) 
and t = — q 2 /Am 2 . 

- the term for vector polarization W^CV): 



WV(V) = — Si{q 2 )ep Vap s a q p + -^S 2 (q 2 )[pi P .e„ aap s a q a pi p - piuSp a apS a q^Pip\ + 

qcrPlp+PlvStiaapSaqaPlp], (7.21) 

m d 

S^q 2 ) = -3m 2 (T-l)\G M d\ 2 , 

S 2 (q 2 ) = 3m 2 d [\G Md \ 2 ~ 2Re(G Cd - \G Qd )G* Md ], 

S 3 (q 2 ) = &m 2 Im{G C d-^G Qd )G* Md . 
- the term for tensor polarization W pu (T): 

Wp V (T) = V 1 (q 2 )Qg fl „ + V 2 (q 2 )-%p lfl p 1 „ + (7.22) 

m d 

+ V 3 (q 2 )(p 1 pQ lJ +p lu Q P ,) + V 4 (q 2 )Q tn , +iV b {q 2 )(pipQ v -pivQp), 

where 

Qp = Qpvqv - Qp,qp = o 

f\ _i_ Qpq^ p. _ q^qa p. _ qpqa p. ?c 

^5 pv pv T" ^ ^ ^2 ^2 p^q^ 

<3 = Qapqaqp- (7.23) 

The tensor structure functions Vi{q 2 ) are combinations of deuteron FFs as follows 

V 1 (q 2 ) = -3\G Md \ 2 , 



V 2 (q 2 ) = 3 



\GM d \ 2 + ^ _ ^ Re(Gc d - 7^GQ d - rGMd)G*Q d , 

= -6r[|GMd| 2 + 2i? e G Q(i G^] , (7.24) 
Viiq 2 ) = -l2m 2 T(l-T)\G M d\ 2 , V 5 (q 2 ) = -12rIm(G Qd G* Md ). 
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Using the definitions of the cross-section (7.15), leptonic (7.16) and hadronic (7.19) tensors, 
one can easily derive the expression for the unpolarized differential cross section in terms of the 
structure functions Wi 2 (after averaging over the spins of the initial particles) 



da v 



T- 1 - 



{u-tf 



Am 2 d q 2 



where t = (fci — pi) 2 , u = (ki — P2) 1 ■ 

In the reaction CMS this expression can be written as 



da v 



a 2 /? 3 
V 



dQ 



D, D = t(1 + cos 2 0)\G Md \ 



3 • 2 

- sin 
2 



\G 



Cd\ 



^ 2 \G Qd \ 



(7.25) 



(7.26) 



where is the angle between the momenta of the deuteron (p) and the electron beam (fc). In- 
tegrating the expression (7.26) with respect to the deuteron angular variables one obtains the 
following formula for the total cross section of the reaction (7. 14) 



<Jtot(e + ( 



dd) = 



2«3 r 



TTO 2 ^ 

3q 2 



3\G Cd \ 2 + iT(\G Md \ 2 + ^T\G Qd \' 



(7.27) 



One can define also an angular asymmetry, R, with respect to the differential cross section 
measured at 6 = n/2, a a 



da v 



dn 



= ct (1 + Rcos 2 9) 1 



where R can be expressed as a function of the deuteron FFs 



R 



2t(\G 



M 



lr\G Q \ 2 ) 



3|G C | 2 



2t(\G m \ 2 + It\G q \ 



+ 3|G C 



(7.28) 



(7.29) 



This observable should be sensitive to the different underlying assumptions on deuteron FFs; 
therefore, a precise measurement of this quantity, which does not require polarized particles, 
would be very interesting. 

One can see from (7.26) that, as in the space-like region, the measurement of the angular 
distribution of the outgoing deuteron determines the modulus of the magnetic FF. The separation 
of the charge and quadrupole FFs requires the measurement of polarization observables to be 
more convenient derived in CMS. When considering the polarization of the final particle, we 
choose a reference system with the z axis along the momentum of this particle (in our case it 
is p). The y axis is normal to the reaction plane in the direction of k x p, x, y and z form a 
right-handed coordinate system. 

The cross section can be written, in the general case, as the sum of unpolarized and polarized 
terms, corresponding to the different polarization states and polarization directions of the incident 
and scattered particles 



da 

dn 



da un 

— [l + P y + X Px + X Pi 



z ~t~ Pzz^zz Pxz^xz ~f~ Pxxi^^xx ^yy) ^~ ^Pyz-Ryz] ■> 
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(7.30) 

where Pi (Pij), i, j — x,y,z are the components of the polarization vector (tensor) of the out- 
going deuteron, Rij, i, j = x,y,z the components of the quadrupole polarization tensor of the 

da un 

outgoing deuteron Q ul/ , in its rest system and is the differential cross-section for the un- 

dil 

polarized case. 

The degree of longitudinal polarization of the electron beam, A, is explicitly indicated, in 
order to stress the origin of the specific polarization observables. 

Now, let us consider the different polarization observables and give their expression in terms 
of the deuteron EM FFs. 

• The vector polarization of the outgoing deuteron, P y , which does not require polarization 
in the initial state is 



Py = -^sin(20)/m [(g w - T -G Qd ) G* Md 



(7.31) 



• The part of the differential cross section that depends on the tensor polarization can be 
written as follows 

dor _ da zz da xz da xx 

dn ~ dn Uzz + dn Uxz + dn [Hxx Uvv) ' 

da zz _ a 2 /3 3 3r 

~dn ~ Aq 2 T 
da xz a 2 ^ 
dn 4<7 
da xx a 2 /3 3 3r 



(1 + cos 2 9)\G Md \ 2 + 8 sin 2 6 i^-\G Qd \ 2 - Re{G cd .G* Qd , 



2 3T 3 / 2 sm(2d)Re(G Qd .G* Md ), 

" sin 2 6\G M d\ 2 , (7.32) 



dn Aq 2 

• Let us consider now the case of a longitudinally polarized electron beam. The other two 
components of the deuteron vector polarization (P x , P z ) require the initial particle polar- 
ization and are 

P x = -3^ sing Re (c cd - T -G Qd ) G* Md , P z - ^ cos 6\G Md \ 2 - (7.33) 

From angular momentum and helicity conservations it follows that the sign of the deuteron 
polarization component P z in the forward direction (9 = 0) must coincide with the sign of the 
electron beam polarization. This requirement is satisfied by Eq. (7.33). 

A possible nonzero phase difference between the deuteron FFs leads to another T-odd polar- 
ization observable proportional to the R yz component of the tensor polarization of the deuteron. 
The part of the differential cross section that depends on the correlation between the longitudinal 
polarization of the electron beam and the deuteron tensor polarization can be written as follows 

l£F = ^ 6t3/2 Sin eim ( G Md.G* Qd )R yz . (7.34) 
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The deuteron FFs in the time-like region are complex functions. In the case of unpolarized 
initial and final particles, the differential cross section depends only on the squared modulus 
\G M d\ 2 and on the combination G = \G C d\ 2 + ^i~ 2 \GQ d \ 2 . So, the measurement of the angular 
distribution allows one to determine |GMd| an d the quantity G, as in the elastic electron-deuteron 
scattering. 

Let us discuss, which information can be obtained by measuring the polarization observables 
derived above. Three relative phases exist for three deuteron EM FFs, which we note as follows: 
«i = a M -ct Q ,a 2 = a M -a c , and a 3 = a Q - a c , where a M = ArgG M d, «c = ArgG C d, 
and uq = ArgGgd- These phases are important characteristics of FFs in the time-like region 
since they result from the strong interaction between final particles. 

Let us consider the ratio of the polarizations P yz (we would like to note that it requires a 
longitudinally polarized electron beam) and P xz (when the electron beam is unpolarized). One 
finds 

Pxz 

R\ = ~— = — cos cot o;i . (7.35) 

Pyz 

So, the measurement of this ratio gives us information about the relative phase a\. The measure- 
ment of another ratio of polarizations, R 2 = P xz /Pxx gives us information about the quantity 

\G Qd \ 

R 2 = ^ = 8^7 cot 6> cos ai 1^4- (7.36) 

Pxx \(*Md\ 

This allows one to obtain the modulus of the charge FF, |Gcr<i|, from the quantity G, known from 
the measurement of the differential cross section. The measurement of a third ratio 

P y sina 2 -rsm ai r \G Qd \ 

R3 = -^- = -cos8 , r=-j— (7.37) 

P x cos a 2 — r cos a\ 3 |Gc<z| 

allows to determine the phase difference a 2 . And at last, if we measure the ratio of the polariza- 
tions P zz and P xx 



p 

R 



P xx sin 2 1 



1 , 2 a , o • 2 a \Gcd\\GQd\ , . 

1 + cos 9 + 8sm z 9 — — ^—(r - cosa 3 ) 

\GMd\ 



(7.38) 



we can obtain information about the third phase difference a 3 . Moreover, one can verify the 
relation: 

0:3 = a 2 — ai. 

Thus, the measurement of the polarization observables in the process (7.14) allows to deter- 
mine the deuteron EM FFs in the time-like region for the firs time. 



7.3 Alternative method of experimental determination of deuteron electromagnetic form 

factors 

In 4.3 we have demonstrated that alternative method of a measurement of the same physical 
quantity, in this case the proton electric FF GE P (t), can lead to very surprising results. To be 
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motivated by means of this example, here we propose [162] to experimentalists a new method 
of experimental measurement of all three deuteron EM FFs in the space-like region for the first 
time. 

Knowing experimentally the FFs (2.19) Fi(t), F 2 (t) and F 3 (t), one can specify experimental 
behaviors of G c (t), G M (t) and G Q {t) by (2.20). 

The FF F 2 (t) can be found by a measurement of the unpolarized differential cross-section in 
the laboratory system 



J Lab 
uu unp 



i 

647T 2 77j3 



l + (2£ e /m d )sin 2 (0/2). 



\M\ 2 unp 



(7.39) 



with 



\M\ 



unp 



e 4 m 2 cos 2 (6>/2) 
~Ef sin 4 (0/2) 



l + {2E e /m d )sm 2 (6/2) 



A(t)+B(t) tan 2 (6</2) 



, (7.40) 



and deuteron elastic structure functions (2.22). 

In order to find Fi(t), we suggest to measure elastic electron-deuteron scattering with only 
vector polarized deuterons, where the incoming (target) deuteron is polarized in the direction ft of 
the three momenta of outgoing deuteron and the vector polarization of the outgoing deuteron has 
the same direction. Then in the laboratory system (the rest frame of the incoming deuteron) the 
four vectors of electron and deuteron momenta {k\, p\, p 2 ) and the deuteron vector polarizations 
£i, £2 can be expressed in components as 



k 1 = (E e ,0,0,E e ); pi 



m\ 



{m d , 0,0,0); p 2 = (E 2 , \p 2 \n); 



6 = (0,n); & = — (\P2\,E 2 n); 

m d 



= 2m dv /^(l + ^); E 2 =m d (l + 2 V ), 



(7.41) 



where the outgoing deuteron vector polarization was obtained by Lorentz boost to be the kine- 
matic frame of the recoil deuteron. As a result one obtains the following relations 



= -1-2775 Kl-P2) 

(6 • h)(£,2 -Pi) = -2m d (E e + m d )ri; 

V 



IP2I 



(6 -Pl) = \P2\ 



Ki • fci)(6 • h) 



l + T) 



{E e +m d )[m d (l + 2r))-E e ]; 



(pi ■ fei) = E e m d ; {pi ■ p 2 ) = m d {\ + 2rj) 



(7.42) 



to be useful in further calculations. 

The absolute value squared for the amplitude of the polarized elastic electron-deuteron scat- 



Polarization phenomena in electromagnetic interactions of hadrons 



127 



tering is 



\M\ 2 



n 



spins ^ 



+ f 3 



(&*■</)(&■«) 



2m 2 

*l(6-£) + *3 



(&■?)(£■«) 



2m 2 



(P1+P2) 



G 2 



(7.43) 



then by using properties of the deuteron polarization vectors £1 , £2. (7-41) and (7.42), it can be 
arranged by the computer program FORM and simplified by Maxima to the form 



..,,9 1 e 4 m 2 , cos 2 (#/2) / . „ , , . 



i4(t)+B(t) tan 2 (0/2) + -^ 1 2 (t) 



(7.44) 



In the latter equation one can identify the absolute value squared of the unpolarized matrix 
element (7.40) and as a result one obtains 

\M\ 2 n = \\M\l p + + (2E e /m d ) f dn a (*/2))*?(t). (7.45) 

Multiplying the relation (7.45) by 
1 



647T 2 ™ 2 \ l + {2E e /m d ) sin 2 ((9/2) 
one comes to the equation 

a 2 cos 2 (#/2) 



</<> 



3 dil 8E 2 sin 4 (6»/2) 1 + (2E e /m d ) sin 2 (<9/2) 

from which the EM FF F\ (t) can be extracted. 

The third FF F3 (t) is given by a solution of the quadratic equation 



F 2 (t). 



(7.46) 



12 



v 2 {i + v 2 



+ 



Fi(t) + 

16 



(7.47) 



3^7(1 + v)Fi(t) + yry 2 (l + V )(Fi(t) - F 2 (t)) 



F 3 (t) + 



Fi(t) + -»y(Fi(t)-F 2 (t)) 



^Fl(t) + -r 1 2 (F 1 (t)-F 2 (t)) 2 -A(t) 



}- 







(with 77 = t/4m^) to be obtained by a substitution of (2.20) into the first relation of (2.22). 

Ones the data on i^i) an d ^H*) ^ known, by using (2.20) one can obtain the 

experimental behavior of Gc(t), Gm(*) and Gq(*). 
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8 Muon anomalous magnetic moment 

The muon is described by the Dirac equation and its magnetic moment is related to the spin by 
means of the expression 




(8.1) 



where the value of gyromagnetic ratio g is predicted (in the absence of the Pauli term) to be 
exactly 2. 

However, interactions existing in nature modify g to be exceeding the value 2 because of the 
emission and absorption of 

• virtual photons (electromagnetic effects), 

• intermediate vector and Higgs bosons (weak interaction effects) 

• vacuum polarization into virtual hadronic states (strong interaction effects). 

In order to describe this modification of g theoretically, the magnetic anomaly was introduced 
by the relation 




(8.2) 



where to every order some Feynman diagrams (see Figs. 8.1-8.3) correspond and 

a = 1/137.03599976(50) 

is the fine structure constant. 

The muon anomalous magnetic moment a M is very interesting object for theoretical investi- 
gations due to the following reasons 



Figure 8.1. The simplest Feynman diagram of an interaction of the muon with an external magnetic field. 
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Figure 8.2. The lowest-order hadronic vacuum-polarization contribution to the anomalous magnetic mo- 
ment of the muon. 




Figure 8.3. The third-order hadronic vacuum-polarization contributions to the anomalous magnetic moment 
of the muon. 



i) it is one of the best measured quantities (BNL E-821 experiment) in physics [163] 

a e * p = (116592040 ± 86) x 10" 11 (8.3) 

ii) its accurate theoretical evaluation provides an extremely clean test of "Electroweak theory" 
and may give hints on possible deviations from Standard Model (SM) 

Hi) moreover, in near future the measurement in BNL is expected to be performed yet with an 
improved accuracy 



A«7 = ±40 x 10~ n 



(8.4) 



i.e. it is aimed at obtaining a factor 2 in a precision above that of the last E-821 measure- 
ments. 
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ABC 



Figure 8.4. Third order hadronic light-by-light scattering contribution to a*° d (.4) and class of pseudoscalar 
meson square loop diagrams (B) and quark square loop diagrams (C) contributing to (.4). 



At the aimed level of the precision (8.4) a sensibility will already exist to contributions 

a (2,3)weak = ^ ± ^ x jq-U ( g 5) 

arising from single- and two-loop weak interaction diagrams. And so, if we compare theoretical 
evaluations of QED contributions up to 8th order 

a QED = (H6584705.7 ±2.9) x 10" 11 [164] 

the single- and two-loop weak contributions 

o ( ( 2,3)t U eofc = (151±4)xl0 _ii [165] 
a (2,3) W eak = ^53 ± 3) x lfJ -ll [166] 

4 2 ' 3)u,eafe = (152 ± 1) x 10- 11 [167] 
strong interaction contributions 



n had 
a fi — 


(7068 ± 172) x 10- 


-11 


[168] 


r .had 

% _ 


(7100 ± 116) x 10" 


-11 


[169] 


r .had 

u fi — 


: (7052 ± 76) x 10" 


11 


[170] 


ri had 

a fi — 


(7024 ± 152) x 10- 


-11 


[171] 


n had 

<v - 


: (7021 ± 76) x 10- 


11 


[172] 



it is straightforward to see that the largest uncertainty is in a^ ad . 

Error is comparable, or in the best case two times smaller than the weak interaction contribu- 
tions. 

So, in order to test the SM predictions for and to look for new physics in comparison with 
BNL E-821 experiment, one has still to improve an evaluation of a^ ad and its error. 

The most critical from all hadronic contributions are the light-by-light (LBL) contributions 
(see Fig. 8.4) to be approximated by a sum of meson pole terms (see Fig. 8.5). 

Therefore, we have recalculated [173] the third-order hadronic LBL contributions to the 
anomalous magnetic moment of the muon a^ ad from the pole terms of the scalar a, a and pseu- 
doscalar 7T°, rj, rj' mesons (M) in the framework of the linearized extended Nambu-Jona-Lasinio 
model 



£qqM = 9Aiq(x) [a(x) + itt(x)j 5 ] q(x) . 
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C 



Figure 8.5. Meson pole terms contribution to LBL diagram 



The reason for the latter are predictions of series of recent papers 

a LBL = (+52 ± 18) x 1CT 11 [174] 
a LBL = (+92 ± 32) x 10" 11 [175] 
a^ BL = (+83+12) x 10- 11 [176] 
a fl BL M = (+58 ± 10) x 1CT 11 [177] 

which differ in the magnitude. 

Moreover, in these papers only the pseudoscalar meson pole contributions were considered. 
We include the scalar meson (a, a ) pole contributions as well. 

Current methods in a description of the 7* — »■ M7* transition form factors are ChPT and the 
vector-meson-dominance (VMD) model. 

Here the corresponding transition form factors by the constituent quark triangle loops (see 
Fig. 8.5) with colorless and flavorless quarks with charge equal to the electron one are repre- 
sented. 

An application of a similar modified constituent quark triangle loop model for a prediction of 
the pion electromagnetic form factor behavior was carried out in [178] where also a comparison 
with the naive VMD model prediction was demonstrated. 

The mass of the quark in the triangle loop is taken to be 

m u = m d = m q = (280 ± 20) MeV 

determined [179] in the framework of the chiral quark model of the Nambu-Jona-Lasinio type by 
exploiting the experimental values of the pion decay constant, the p-meson decay into two-pions 
constant, the masses of pion and kaon and the mass difference of r\ and rj mesons. 
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The unknown strong coupling constants of ir°, i], rj and ao mesons with quarks are evaluated 
in a comparison of the corresponding theoretical two-photon widths with experimental ones. 

The cr-meson coupling constant is taken to be equal to 7r°-meson coupling constant as it 
follows from the corresponding Lagrangian. 

The cr-meson mass is taken to be m CT =(496 ± 47) MeV as an average of the values recently 
obtained experimentally from the decay D + — > 7r~7r + 7r + [180] and excited T decay [181] 
processes. 

As a result we present explicit formulas for a LBL (M) (M — 7r°, tj, rj , <r, a ) in terms of 
Feynman parametric integrals of 10-dimensional order, which subsequently are calculated by 
MIKOR method. 

Finally, one finds 

(81.83 ± 16.50) x 10^ n (8.6) 
(5.62 ± 1.25) x 10" 11 
(8.00 ± 1.74) x 10 -11 
(11.67 ±2.38) x 10" 11 
(0.62 ±0.24) x 10" 11 . 

So, the total contribution of meson poles in LBL is 



a L » BL (v) 
a LBL (ao) 



% (M) = (107.74 ± 16.81) x 10" , (8.7) 

where the resultant error is the addition in quadrature of all partial errors of (8.6). 

Together with the contributions of the pseudoscalar meson (n ± ,K ± ) square loops and con- 
stituent quark square loops (Fig. 8.4) it gives 



a LBL (total) = (111.20 ± 16.81) x 10" n . (8.8) 

The others 3-loop hadronic contributions derived from the hadronic vacuum polarizations 
(VP) were most recently evaluated by Krause [182] 

a (3)V-P = ( _ 101±6)xl() -ll 

Then the total 3-loop hadronic correction is 

a (3)had = q^bl {total) + af VP = (10.20 ± 17.28) x 10~ n (8.9) 

where the errors have been again added in quadratures. 

If we take into account the most recent evaluation [172] of the lowest-order hadronic vacuum- 
polarization contribution to the anomalous magnetic moment of the muon 



a (2)had = ( 7 q 21 ± yg) x 1Q - 
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together with the pure QED contribution up to 8th order and the single- and two-loop weak 
interaction contribution, finally one gets the SM theoretical prediction of the muon anomalous 
magnetic moment value to be 

af = (116591888.9 ± 78.1) x HT 11 . (8.10) 

Comparing this theoretical result with experimental, one finds 



- af = (151 ± 116) x 1CT 11 (8.11) 

which implies a reasonable consistency of the SM prediction for the anomalous magnetic mo- 
ment of the muon with experiment. 

However, one expects in near future a two times lowering of the error in BNL E-821 experi- 
ment and then there can appear still a room for a new physics beyond the SM. 

Therefore, one has to think on further improvement of the central value of the muon anoma- 
lous magnetic moment and especially of the lowering of its total theoretical error. 



8.1 Remarkable suppression of the e+e -> tt + -k contribution error into muon g — 2 

The first improvements we see still in the lowest-order hadronic vacuum-polarization diagram 
contributions (Fig. 8.2) to be dominant among all other hadronic contributions, which can be 
expressed by the integral 

1 f°° 

a( ? )had = ^J 4 2 V h (s)K^ S )d S ; (8.12) 
where <r h (s) stands for the total cross section cr(e+e~ — > had) and 

" [S) ~J x2 + (l-x) S /mf X - 

Moreover, the two-pion channel a!^ ad ' LO {e + e~ — > ir + n~), is dominant among all the lowest 
hadronic vacuum-polarization contributions. The data are used to evaluate the integral up to few 
GeV, above a perturbative calculation is possible. 

The expression (8.12) was always evaluated [183, 184] only by the direct integration of the 
0lo (ti + it~) data (with small exceptions, see [185]). 

In this paragraph we demonstrate [186] on a^ ad,io (e + e~ — > 7r + 7r~), that by using the 
USzA model of the pion EM structure (see 3.3) one can achieve remarkable suppression of the 
contribution error in comparison with a direct integration over existing experimental data. In 
order to compare our results with results of other authors obtained by integrating over the existing 
experimental data, we chose three different limits for the upper integration limit (3.24 GeV 2 , 
2.0449 GeV 2 and 0.8 GeV 2 ). 

Two approaches were used for the error evaluation. 
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The first one was a Monte Carlo method based on a random number generator with the 
assumption of the Gaussian distribution for the uncertainties of the published data points. For 
each point a new one was randomly generated using the Gaussian probability density function 
with the mean identical to the original point and a equal to the published error. Doing this for 
each point, new random data set was obtained. This data set was then fitted and the values of the 
parameters of the model pi as well as the value of a ^f d - LO (-k + it~) were extracted. Repeating 
the whole procedure 4000 times, we reached statistics high enough to allow us for a reliable 

error calculation. The mean a^ ad ' LO (ir + ir~) and the a were calculated from the 4000 values 
and since the mean is not, in general, identical with the optimal-fit value a® 'p T (ir + Tr~) we 
present asymmetric uncertainties 

<ad ,Lo (7r+7r - } = $£? 

where 

< had,LO ( + — \ had.LO I + — \ 

A = 17 + a n,OPT l 71 " 71 ) - a M ^ \JT 7T ) 

and 

r, had.LO I + — \ had.LO / + — \ 

B = a + a M yK it ) - a M;0 ' PT {tt^tt ) . 

In the second approach the program MINUIT was used to establish the uncertainties of the 
model parameters. Then, taking the numerical derivatives for ^a^ l ad ' LO (ir + n~), the uncer- 
tainty was propagated to a ^f d - LO (7r + 7r~) using the covariance matrix. In this method the errors 
are symmetric. 

In addition to the integration of the model, we also performed a direct integration of the data 
points based on the trapezoidal rule, so as to cross-check our compatibility with other authors. 
Our results and some results from other authors [183-185] are summarized in Table 8.1. 

The use of the U&A model dramatically reduces the error on a ^ ad ' LO (tt + tt^). This is not 
an arbitrary feature of the model but originates from model-independent information which is 
additional to the data in the integration region and which can be taken into account when the 
model is used. 

The most important sources contributing to error reduction are 



Table 8.1. Our results and results of other authors [183-185]. 

a» aci < LU (e + e- tt+tt') x 10 11 

Interval [GeV 2 ] Amj < t < 3.24 Amj < t < 2.0449 Ami <t<0.8 

XJhA Model, Meth.l 5132.36+g;^ 5128.22"^ 4870.24;°|° 

U&A Model, Meth.2 5132.37 ± 3.00 5128.25 ± 2.86 4870.44 ±2.64 

Integration over Data 5035.33;2l'i 5031.22^;^ 4756.77;2?.'55 

Davier 5040.00 ± 31.05 

Hagiwara et al. 5008.2 ± 28.70 

Yndurdin et al. 4715 ± 33.53 
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• Expected smoothness of the F w (t) at small scale At: The model provides a function be- 
having smoothly at small At. 

• Experimental data outside the integration region: The fit is done not only to the data inside, 
but also to the data outside the integration region. 

• Theoretical knowledge on F n (t): The model respects all known properties of the pion EM 
FF. 

Especially the first point plays an important role. The new precise data tend to lie above 
older, less precise data and, in some regions, the vertical spread of the data is very important, at 
the limit of inconsistency. If the calculation of the integral is based directly on data, then less 
precise data shift the mean value of the integral and enlarge the uncertainty. 

When the USzA model is used, the predicted behavior of (t) as given by the result of the 
fit is mostly determined by precisely measured points and is only little influenced by data with 
important uncertainties. This leads to more appropriate mean value and smaller errors. 

As a result, one arrives to the mean value of a '^ ad - LO (ir + ir~) which is higher than what is 
obtained by the direct data integration and to much reduced uncertainty. The shift in the mean 
value goes in the right direction and brings the theoretical value closer to the experimental one. 

The error estimates from the two used methods are not fully compatible, the first method gives 
smaller errors. This might be related to statistical fluctuations (1st method) and to approximations 
as numerical derivatives and linearization (2nd method). 

Here we have presented the calculation of a ^ ad > LO (7r+7r _ ) based on the XJhA model of the 
pion EM structure. This approach allows for important error reduction. It can be extended also to 
an evaluation of the contributions of other channels to a h ^ d ' LO with two-particles in final states. 
As the important multi-particle states can be mostly reduced to two-resonance final states, the 
method exploiting the UhA models can be also applied to them. 
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9 Sum rules for hadron photo-production on hadrons and photon 

Under the sum rules commonly one understands the expressions bringing into relation some 
physical quantities with other physical quantities. There exist various approaches for a derivation 
of a complex of sum rules as it is fully demonstrated e.g. in [187]. In this paper we are concerned 
of the sum rules for hadron photo-production on hadrons and photon. 

Historically there was only one attempt, by Kurt Gottfried [188], to derive sum rule for 
hadron photo-production, specially on proton target, relating the proton mean squared charge 
radius < r\ > and the proton magnetic moment /j, p = 1 + n p to the integral over the total 
hadron photo-production cross-section on proton, considering very high-energy electron-proton 
scattering and the non-relativistic quark model of hadrons. However, the corresponding integral 
is divergent and the Gottfried sum rule practically cannot be satisfied. 

Here we are interested in a derivation of hadron photo-production sum rules, bringing into 
relation static properties of pseudoscalar meson nonet, 1/2+ ground state octet baryons and 
quarks with the convergent integral over the total cross-sections of hadron photo-production on 
pseudoscalar mesons, on octet baryons and over 77 — > 2jets cross-section, respectively, ei- 
ther by exploiting analytic properties of "the retarded forward Compton scattering amplitude on 
hadron" [189] A(s\, q) in si-plane (such amplitude represents only a class of diagrams in which 
the initial state photon is first absorbed by a hadron line and then emitted by the scattered hadron) 
for meson [190] and baryon [191] targets, or by its explicit calculation in quark loops approxima- 
tion for the case of photon [192] target. The variable s\ is the cm. energy squared of the virtual 
Compton scattering process and q is the third component of the transversal part q± = (0, 0, q) 
of the virtual photon four-momentum q. 

This new approach gives sum rules with convergent integrals and specially the proton-neutron 
sum rule [105] was tested by using existing data and it is fulfilled with high precision. The latter 
fact gives us confidence also to all other derived sum rules by a similar approach. 



9.1 q 2 -dependent meson sum rules 

The q 2 -dependent meson sum rules are derived by investigating the analytic properties of the 
retarded Compton scattering amplitude A h (si,q) in s\- plane as presented in Fig. 9.1a, then 
defining the integral / over the path C (for more detail see [193]) in the si-plane 

/ = /^i^(i^(si,q)-<,(si,q)) (9.1) 

c 

from the gauge invariant light-cone projection p^p'(Aj tv (si, q) of the amplitude A h (s\, q) and 
once closing the contour C to upper half-plane, another one to lower half-plane (see Fig. 9.1b). 
As a result the following sum rule appears 

00 

-^[ImA h (s u cO- ImA h { Sl ,q)]. (9.2) 

r.h. 

The left-hand cut contributions expressed by an integral over the difference [ImA h {s\, q) — 
ImA h (si,q)] are assumed to be mutually annulated. 
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Now, one has to take into account the corresponding residuum of the intermediate state pole 
(see Fig. 9.1). 

As the electromagnetic structure of mesons is described by one charge form factor, the 
residuum takes the form 

iW M) = 27raF| f (-q 2 ), (9.3) 

where an averaging over the initial photon spin is performed. 

Then, substituting (9.3) into (9.2) and taking into account [190] 

,da e ~ h ^ e ~ x (s,q) da e ~ h '^ e ~ x ' '(s,q)\ _ 
I ¥q ¥q ) ~ 

OO 

f ^[ImAHsi.q) -Jm/*' («i,q)]. (9.4) 
J s 1 



47r 2 q 2 



with cPq = ndq 2 , one comes to the q 2 - dependent meson sum rule 

[^|,(-q 2 ) - F%{0)] - [F|(-q 2 ) - F|(0)] = 
9 /da eP ^ eX dn eP '^ eX '\ 

= ^> 2 (^P ^^-Y 



where the left-hand side was re-normalized in order to separate the pure strong interactions from 
electromagnetic ones. 



138 



Analyticity in a phenomenology of electro-weak structure of hadrons 



9.2 Universal static sum rules for total hadron photo-production cross-sections on 

mesons 



Now, employing the Weicsacker- Williams like relation [190] 

j( dcr e ~ p ^ e ~ x _ da e ~ p '^ e ~ x ' 

dq 2 d<f J\^^ 



,(do~ ^ da" " \ 



= -J —\- a tot (*l)-Otot («l)] 



for mesons, taking a derivative according to q 2 of both sides in q 2 -dependent meson sum rule for 
q 2 — > and using the laboratory reference frame by si = 2m B uj, one comes to the new universal 
meson sum rule relating meson mean square charge radii to the integral over a difference of the 
corresponding total photo-production cross-sections on mesons 

\{{r%)-{r 2 P ))= (9.7) 



3 



dcu 



r jP^rX, \ _ -(P'^X 



■ (J. 



Lot 



(w) - o~L ~* (w)] . 



in which just a mutual cancelation of the rise of the latter cross sections for u — > oo is achieved. 

According to the SU(3) classification of existing hadrons the following ground state pseu- 
doscalar meson nonet 7r~, ir°, tt + , K~ , K°, K°, K + , rj, rf exists. However, in consequence of 
CPT invariance the meson electromagnetic form factors Fp(— q 2 ) hold the following relation 

F P (-q 2 ) = -Fp(-q 2 ), (9.8) 

where P means antiparticle. 

Since n , r\ and rf are true neutral particles, their electromagnetic form factors are owing to 
the (9.8) zero in the whole region of a definition and therefore we exclude them from further 
considerations. 

If one considers couples of particle-antiparticle like ir^ , K ± and K a , K°, the left hand side 
of (9.5) is owing to the relation (9.8) equal zero and we exclude couples ^ , K ± and K°, K° 
from further considerations as well. 

Considering a couple of the iso-doublet of kaons K + , K° and K~ , K°, the following Cabibbo- 
Radicati like sum rules [194] for kaons can be written 

\it 2 a(rj c+ ) = (9.9) 



6 

f 



duo 



Hot M - °tot M 



w 2 a(-l)(r 2 K -) = (9.10) 



f°° du> r lK -^ x , \ -yK°->-X, \ 
= / 77 ®tot M-^/ot (W) 
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in which the relation {r 2 K + ) 



— {r 2 K -) for kaon mean squared charge radii, following di- 
rectly from (9.8), holds and divergence of the integrals, due to an increase of the total cross- 
sections a]^ ~ > ' X (uj) for large values of u>, is taken off by the increase of total cross-sections 



cr. 



(w) and a 



(uj), respectively. If besides the latter, also the relations 



°tot M = °Lt M 



following from C invariance of the electromagnetic interactions, are taken into account, one 
can see the sum rule (9.10), as well as all other possible sum rules obtained by combinations 
K+K°, R-R", to be contained already in (9.9). 

The last possibility is a consideration of a couple of mesons taken from the isomultiplet of 
pions and the isomultiplet of kaons leading to the following sum rules 



1 



(9.11) 



f 



dbj 
u 



~ f K ± ^X 



(9.12) 



f°° [ 77r ±^X/ \ ~fK a ^rX t 



Now taking the experimental values [12] 

(±l)(r*±) = +0.4516 ± 0.0108 [fm 2 ] 

{±l)(r 2 K± ) = +0.3136 ± 0.0347 [fm 2 ] 
one comes to the conclusion that in average 



77r ± ->X 



H-a7 f ^M]>0 



-~ t K-^X 



(w) - cr; 



7 K°^X 



(w)] > 0, 



(9.13) 



from where the following chain of inequalities for finite values of uj in average follow 

irK^M > ^» > *Zf ">» > 0. (9.14) 



Subtracting up (9.9) or (9.10) from the relation (9.12), the sum rule (9.11) is obtained, what 
demonstrates a mutual consistency of all considered sum rules. 
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9.3 q 2 -dependent octet baryon sum rules 

Similarly to the meson case, as the electromagnetic structure of octet baryons is described by 
Dirac and Pauli form factors, the residuum takes the form 

Res B = 2na(F 2 B + A-F% B ), (9.15) 

where an averaging over the initial baryon and photon spins is performed. Finally one obtains 
the q 2 -dependent baryon sum rule in the form [191] 

[Fl B ,{^) - F 2 B ,(0)] [^(-q 2 ) - F 2 B (0)] + 

, n2 r ^(-q 2 ) f 2 b {~<?) 

+ q [ 4m|, 4m| J " 

= ^< q > (— ^ (9 - 16) 

where the left-hand side was again re-normalized in order to separate the pure strong interactions 
from electromagnetic ones. 



9.4 Universal static sum rules for total hadron photo-production cross-sections on 

baryons 

Employing in (9.16) the Weicsacker- Williams relation for baryons, taking a derivative according 
to q 2 of both sides in (^-dependent baryon sum rule for q 2 — > and using the laboratory 
reference frame by si = 2m B u) , one comes to the new universal static baryon sum rule 

I[F, B (0)(^>-^,(0)« B ,)]-[^-^-] = 

OO 

relating Dirac baryon mean square radii {r\ B ) and baryon anomalous magnetic moments n B 
to the convergent integral, in which a mutual cancelation of the rise of the corresponding total 
cross-sections for u> — > oo is achieved. 

According to the SU(3) classification of existing hadrons, there are known the following 
members of the ground state 1/2+ baryon octet (p, n, A , S + , S°, S~, S°, S~). As a result, 
by using the universal expression (9.17) one can write down 28 different sum rules for total 
cross-sections of hadron photo-production on ground state 1/2+ octet baryons. 

In order to evaluate their left hand sides and to draw out some phenomenological conse- 
quences, one needs the reliable values of Dirac baryon mean square radii (r 2 B ) and baryon 
anomalous magnetic moments k b . 

The latter are known (besides TP, which is found from the well known relation k s + + 
k s -=2k s o) experimentally [12]. However, to calculate (r 2 B ) by means of the difference of 
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the baryon electric mean square radius (r 2 EB ) and Foldy term, well known for all ground state 
octet baryons from the experimental information on the magnetic moments given by Review of 
Particle Physics 

(r! B ) = (r 2 EB )-^, (9-18) 

B 

we are in need of the reliable values of (r EB ). 

They are known experimentally only for the proton, neutron and X~ -hyperon. 

Fortunately there are recent results [195] to fourth order in relativistic baryon chiral per- 
turbation theory (giving predictions for the E~ charge radius and the A-E° transition moment 
in excellent agreement with the available experimental information), which solve our problem 
completely. 

Calculating the left-hand sides of all sum rules one finds 

oo 

OO 

_2_J ^[^ (u) _^ (w)]=2 ^ ^ ^>)>*&°->>) 

(9.20) 



OO 

i / ^Kf^M--7 f^M]= 4.2654m & , => ">» > ">» 

(9.21) 



OO 

(9.22) 



OO 

(9.23) 
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i/^ - <T^M] - 1.6673m6, =► > <^°^M 

(9.24) 

OO 

^ / ^[<TlV X n-^ + = -0.4158mfe, =► *2T» < ">» 



OO 

2 fdLJr , 
~2~~ / I* 7 '' 

-K z a J uj 1 



(9.25) 
(9.26) 



OO 



7r X (-) - ^"^M] - 3.8496m6, => > *£T~^M 

(9.27) 



_2_ | ^ [<72r * (w) _ ^ = 1.7259m6, > ^» 

(9.28) 



OO 

^ / V Kr*H - 'ZT^M] = 3.3180m6, => > *7«T 

(9.29) 



OO 

_2_ J ^ [ar(w) „ ff /,x H] = _ 03260m6; ^ - air x [w) < ^a^x (w) 
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OO 

_2_ J ^[^.^^.^ 



(9.30) 



-7T1- >X / \ , - 

°tot M < a 



(9.31) 



OO 

_2_ | ^[^M-^f^M] =-0.3747m6, 



(w) < <t (w) 
(9.32) 



ir^a J u) 



dcu 



K< ^M] = 1.8082m6, 



(9.33) 



OO 



(9.34) 



'tot 



(uj)] = 1.2766mfe, 



(9.35) 



OO 

7r z a J UJ 



(9.36) 



OO 

_2_ | ^[ CT7 A°^ M _ (7S o ^ M ] = _0.0006m6, 



- 7 A°-s-X 
T tot 



M « Cr/ ot (w) 
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(9.37) 



2 fdu ^ x {u) _ a g-^x M]= 2 . 1823mb; ^ a^ x ( U )>a^ x {u) 

(9.38) 

OO 

4- /^Kf^H-^ X H]= 0.0586mb, => ^ X H>^ X H 
7r z a J UJ 

(9.39) 

OO 

(9.40) 



OO 

Tr^a 7 77^ tot (w)J = 2.1417mb, <7 t 7 ot (w) > < t (w) 

(9.41) 



OO 

(9.42) 



OO 

^/^[^°^M-^ x M]=o.ii68mb, => *2^>)>*2?->>) 

(9.43) 



OO 
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k k k k 

Figure 9.2. Feynman diagram of cy — > e^f scattering with LBL mechanism to be realized by quark-loops 

(9.44) 



OO 

(9.45) 



OO 

' ' /W [a7 o f- X M-a7 o f^H]=-0.5316mb, a£~ (to) < o& ^ X (to), 

(9.46) 

from where one gets the following chain of inequalities 

0-*oi M > °tot M > °tot M ~ °tot M > 

for total cross-sections of hadron photo-production on ground state 1/2+ octet baryons to be 
valid in average for finite values of to. 



it 2 a J to 



9.5 Sum rule for photon target 

Let us consider the two photon exchange electron-photon zero angle scattering amplitude in the 
process 

e(p, A) + 7 (fc, e) -> e(p, A) + 7 (fc, e), (9.47) 

in two-loop (a 3 ) approximation as presented in Fig. 9.2, with p 2 = m 2 e , k 2 = and assuming 
that s = 2p.k ^> m 2 . 
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Averaging over the initial electron and photon spin states (initial and final spin states are 
supposed to coincide) one can write down the amplitude of the process (9.47) in the following 
form 



A^(s,t = 0) = s^j ^ds^All^ 1 ^-, (9.48) 
where the light-cone projection of the light-by-light (LBL) scattering tensor takes the form 



5*i S2 S3 

~D[ D~2 D~3 



with 



Si 


_ (l/4)Trp{q- + m q )p(q- 
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(9.49) 

,(9.50) 
,(9.51) 
,(9.52) 

where g_ means the quark four-momentum in the quark loop of the process 77 — > 77, N c is the 
number of colors in QCD and Q q is the charge of the quark q in electron charge units. 

Now, taking a derivative of the relation (9.48) according to d 2 q and investigating the analytic 
properties of the obtained expression in the si -plane one gets the configuration as presented in 
Fig. 9.3, where also the path C of the integral expression 

I= f d ^J^ML (9 . 53) 

Jc d 2 q 

is drawn. When the integration contour is closed to the right (on s-channel cut) and to the left 
(on the u-channel cut) one comes to the relation 

r ^, d -^gM> ]ult= r ^ ^m ^, 9S4 > 

7-4m 2 -q 2 « Q J4m 2 +q 2 " Q 
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where the right s-cannel discontinuity by means of Eq.(9.48) is related, due to optical theorem in 
a differential form 

As dR - 2S ^q~' (9 ' 55) 

to the Q 2 =q 2 =— q 2 dependent differential cross-section of qq pair creation by electron on photon, 
to be well known in the framework of QED [196] for l + l~ pair creation 

«s» = (q3 - m » )J+1 ' 

4 ln[ v /q 2 /(4m 2 ) + Jl + q 2 /(4m 2 )]. 

q 2 (q 2 +4m 2 ) 

But the right hand cut concerns of two real quark production for si > 4m 2 , which is associated 
wit 2 jets production. 

The left-hand cut contribution has the same form as in QED case with constituent quark 
masses and as a result one obtains 

Finally, for the case of small q 2 and applying the Weizsacker- Williams relation, one comes to 
the sum rule for photon target [192] as follows 

TE% = —2 I —-ir 2 ^). (9.58) 

3 ' HI 2 -KOe Z J Si 
q Am^ 

The quantity a]J t ~^ 2: ' ets (s±) is assumed to degrease for large values of s\. It corresponds to 
the events in 77 collisions with creation of two jets, which are not separated by rapidity gaps and 
for which until the present days there is no experimental information. The latter complicates a 
verification of the obtained sum rule for photon target. 
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10 Conclusions 

In this review we have formulated the main principles of a construction of the universal Unitary 
and Analytic model of the hadron electro-weak structure. Before the problem of inconsistency of 
the asymptotic behavior of VMD model with the asymptotic behavior of form factors of baryons 
and nuclei is solved, leading to a formulation of the asymptotic conditions and their general 
solutions giving explicitly all three possible forms of VMD form factor representations to be 
automatically normalized and possessing the right asymptotic behavior as required by the quark 
model of strongly interacting particles. Also a general approach for determination of the lowest 
normal and anomalous singularities of form factors from the corresponding Feynman diagrams 
is formulated. Finally, a number of concrete results to be obtained by making use of the analytic 
properties of the dynamical physical quantities, like electro-weak form factors of hadrons or am- 
plitudes of various electromagnetic processes of hadrons is presented. Some of them are already 
successfully verified in practice, others are expecting for experimental confirmation. All this 
together demonstrates the analyticity to be an unavoidable powerful tool of the present particle 
physics phenomenology. 
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